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CHAPTER TWO 


Static Aeroelasticity — Unswept wing structural loads and performance 


2.1 Background 


All structures deform when external loads are 
applied although the deflections may be barely 
discernible. In most cases, the external and 
internal loads do not depend on the structural 
deformation. From an analysis perspective this 
means we can compute the internal loads and 
the external deflections independently. 


\ Structural 
\ flexibility 


Figure 2.1.1 — Static aeroelasticity encompasses 
problems involving the intersection between 
steady-state aerodynamic and structural 
deformation interactions. 


These structural analysis problems are called 
statically determinate and include structural 
stability problems such as column buckling. 
However, if the loads and structural deflection interact the structural analysis problem becomes very 
different, both physically and computationally, because the problem is statically indeterminate. Both 
loads and deflections must be determined simultaneously. This load/deflection interaction is 
represented graphically by the Venn diagram in Figure 2.1.1 in which the overlapping orange area 
represents the statically indeterminate problem area. 


The purpose of this chapter is to illustrate the effect of aerodynamic load/structural deflection 
interactions on aircraft operation and performance. These problems are statically indeterminate and 
also involve considerations of structural stability. We will use simple models to reveal the causes of 
the aeroelastic phenomena and suggest cures for these problems. Chapter 2 contains homework 
problems to reinforce understanding these concepts. 


2.1.2 Scope and purpose 


As indicated in Figure 2.1.2, this 
chapter begins with a discussion 


of aeroelastic models and the 
introduction to special 
terminology required to define 
the features of these models. This 
includes a brief discussion of 
structural analysis matrix 
methods and concepts such as the 
shear center, aerodynamic 
coefficients and aerodynamic 
center of pressure. 


A primary goal to Chapter 2 is the 
development of a typical section 


aeroelastic model to illustrate 
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System modeling 
Matrix structural analysis 


Terminology - 


shear center, aerodynamic coefficients, 
lift curve slope... 


The typical section model 


Lift effectiveness Divergence — stability and system 
eigenvalues 
System modification — 
active control 


Figure 2.1.2 - Chapter 2 topical outline 


Control effectiveness 
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interactions between the structural internal forces and the external aerodynamic forces. The typical 
section model development is followed by the development of a multi-degree of freedom model to 
illustrate the mathematical criteria for static stability of an aeroelastic system. This model is then 
used to introduce the concepts of lift effectiveness and structural stability, namely the aeroelastic 
divergence phenomenon. The divergence discussion includes the eigenvalue problem that 
characterizes the stability of multi-degree-of-freedom systems. In addition to illustrating aeroelastic 
effects on wing load distribution and internal stresses, control effectiveness and aileron reversal are 
described using the simple typical section model. 


Finally, the modification of static aeroelastic response using ailerons as part of a feedback control 
system to adapt and adjust the lift and moments applied to flexible wings is also illustrated. This 
discussion will show how it is possible to modify the static stability and response of an aeroelastic 
system with feedback control. Although feedback control is seldom, if ever, used for static 
aeroelasticity this discussion shows the importance of controlling system interactions to aeroelastic 
response. 


2.1.2 The structures enterprise and its relation to aeroelasticity 


Every aircraft company has a large engineering division with a name such as “Structures 
Technology.” The purpose of the structures organization is to create an airplane flight structure with 
structural integrity. This organization also has the responsibility for determining and fulfilling 
structural design objectives and structural certification of production aircraft. In addition, the 
organization conducts research and develops or identifies new materials, techniques and information 
that will lead to new aircraft or improvements in existing aircraft. 


The structures group has primary responsibility for loads prediction, component strength analysis 
and structural component stability prediction. There is strong representation within this group of 
people with expertise in structural mechanics, metallurgy, aerodynamics and academic disciplines 
such as civil engineering, mechanical engineering, chemical engineering, and engineering 
mechanics, as well as the essential aeronautical engineering representation. 


The structural design process begins with very general, sometimes “fuzzy,” customer requirements 
that lead to clearly stated engineering design criteria with numbers attached. A summary of these 
general design criteria is shown in Figure 2.1.3. 


Beginning at the top of the “wheel” we Design loads 
have design loads. These loads include Environment and 
. . discrete events Static strength 
airframe loads encountered during 
landing and take-off, launch and Maintainability ; 
. . + Repairability Structures Stiffness and flutter 
deployment as well as in-flight loads + Inspectability Design 
i i gee Durabilit 
and other operational loadings. There Producibility Criteria ara 
are thousands of such “load sets.” * Corrosion 
: sg + Impact 
Once these load sets are identified, Crashworthiness 
: ‘ : : Fail safety 
there are at least nine design criteria 
‘ Damage tolerance 
that must be taken into account. On the and safe life 
wheel in Figure 2.1.3 stiffness and 
flutter are one important set of criteria Figure 2.1.3 - Structural design requirements 
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that must be addressed. 


The traditional airframe design and development process can be viewed as six interconnected blocks, 
shown in Figure 2.1.4. During Block 1 the external shape is chosen with system performance 
objectives in mind (e.g. range, lift and drag). Initial estimates of aircraft component weights use 
empirical data gathered from past experience. On the other hand, if the designs considered at this 
early stage have radical new forms, these estimates may be in error; but these errors will only be 
discovered later. 


Preliminary ° 
‘ : Production & follow-on 
— 
alias <———".. Detail design development 
Change decisions Errors discovered 
1-External 4-Stress 


shape, internal 2 - Choice of 3 - Detail analysis, 5- Full scale pc Sia ck 
ee Sr ap rf crawings ane and production 
° structure ore : i ’ 
j : writing and analysis, component ificati 
components, configuration posta hota certification aa modifications 
loads (skeleton) 


development 


Small group of specialists Large group of Preparation for certification by highly 


with broad experience. designers. Main qualified group of analysts. Adequate 
mathematical modeling. 


Main decisions based on decisions based on the : 
experience and intuition simplest mathematical PONG Male ROMemnes eGISS MOM: 
models. 


Figure 2.1.4 -— The structural design and development process requires testing, analysis and feedback 


Time is extremely important; the analysis scope and complexity during Block 1 is limited to very 
simple models despite the strong correlation between internal structural arrangement, materials 
selection, weight, manufacturability and cost. The ultimate success of the design effort depends 
heavily upon the training and real-world experience of the structural design team assigned to the 
Block 1 effort. Details that might escape the eye of a novice will be readily apparent to an 
experienced person. During Block 1, errors in judgment, caused primarily by lack of information, 
not lack of skill, are introduced, discovered in ensuing Blocks, and finally repaired. 


Block 2 is concerned with selecting the internal structural layout. This delay in working out the 
details of the structure may seem strange since a great deal of effort is expended during Block 1 
getting the external geometry (called the 
“outer mold lines”) right. For a wing, the 
structural effort includes the choice of the 
number of spars and ribs and the location of 
critical components. The activities in Block 2 
are characterized by low level structural 
analysis using simple beam models such as 
that shown in Figure 2.1.5 or other simple 
idealizations for element sizing and response 
estimation.’ Seldom will the efforts in Block 2 
revisit decisions made in Block 1, particularly 


those involving planform shapes. Figure 2.1.5 —Beam/stick models are used for 
preliminary flutter analysis (Reference 2.1) 
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Block 3 includes activities involving 
detailed drawings of frames, ribs, spars 
and attachments, as well as 
documentation. Here, extensive details 
emerge and if the manpower is available, 
higher fidelity models of critical 
components are created. 


Block 4, the certification activity, 
requires extensive stress analysis. Stress 
analysis is labor intensive and involves 
large, analytical models developed from 
the drawings generated in Block 3. 


Usually these complex models involve Figure 2.1.6 — Fighter aircraft finite element model 
showing interconnected aerodynamic and structural panels 


discretized finite element models of the 
type shown in Figure 2.1.6. Next, Block 
5 includes full scale and component tests including fatigue and static loads tests. Block 6 is the 
beginning of production and delivery to the customer. 


Note that Blocks 2 through 6 involve the opportunity and need to re-visit decisions made in earlier 
Blocks. Analytical modeling furnishes additional information that allows discovery and knowledge 
generation as the effort progresses. On the other hand, you cannot analyze what you have not 
defined, so surprises may visit the development team often. 


2.1.3. The evolution of aircraft wing structures-form follows function 


The structural mathematical models used in this book are simplistic and look very little like real 
aircraft structures such as the one shown in Figure 2.1.7. On the other hand, simplistic models led 
aeroelasticians and designers out of darkness into the light of understanding the inner workings of 
aeroelasticity and provided critical design guidance in the early years of aircraft development. To 
appreciate how analytical models are developed, we need some historical perspective on aircraft 
structural design evolution. A more 
complete perspective is contained 
in several books (References 2, 3, 
4, 5 and 6) and references listed at 
the end of this chapter . 


The geometry of modern high 
speed/high performance wing 
structures such as the F-22 are 
anything but simple, consisting of 
complex geometric arrays with 
large numbers’ of different 
structural elements. On the other 


Figure 2.1.7-The F-22 structure is a complex geometric layout Hany CO pate Uae a OGeEn TOW 


of basic structural elements created to transfer and resist speed, long endurance, high 
external and internal loads. altitude Global Hawk with its 


FIGHT. 
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slender, high aspect wings to the ancient Hawker Hurricane shown below; notice how little (with the 
exception of materials and manufacturing) the basic structural concepts have changed in 60 or 70 
years. 


Figure 2.1.8-The Global Hawk unmanned air vehicle began in 1995. The wing material is advanced 
composite carbon fiber. Sir Sydney Camm’s Hawker Hurricane fighter first flew in 1935. The original 
aircraft had fabric covered wings. Metal stressed skin wings were introduced in 1939. Both vehicles use a 

standard stressed skin, multi-spar wing layout. 


Wing structural design history is divided into two different eras; 1) the biplane/early monoplane era 
in which thin lifting surfaces were constructed of frame-like structures to which was attached a non- 
load-bearing cloth (or other material) material to give the surface its aerodynamic shape; and, 2) the 
reinforced, thick-wing, semi-monocoque wing era in which the entire surface, including the skin, 
carries structural loads and fuel and other components are housed in the wing. No matter what era 
we look at, the structures themselves are constructed of many different types of components 
arranged in geometries that provide aerodynamic shape and provide structural integrity. 


Figure 2.1.9-“Our invention relates to that class of flying-machines in which the weight is sustained by the 
reactions resulting when one or more aeroplanes are moved through the air edgewise at a small angle of 
incidence, either by the application of mechanical power or by the utilization of the force of gravity.” 
Wright Brothers Patent Application, March 1903. The 1903 Flyer airframe weighed 405 Ib. — the engine 
weighed about 200 lb. complete with radiator, water, fuel and accessories. Take-off weight was about 745 Ib. 


The first Wright Brothers airplane illustrates the type of biplane construction common to most 
aircraft during the first two decades of powered flight. While the Wrights added considerable 
technology to this aircraft, structural design was a settled subject, mostly due to previous efforts by 
Octave Chanute. 
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Octave Chanute was a distinguished railroad bridge engineer and self-made aeronautical hang-glider 
constructor who developed an interest in aviation at age sixty. He pioneered the biplane concept in 
the 1890’s. His hang-glider flights in the Indiana sand dunes near Lake Michigan convinced him 
that biplanes were the most practical solution for combining strong practical structures and 
aerodynamics. Given his background, it is no coincidence that these biplanes resembled railroad 
bridge truss construction. The cloth skin attached to the structure generated aerodynamic loads but 
was incapable of supporting loads (and not intended to). Instead, the aerodynamic loads were 
transferred to wing ribs and spanwise structural members near the leading and trailing edges and 
then conveyed to the central part of the glider to support the pilot. 


While Chanute and the Wright Brothers had solid scientific and engineering reasons for their biplane 
aerodynamic and structural designs, most early aviation experimenters neglected theory in favor of 
intuition and the artisan’s experience. In a brief description of early design procedures in 1916, 
Arthur W. Judge wrote: "In the earlier type of aeroplane body it was the usual practice to obtain the 
sizes of the different members by trial and error methods, or to make chance shots at the dimensions, 
and to trust to luck whether the resulting body has any margin of safety or not." Another article 
observed that "Faulty construction is the most fertile source of accidents, and always will be until 
constructors put first-class engineering knowledge into their work.”’ With few exceptions, 
analytical modeling of structures was not regarded as a primary study topic for most aeronautical 
engineers. 


As a result, in the beginning, aircraft structures were little more than cloth covered trusses and 
frames. It was the French, with their interest in monoplanes and fast, light-weight airplanes who 
changed this, at least for fuselage construction. In 1911, Louis Bechereau expropriated an idea from 
Eugene Ruchonnet, a Swiss engineer who had worked as a shop foreman for the French Antoinette 
Airplane Company. 


Ruchonnet focused his attention on the fuselage, which up to that time was constructed as a truss 
frame with cloth covering, much like the wing. Rather than simply attaching fabric to a truss 
structure, Ruchonnet’s idea was to use the fuselage skin carry the structural load and eliminate the 
heavy truss. As a result, he built a laminated wood fuselage shell to create a streamlined shape. 
These wood layers were glued together with their grains running in different directions to strengthen 
the skin. Ruchonnet called the new technique “monocoque” or single-shell construction. 


Bechereau used the monocoque fuselage idea, together with a frame/cloth skin monoplane design to 
build the Deperdussin streamlined, externally braced, mid-wing monoplane. This airplane won the 
Gordon Bennett Cup in 1913 with a speed of 125 mile per hour — only a decade after the Wright’s 
first flight! 


Later, as airplane size increased, the unreinforced, monocoque shell fuselage design concept was not 
suitable because the highly loaded compressive skin buckled under high loading. A variation of the 
monocoque concept, a reinforced or “semi-monocoque” structure was developed to create a light- 


* 


A. W. Judge, The Design of Aeroplanes, Isaac Pitman and Sons, Ltd., 1917, p. 156. This 242 page 
textbook by Judge is one of the early textbooks on airplane design. 
Aeronautical Journal, July 1911, p. 125 
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structure is still used today, both for fuselage 


Figure 2.1.10-Ruchonnet’s monocoque and semi-monocoque shell structures introduced yet another 
structural design concept into structural analysis and design. Deperdussin Racer “Le Monocoque” with 
monocoque shell structure fuselage, 1912. Compare this design to the Wright Brothers airplane only nine 
years earlier. 

Thin-wing monoplane structures were not as strong or as stiff as biplane arrangements so they 


required external wire bracing (see Figure 2.1.11). To create greater strength, wings began to use a 
substantial load-carrying member, the wing spar, located near the front of the wing and extending 
from the fuselage to the wing tip. The wing spars were usually conventional beam designs lifted 
from civil engineering structures practice. Originally the spar was simply a heavy beam member, 
but developed into a hollow box design that also resisted wing torsional loads and rotation. 


By the end of World War I, bi-plane wing design had matured to the extent that there were several 


variations of a common spar/rib configuration 
structure that supported a thin wing biplane or 
monoplane. This design used two spars with 
numerous ribs and spanwise stringers to help 
shape the wing. 


Early monoplanes still used thin wings that 
required external wire bracing to keep the wings 
from folding in flight. In Germany, Hugo Junkers 
recognized that streamlining the airplane shape 
was the key to increasing speeds with limited 
engine power. Junkers began by studying 
thermodynamics and became independently 
wealthy after inventing an automatic water heater. 
Junkers was also Professor of Thermodynamics at 
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Figure 2.1.11-The externally wire-braced 
Bleriot XI flew between England and France in 
1909. 
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the Technical University of Aachen (Germany). 
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Figure 2.1.12-Arthur Judge discussed the latest (circa 1918) wing structural designs with two different types 
of ribs. Most wings were biplanes. Immediately after World War I aircraft structures had advanced to the 
point where there were at least six different types of wooden spars, including one with laminated shear 
webs. The spar in E is intended to be a hollow leading edge and was used by Farman and Nieuport designs. 


Junkers found that thick airfoil sections performed better than thin sections in use at the time. He 
developed an internally braced cantilever monoplane design, made possible by thicker wings. Fixing 
the spar securely to the fuselage to create a cantilever beam arrangement allowed designers to 
eliminate the external bracing wires. Junkers’ monoplane used combinations of steel beams 
interconnected by load-bearing skin to create a semi-monocoque wing structure. Junkers founded a 
successful company and later produced the Junkers Ju 52-3m tri-motor which for a time was the 
favored transport aircraft in Europe. A trademark of this design was its corrugated metal stressed- 
skin. 


While cantilevered wings still used ribs and fabric to create the wing aerodynamic shape, designers 
began to use thin wooden veneers for skin. These wooden-skinned wings became known as 
stressed-skin wings because the veneers added strength while they shaped the wing surface itself. 
The structural material was nearly always wood, but Judge notes that “...steel and other metals (are) 
replacing these nonhomogeneous materials (wood) and successful all-metal wings and indeed 
complete machines have been built.” (Judge, p.68). 


In addition to Junkers, two others contributed to the development of revolutionary aircraft structures 
and airplanes with these structures before and after World War I. Like Junkers, these two were 
German. Claude Dornier and Adolf Rohrbach began their careers as structural engineers for the 
Graf Zeppelin Company. Both had extensive experience designing the internal metal framework for 
the large airships developed before World War I. Rohrbach’s academic training was in shipbuilding. 
Both eventually went on to create their own aircraft companies. 


Dornier worked for Count Zeppelin who, in 1914, asked him to develop flying boats. Dornier flying 
boats were famous and helped pioneer overseas routes for airlines. His structural designs were 
characterized by smooth wing skin that carried loads. 
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Figure 2.1.13-Three similar but different types of early stressed skin construction, left to right: Dornier, 
Junkers, Rohrbach. Dornier’s design is the result of his experience building Zeppelin airship structures. 
Junkers’ concept uses extensive trusswork. Rohrbach’s concept is cellular. The basis for Rohrbach's wing 
design was a central box section girder of thick duralumin sheet stiffened by fore and aft bulkheads 


Rohrbach worked for Dornier at the Zeppelin flying boat plant. After World War I he founded a 
company in Denmark to get around punitive features of the Versailles Treaty that forbade Germany 
from developing military aircraft. He significantly influenced aircraft design in Japan and the 
United States. Mitsubishi sent Japanese engineers to Germany in December 1921. Mitsubishi- 
Rohrbach GmbH, was formed in Berlin in June 1925. Aircraft were assembled in Copenhagen in a 
new factory built with the assistance of the Imperial Japanese Navy. The personnel from Japan 
learned aircraft structural design from Rohrbach. One result was the Japanese A6M1 fighter, the 
infamous “Zeke” or “Zero” that was the scourge of the Pacific during the first half of World War IL. 


Rohrbach was probably the first to use the term “stressed skin.” He had a major influence on 
American structural design because of a visit to the United States in 1926. He presented at least one 
lecture in Los Angeles and is credited with inspiring designers such as John K. Northrop to use 
stressed skin design in new airplanes, notably the Lockheed Vega. 


Figure 2.1.14-Northrop’s cellular wing design for his Northrop Alpha was the latest in wing design in the 
early 1930’s.(sales brochure) It was much lighter than other stressed skin designs of the time. 
(Introduction to Airplanes - Navy Training Courses, 1944). 


The development of advanced aircraft structures came to the United States from Europe in the 
1930’s. The NACA championed the use of metals and stressed-skin designs. Boeing, with its model 
247 pioneered advanced design and construction. The DC-3 with its clean smooth exterior 
aluminum stressed skin added to progress in the area. 


Once designers settled on thick wings, large internal volumes became available for storage. This led 
to wings with the ability to carry fuel (“wet wings”) and internal gun mounting, not to mention 
landing gear storage. By the mid-1930 the basic concept of stressed-skin construction was the fiat 
accompli for airplane wing design. 
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Figure 2.1.15-By the mid-1930’s the basic concept of stressed-skin construction had evolved to the point 
where it was the choice of designers of high speed, large aircraft. 


The quest for high speeds and light weight structures meant that the surfaces were flexible. 
Structural analysis attack two different problems, those related to stress analysis and those related to 
vibration and aeroelastic analysis. Stress analysis, being concerned with localized effects of 
loading, requires a very different model than aeroelastic analysis. Stress models are very detailed 
and pay attention to small internal details since they must model critical design details such as 
structural attachments. Aeroelastic models must only capture “macro” as opposed to “micro” 
structural features. These features include surface shape, both static and dynamic. Analytical 
odeling is the subject of the next section. 


2.2 Analytical modeling 


Analysis differs from synthesis. Synthesis involves putting things together. Analysis consists of 
breaking the system into parts to determine their interrelationships and interactions. Analysis 
requires identifying essential system details, separating the essential from the unessential, definition 
of models to reflect these details, and then connecting the details.’ 


Mathematical models, essential to all engineering projects, allow us to acquire information and then 
transform it into knowledge. Models predict outcomes and help us to decide among alternatives. 
Models reduce risk in making decisions and help us size critical dimensions as well as provide 
details for the airplane. 


Aeroelastic models must be predictive; they must provide details required to allow us to understand 
what happens to the structure when we fly at different speeds and altitudes. Given an input or set of 
inputs, models provide an output. Input might be airspeed, V, while output might be wing tip 
deflection or load distribution. In addition to input/output predictions, analytical models provide us 
with understanding of how different components of a system interact. Models provide similarity 
parameters that show how physical parameters couple together. 


A model may be physical, a wind tunnel model for instance, or the model may be analytical and 
simply consist of a set of equations. The model may resemble the actual engineering article very 
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closely or it may leave out 
unessential features. 
Figures 2.1.4 and 2.1.5 are 
examples of analytical 
models that represent the 
dynamic behavior of a 
flexible beam flexible aircraft, but one 


6 degree of 
freedom model 
(3 masses with 
longitudinal and 

rotational freedom 
linear translational 
and rotational 


springs) model with ee FT 
translation; very Closely resembles an 
extensional and airplane while the other 
bending degrees of 
ee does not. 


Models have objectives. 
What are we looking for? 
What do we want to know? 
Figure 2.2.1 — Analytical models come in a wide variety of forms, The answers to these 

applying symbolism to real configurations. questions determine the 
complexity of the model. 
Models also have assumptions that simplify their development but may limit or even destroy their 
validity in certain circumstances.* Everyone remembers D’ Alembert’s paradox in fluid mechanics in 
which the exclusion of viscosity and the boundary layer gives a beautiful, simple, but incorrect result 
for the drag force. 


Models are virtual, not real because some unnecessary details are left out of models. An example is 
the idealization or abstraction of a wing as a simple beam, even though the actual wing is a thin 
shell-like structure. This abstraction may be appropriate for estimation of wing natural frequencies 
but wholly inadequate for estimation of stresses in ribs. In our aeroelastic models we will begin with 
highly idealized abstractions that substitute springs for actual structural elements and simple wing 
planforms to simulate the aerodynamic features we are studying. 


The two most commonly used structural models are: 1) differential equation representations of static 
equilibrium, either in terms of internal forces or displacements; or 2) lumped parameter models with 
idealizations that lead to a finite number of degrees of freedom. Differential equation models, if 
solved numerically, for instance by finite differences, lead to a matrix representation of structural 
and aerodynamic effects. 


2.2.1 The typical section, the flying door and Rayleigh-Ritz idealizations 


With one exception, the models in this chapter are discrete element models such as that shown in 
Figure 2.2.1. We will also assume that the behavior of our systems is linear, that is, output is 
proportional to the input. 


All mathematical descriptions of aeroelastic phenomena use idealizations in which some things are 
emphasized and other factors are left out. The most common and popular idealization for aeroelastic 
phenomena are segmented models. These models, shown in two forms in Figure 2.2.2, resemble 
single or interconnected wing segments. 
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In its simplest form, the single “typical section” model resembles a wing with constant chord 
sandwiched between the walls of a wind tunnel. This typical section has pitching and “bending” 
displacement degrees of freedom; its angle of attack changes with the applied load and also deflects 
upward or downward. These displacements are resisted by two types of springs that are surrogates 
for the wing structure. The typical section model is extended to consist of a series of interconnected 
wing segments, elastically restrained as shown in Figure 2.2.2. 


When the typical section consists only of a 


single wing segment the chordwise dimension is Li ft | 


taken as the wing dimension at a position %4 of sin -__-h (plunge at 
the way out the semi-span. Similarly, other (oe = _ Shear center) 
critical dimensions such as the position of a an i 


point called the shear center and_ the 
aerodynamic center are taken from the typical 


section at the % chord. <= 
shear 
centers 

A second model used for static aeroelasticity <_ 

. . . . . A 

studies is a variation of the unswept typical ro 

section. The so-called “semi-rigid” model, 

shown in Figure 2.2.3, appears in research t 

reports from the Royal Aircraft Establishment at : 


Farnborough in England and includes sweep as a Pinu DD = Teewenh pica secon medals 


patameter, Ths model is used in Chapter 3 showing pitch and plunge degrees of freedom to 
when we discuss the effects of sweep of simulate wing torsion and bending. 


aeroelastic behavior. The swept wing model 
resembles a “flying door’” restrained at the wing root by elastic springs, but able to twist and rotate 
to simulate bending. 


The Rayleigh-Ritz model is a third type of simple analytical model. Like the other two models the 
Rayleigh-Ritz model has a limited number of degrees of freedom, but these degrees of freedom are 


FUSELAGE CENTERLINE 


| HINGE LINE 


TWISTING HINGE LINE. 
{ELASTIC AXIS} 


XY streamwise 
direction 


Figure 2.2.3 - Swept wing models with two degrees of freedom to simulate swept wing aeroelastic behavior. 
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deformation shapes with unknown amplitude. These unknown amplitudes become the degrees of 
freedom for the model. Unlike the other models Rayleigh-Ritz models have distributed stiffness. 
Rayleigh-Ritz models begin as continuous structures such as beams and plate structures, with 
theoretically infinite numbers of degrees of freedom that are then reduced, through the use of energy 
methods, to a finite number of degrees of freedom. This model is discussed in an ensuing section. 


2.2.2 — Functional diagrams and operators — modeling the aeroelastic feedback 
process 


One important outcome of our aeroelasticity studies is the recognition that aeroelasticity is a 
feedback process in which an initial aerodynamic form, a wing surface for example, is loaded 
initially by placing it at an angle of attack. As a result, it deforms to produce deformation that in 


turn produces new loads that 

then produce additional loads. 

This load/displacement é 
feedback process may or may Wing 
not converge. The process 


may be viewed as a simple 
block diagram, as indicated in 
Figure 2.2.4. 


Figure 2.2.4 — Wing functional diagram showing inputs and output 


In addition to the inputs, the wing lift depends on the wing deformation. As a result we view the 


model as _ having two 
interconnected elements: 1) a 
wing lift generating element 

that takes information about | pressures 

calculates deflected shape. 

This wing may be represented deflections 

functionally as a block diagram _ Figure 2.2.5 — Block diagram showing feedback connections between 
shown in Figure 2.2.5. aerodynamics model and structural model 


the wing shape; and, a flexible 
structure element that takes 
loads information and 


The calculation of aerodynamic pressures (and forces) requires that the wing shape and incidence are 
known inputs. The input angle of attack includes wing deformation. The calculation of the 
deformation is done in the flexible structure module and then the wing shape, incidence and 
aerodynamic pressures are updated. This process continues until the process converges. 


In practice, this iteration is seldom necessary. We can represent what goes on in each of the modules 
by mathematical operators. In our case, these operators are matrix relationships, but we are getting 
ahead of ourselves. The point is that we need to connect aerodynamic and structural elements to 
create an aeroelastic feedback system. Whether the system is operating subsonically, transonically 
or supersonically, there is a special operator or module for the aerodynamic pressures and loads. 
Similarly, whether we choose to represent the wing internal structure in simple or exquisite detail, 
there is a structures module. 
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Both modules connect loads and deformations. For instance, the aerodynamic loads distributed over 
a surface depend on the local incidence of the surface. This incidence is due to both fixed loads 


{P } due to known angles of attack and loads due to flexible deformations {x Ne The loads P; are 
due to initial aerodynamic loading such as an angle of attack of the aircraft input by the pilot. The 


expression for these forces and moments distributed over the wing, eT is written as 


{Le"}=[ A, |{x }+{P} where the matrix [ A; | is a collection of aerodynamic influence 
coefficients. The relationship between the external structural forces and moments and the same 
deformations {x i is written as {ie = LK, {x i} where the matrix | K a is the so-called 


stiffness matrix, a matrix of structural influence coefficients. 


Connectivity between the structures and aerodynamic modules is accomplished when we write 
equations of motion or, in the case of statics, equations of equilibrium, an equation that sums internal 
and external forces and moments and sets them to zero. This provides an equation with the form 


tL }-1L}=0= LX, l{x}-L4Ay Hes} -{23 =0 
Solving for the deformations {x,} we have Sia = [ K i |tx i} . The load distribution found to 
be {Li} = [4 |[[y ]- cule {P}+{P}. 


The aeroelastic problem is called a statically indeterminate problem because structural deformation 
creates loads that enter into the equilibrium equations. The feedback process is contained in matrix 
operations such as addition and inversion so we do not need to iterate. The solution procedure is 
similar if we have differential equations to represent the structure and the aerodynamics although the 
operation for inversion or “inverse operators” is somewhat different. 


The next few sections discuss and develop simple structures and aerodynamic modules for 
aeroelasticity models. Refer back to the present section to see how the general process just outlined 
is applicable to these simple models. 


2.3 Matrix structural analysis — stiffness matrices and strain energy 

The demands for performance, efficiency and high speed for modern aircraft would be difficult were 
it not for high-speed computers and modern automated analysis procedures such as finite element 
(FEM) structural methods and aerodynamic paneling methods and computational fluid dynamics 
(CFD) efforts. These advanced analysis methods provide high fidelity analytical models to aid 
design decision-making and shorten the time required to make decisions. This section presents an 
overview of a structural analysis approach known as matrix methods and shows an example of how 
we develop them; it also introduces the shear center and elastic axis concepts. 


As discussed in the previous section, the external forces and moments applied to an aircraft {F;} 
and displacements {u;} of the aircraft structure are related to each other. If the structural behavior is 
static and assumed to be linear elastic, this relationship is written as 
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{Fi} =1K,]{u;} (2.3.1) 


where [K i] is the n x n system stiffness matrix. The elements, K,,, are called stiffness influence 


i ’ 


coefficients. 


To develop the relationships required for the calculation of the stiffness matrix we use the finite 
element method. Finite element models, such as that shown in Figure 2.1.5, are collections of 
mathematical elements joined together at “node points” where displacements are computed and 
forces and moments are applied. A node point is a junction point connecting two or more elements. 
Finite element node points may or may not correspond to physical joints on the actual structure. For 
simple trusses, the physical node points and the mathematical node points are identical. The 
relationship between the structural forces and displacements on each of these elements is expressed 
as a Set of linear (or nonlinear) algebraic equations whose coefficients are referred to as the stiffness 
matrix, a matrix of structural influence coefficients relating external forces and moments to the 
displacements at node points.. 


In reality, aircraft structural systems contain groups of discrete elements connected by rivets, welds 
and joints. The complex geometry and the wide variety of loads applied makes the analysis task 
challenging. Differential equation models with smooth continuous features are not effective for 
complex structures since a fundamental assumption for such equations is that the system is 
continuous, at least in discrete areas. 


Finite element methods belong to a class of structural models called matrix models. Matrix method 
analysis development began in the aircraft industry with flutter and vibration clearance efforts. 
About 60 years ago, the need to accurately represent aeroelastic and vibratory interactions of 
geometrically complex flexible airplane structural elements at the earliest stages of design (and the 
ability to use electronic computers) led to the development of the finite element method (or Direct 
Stiffness Method) of structural analysis. 


The finite element method began at the Boeing Company in the 1950’s as the result of the quest for 
accurate analytical procedures for vibration and stress analysis on complex airplane airframes.” 
Within a few years, this idea had exploded onto the analytical scene and became the mainstay of 
modern structural analysis. 


The finite element method has evolved rapidly over the past thirty years and led to accurate elements 
and large scale computer codes such as NASTRAN and ANSYS. These codes are widely used by 
both aerospace and non-aerospace companies. 


The structural idealization of an aircraft structure and the definition of the loads on the structure 
require that we define degrees of freedom of structural motion. These degrees of freedom may be 
actual displacements (including rotation) of the structure or they may be amplitudes of assumed 
continuous displacement functions (shapes) or vectors that describe the structural deformation. The 
final result is that the actual structure, with an infinite number of displacement degrees of freedom, is 
replaced by an idealized model with a manageable number of unknown displacement quantities or 
degrees of freedom such as shown in Figure 2.1.5. 
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Aircraft loads are created by combinations of distributed pressures and inertia and relatively 
concentrated forces or inertia (such as engines). For matrix analysis, both the distributed and 
concentrated loads are converted into a set of equivalent discrete forces and moments associated 
with each of the degrees of freedom. The result is a set of generalized forces (the term generalized 


force can refer to either a force or a moment) denoted as {F;}. Each element F;, is associated with a 
corresponding deflection vector element u; (these generalized displacements can also be either a 


deflection or a rotation). 


In a system without dissipative outlets such 

as friction, total system energy is conserved. origin 
Therefore, all of the work done by the — J ; yr 
external forces must be stored within the 
structure. If the structure is restrained 
against rigid body motion and the forces are 
applied slowly, the work done is stored as 
strain energy or potential energy. It can be 
recovered when the external loads are 
removed. 


Figure 2.3.1 - Spring-mass-force system 
Consider the simple spring element shown in 


Figure 2.3.1. If the external force is applied very slowly so that the acceleration is so small that 
static equilibrium is maintained approximately throughout the loading process, then the work done, 
W, by the force F' in moving through displacement u is equal to the internal strain energy stored. 


en nen re =U (2.3.2) 
2 2 


In the more general case where there are many displacements and many external forces F; located at 
the node points where displacements u; are measured, the work done by the external forces is 


ae ariel? 
Wehbe) = a te) (2.3.3) 


The factor % appears in Eqns. 2.3.2 and 2.3.3 because the relationship between the loads and 


deflections is linear and the loads build slowly from zero to their final values, F; and u;. 


Substitution of Eqn. 2.3.1 into 2.3.3 provides an expression for work and strain energy done in terms 
of the structural displacements and the system stiffness matrix in Eqn. 2.3.1. 


W =(u,IK, Ha) (2.3.4) 


The strain energy, U, is: 
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1 
U =F lu, VLK, lu} (2.3.5) 
When differentiated with respect to u; the strain energy U provides force equilibrium equations 
“= SK i =F. (2.3.6) 


Differentiating Eqn. 2.3.6 with respect to u |; We recover the stiffness matrix elements for the rows 


and columns representing the equations of motion. 


2 2 
PU _, _ FU 


——_=K, =——_=K. 231 
Juou;  " Ouou, ” 


since the order of differentiation in Eqn. 2.3.7 may be interchanged. The structural stiffness matrix 
is symmetric. The stiffness matrix elements may be positive, negative or zero when i # j , but the 


main diagonal elements must be positive. The determinant of the stiffness matrix is also positive. 
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2.4 An example - Construction of a structural stiffness matrix — the shear center 
concept 

The idealized two-dimensional "structure" shown in Figure 2.4.1 has two linear elastic springs with 
stiffness K, and K, (with units lb/in or N/m). The deformation state of this model is defined by the 
downward deflection, h, at an arbitrary reference point chosen by the analyst, and a pitch rotation, 6, 
about the reference point, measured positive counterclockwise. The external force, P, computed at 
the reference point, is the resultant of all forces and pressures acting on the structure, while M is the 
resultant of the moments, about the reference point, of all forces and pressures acting on the 
structure. We will compute the stiffness matrix that relates P and M to hand @. 


Reference datum level 


Reference datum level 


downward force = P 


b 


Figure 2.4.1 - Two dimensional bar and F,=K(h+b0) F=K(h-a) 
spring system Figure 2.4.2 - Two dimensional bar free body 
diagram 

We will use force and moment equilibrium equations to derive matrix equations of static 
equilibrium. Consider the free-body diagram of the deformed system (shown in Figure 2.4.2), 
including the forces internal to the springs. When @ is small (sin @ = @), the compression of the 
right hand spring K; is h —a@, while the compression of the left spring Kz is h+b@. Note that the 
dimensions a and b measure the spring locations relative to the reference point. 


The spring restoring forces are proportional to K,,K> , h—a@ and h+b@. Summing forces and 
moments about the reference point we have two equations for static equilibrium. 


> F =P-K,(h-a0)-K,(h+b0)=0 
>| M =M -K,b(h+b0)+ K,a(h—aé) =0 


These two equations are written in matrix form as: 


| | (K,+K,)  (K,b-K,a) (| 
= (2.4.1) 
M) |(K,b-K,a) (K,b’ + K,a’) || 


Equation 2.4.1 relates the applied external loads (a force P and moment M) at the reference point 
to the deflection and rotation at this point. Elastic or static coupling is indicated by the non-zero off- 
diagonal terms (K,b— K,a); in Eqn. 2.4.1. The force P, produces both a @ rotation and a 
displacement, h. The moment, M also produces both types of displacement. 
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We can also use an energy approach to derive this matrix also. 
o>, are functions of h and @ 

6, =h-asind =h-aé 

6, =h+bsind =h+b@ 


The strain energies stored in each spring are 


1 1 

U; =5 KF na (h—ad) 
1 1 

U, => 6; = 5K, (h+b0) 


The deflections at each spring, 6, and 


The total strain energy is U =U,+U,. Differentiate the strain energy expressions with respect to h 


and @to get: 
0U 0U 
“= Ki lh 29) 59 7 aK (A— 48) 
aU, aU 
= K,(h+b0 “ 
Oo = K, (+60) 


Differentiating the strain energy for a second time we have: 


aU, 07U, | 
oh? d@h 

OU _ é OU _ 
oh - onde ? 


=K, —aK, 


5 = bK,(h+bé) 


2 
= 
o°-U 
ag KP 


Using these expressions with Eqn. 2.4.1 we see that the stiffness elements generated using the 
energy approach are identical to those generated using free body diagrams and static equilibrium 


equations in Eqn. 2.4.1. 


Elastic coupling (also called static coupling) depends on the coordinate system chosen. If we 


change the location of the reference point, the dimensions a 


and b will change. The size of the 


resultant force P (it is the sum of all the vertical forces acting on the structure) remains the same, 


but the size and the sign of M will change, since it represents the net external moment of a 


combination of external forces about the reference point. From Equation 2.4.1, if K,a is equal to 


Kyb, the elastic coupling term (K,b—K,q) is zero. Elastic coupling disappears if the ratio of 


dimensions a and b is: 
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Equation 2.4.2 shows that there is a unique point on the idealized structure where a concentrated 
force, P, causes only a displacement, h, while a torsional moment, M, causes only rotation about that 
point. The first condition defines a point called the shear center. The second condition defines a 
point called the center of twist. For a linear elastic structure the position of the shear center and 
center of twist are identical. 


The shear center is defined as “a point on a structural configuration where loading by a concentrated 
force will create only translational displacement but no rotation or twist.” If K, = K,=K, Eqn. 
2.4.2 shows that the center point of the bar (a= b) fits our shear center definition. If the two 


springs have different stiffness values, the shear center is located away from the center of the 
structure. 


We can also solve for the shear center position by first using Eqn. 2.4.1 to solve for h and @ when 


P and M are applied. 
hl. (K,b° + K,a°) —(K,b—K,a) |{ P (2.43) 
0) <A| -(K,b-K,a) (K,+K,) |(|M@ 


where A= K,K, (a? +b° + 2K, K, (ab) is the determinant of the structural stiffness matrix. Note 
that the determinant is always positive. 


To find the shear center we apply only the load P (we set M=0) and then find that 


2 2 
{i _P +K,a ’ eae 
0) A|-(K,b-K,a) 


From Eqn. 2.4.4, when K,b = K,a then @=0, no matter what the size of the load P. This is the 
same result found in Eqn. 2.4.2. The shear center does not depend on the applied load magnitude. 


2.5 Subsonic aerodynamics - fundamentals 


This section reviews the fundamentals of low speed aerodynamics to the extent necessary to develop 
basic aeroelastic models of the types shown in Figures 2.2.2 and 2.2.3. We will review terminology 
and definitions only to the extent necessary to develop these models in this chapter. The review 
follows closely the material contained in Chapter 3 of A.C. Kermode’s Mechanics of Flight."° In 
this section we will discuss the mathematical representation of lift and reference points such as the 
center of pressure and the aerodynamic center. 


Analytical aerodynamic development occurred primarily in Britain and Germany before, during and 
shortly after World War I. Strangely the first two decades produced a great deal of discussion about 
the origins of lifting forces.'' Without this understanding it was difficult to develop analytical 
theories. As a result the engineering world was divided into two groups, often referred to, in Britain 
at least, as the “mathematicians” and the practical men.”* In this section we will take the view of the 


* David Bloor, (The Enigma of the Airfoil: Rival Theories in Aerodynamics 1909-1930, University of Chicago 
Press, 2011) presents the history of the struggles between researchers in Britain and Germany as they sought to 
unravel the mysteries of aerodynamic lift. 
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“practical man” (with apologies to the young women who contribute to the field — much has changed 
since those early days when half the population was sociologically barred from entering engineering) 
and approach the review of aerodynamics from an empirical viewpoint. 


Figure 2.2.2 shows a mathematical model of a wing composed of two different sections or segments. 
In our development, the term “wing segment” is used to mean a “shaped surface, such as a three- 
dimensional airplane wing, tail, or propeller blade that produces lift and drag when moved through 
the air.”’’ Figure 2.2.2 also has a cross-section of a wing. This cross-section is referred to as an 
airfoil. The term airfoil (or aerofoil) is defined as “the cross section of a body that is placed in an 
airstream ... to produce a useful aerodynamic force in the most efficient manner possible.”'® An 
airfoil may be thought of as a wing with infinite span. In this chapter and others to follow, the term 
wing, lifting surface or wing segment refers to three dimensional surfaces while the term airfoil is 
used to denote the two-dimensional wing cross-section. 


5 a A mean 
Parameters used to define wing cross- lias camber line 


sectional (airfoil) features and ; - 
A : ‘ leading trailing 
terminology are shown in Figure 2.5.1. edge erg 


The straight chord line connects the 
rounded leading edge with the sharp 


trailing edge of an airfoil designed for chord c 
subsonic flow. The length of this line is Figure 2.5.1 — Airfoil cross- section with geometrical 
denoted as the chord dimension c. The definitions 


mid-point between the upper and lower 
wing segment boundaries is called the mean camber line. Wing camber is measured by calculating 
the ratio between the maximum distance between the mean camber line and the chord line, measured 
perpendicular to the chord line, and then 

dividing by the chord and multiplying 

by 100 to give camber as a percentage 

such as “5% camber”. For symmetrical 

wing segments, the chord line and the a 
mean camber line are identical so / 
camber is zero. 


The airfoil angle of attack, a is shown : 
: : : : Veo \ Relative Wind SS ti 

in Figure 2.5.2, together with the lift and _ = ss 

drag. Note that by convention the lift Figure 2.5.2 — Wing cross-section showing lift and drag 
always acts perpendicular to the flow forces and angle of attack as measured relative to the flow. 
direction. Since a wing may be 

composed of many different wing sections, the angles of attack of individual wing cross-sections 
may not be identical. 


The angle of attack is sometimes confused with the aircraft pitch angle. As indicated in Figure 
2.5.3, the pitch angle is measured from the same reference line as angle of attack (the chord line), but 
with respect to the horizon. Only when the aircraft is in straight and level flight are the aircraft pitch 
angle and wing angle of attack identical. 
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In subsonic flow, deflecting 
moving air around a_ surface 
traveling with respect to the air 
creates aerodynamic lift and a 
pitching moment, as shown in 
Figure 2.5.4. The aerodynamic Figure 2.5.3 — Pitch angle (shown as 8) and wing angle of attack , 
lift, drag and pitching moments @ are not the same. http:/www.aerospaceweb.org/ 

are functions of the aerodynamic 

pressures acting on the lifting surface. These pressures are, in turn, functions of the wing segment 
angle of attack. 


—_———— 
— — d — 
— Increased ~ 
— —— — ~~ — 
: F —_— ——_ speed — ee, 
Lift is created b a eo : - Se 
y os ret airflow~—=— SS 


phenomenon known as 
vorticity, a topic that we will 


_ — Sa ee — 

Pie Ra eee —— 
n ~ ~ 
oe : ~~. Le * 
VOL —-—" Decreased. — ate 


Tet 
WA 


a 
discuss in Chapter 4. The key _—— __speed Do 
features of Figure 2.5.4 are: 1) en of airtiow WE 

. —_——, Se meee 
the stagnation point on the a ama cia a 


leading edge where flow slows 
to zero; 2) the upper surface Figure 2.5.4 - Airflow around a wing section in subsonic flow. 
region in which flow speeds up; (Kermode, p. 72). 

3) the lower surface region 

where flow slows; and 4) the trailing edge where the two flows come back together and have the 
same velocity. 


: 1 : 
Bernoulli's law requires that the sum of the static pressure and dynamic pressure SG pv , where / is 


the air density and V is the airspeed) be constant in an incompressible fluid flow. As the air speeds 
up over the upper surface the fluid static pressure must decrease so that a low pressure region forms 
there. The fluid speed is highest and the static pressure lowest near the upper surface of the forward 
part of the wing cross-section. On the bottom surface, the air slows, the static pressure increases and 
the net force upward increases to create lift. 


Wing lift and drag and pitching moment depend upon: wing cross-section shape; wing planform 
shape, and, dynamic pressure, /2 pV’. The lift, drag, and pitching moment are denoted as L, D, and 
M, respectively; these three quantities are defined in terms of wing surface planform area, S, 
dynamic pressure, g , and wing mean aerodynamic chord C and aerodynamic coefficients as: 


L=C,qS (2.5.1a) 
D=C,qS (2.5.1b) 
M =C,,qSt (2.5.1¢) 


The force and moment per unit length or per unit span and are defined as: 
l=aqce (2.5.2a) 
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d=c,qce (2.5.2b) 
m=c,qce- (25.26) 


These aerodynamic coefficients are functions of the wing angle of attack, Reynolds number, and 
Mach number. 


The lift and drag coefficients for an example airfoil are shown in Figure 2.5.5. At angles up to 12° 
this plot of the lift coefficient is almost a 

straight line. At about 15° the lift coefficient pe 

reaches a maximum and then declines with 
further increases in angle of attack as the HHH 
wing segment loses lift or stalls. 


Stall is associated with the breakdown or 
detachment of the streamlines around the 
wing. The shape of the cross-section makes 
little difference in the angle at which stall 


42m@H= TION 45er 
+am-HeAsROn A> ZO 


occurs, although it does affect the value of 
the maximum value of C;, (or c;) achieved 


when stall occurs. Stalling is a function of 
angle of attack, not airspeed. 

Figure 2.5.5 — Lift and drag coefficients c, and cp as 
The linear portions of the airfoil lift and functions of angle of attack for an example airfoil. 
pitching moment curves are approximated 
by the expressions 


c,=c, +| — |@=c, +c,a@ 2.5.3 
— = C. 
i is Ae i L (2.5.3) 
dc, 
Cn = Cm + A=c,, +C,, a (2.5.4) 
° 04 


The notation ( a indicates partial differentiation with respect to angle of attack, @. The equation 


for C,, versus 0 depends upon the wing reference point location chosen to measuring the pitching 


moment. 


Wing drag is due to friction drag on the exposed or wetted area of the surface (this drag is also 
called parasite drag). Drag on a wing (as opposed to an airfoil section) has an additional component 
called induced drag that depends on the lift coefficient magnitude. The analytical behavior of the 
drag coefficient is approximately quadratic for angles of attack within the linear region of the lift 


curve. This relationship is: 
AC 
C=O eat a (2.5.5) 
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where the angle a is measured from the position of zero lift. The term C,, is called the parasite or 


zero lift drag coefficient. Eqn. 2.5.5 is also written as: 
C=C, OC; (2.5.6) 


The second term in Eqns. 2.5.5 and 2.5.6 is the induced drag or drag due to lift. This induced drag 
only appears for wings with a finite aspect ratio (aspect ratio is defined as the span squared divided 
by the wing planform area) and becomes very small as the wing planform becomes more slender. If 
the aspect ratio is infinite, the induced drag is zero. This corresponds to the results for a two 
dimensional airfoil. 


2.5.1 Reference points — the center of pressure 


The wing segment or airfoil pitching moment is an integrated product of pressure and distance from 
an arbitrary reference point. Pitching moment depends on angle of attack and the point at which we 
measure the pitching moment. One position for measuring pitching moment is the center of 
pressure. At this point the resultant moment is zero. The pitching moment coefficient at any point 


x (x is measured aft of the wing leading edge) is C,, = C,x+ Cy, 


leading 
edge 


The condition Cy, (x) = 0 defines the center of pressure position. Measured as a fraction of the wing 
segment chord aft of the leading edge, the center of pressure location is denoted as Xcp and given by 


the following equation, 
—-Cy 


x leading 
x, =P = 2.5.7) 
Cc Cc 


L 


Figure 2.5.6 plots the center of 
pressure position on a cambered 
airfoil as angle of attack changes. 
This plot indicates that at negative 
(nose down) angles of attack the 
center of pressure is near the rear 
of the airfoil. As the airfoil angle 
of attack increases, the center of 
pressure moves forward to near 
the 0.30 chord position and 
remains there until near stall 
when the center of pressure 


9 angle 


£ 
a 
~ 
“ 
I 


i fe) 
moves aft again as angle of attack lace ne = r = ae ag 
increases. Figure 2.5.6 also edge Angie of attack 
shows data for a _ different Figure 2.5.6 — Airfoil center of pressure & pitching moment 
reference point, the aerodynamic coefficient at the aerodynamic center (Kermode) 
center. 
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2.5.2 A different reference point - the aerodynamic center 

There is a special point on an airfoil where the aerodynamic pitching moment does not change with 
angle of attack. This point is called the aerodynamic center. The aerodynamic center is the point on 
the wing chord where the pitching moment does not change with angle of attack. We write this 
relationship as: 


ICy 
acy OC, x oe 
— t= C =0= L AC 4 pc 2.5.8 
OO Me Oa Cc OOL ee”) 


Solving Eqn, 2.5.8 for x, = wf. the aerodynamic center location aft (downstream) from the 


leading edge, as a fraction of the airfoil chord, 


ECy 
ao 
om C,, (leading edge 
Rye = = - Ma a add (2.5.9) 
L 
0a "s 
Hit , ke ee 
Both C, =Cy, and ——~ are negative, so Xj. iS positive. 
on Oa 
edge 


To understand the difference between the center of pressure and the aerodynamic center, consider 
Figure 2.5.7 which plots the pitching moment coefficient on a cambered airfoil as its angle of attack 
increases. The pitching moment coefficient measured at the leading edge (the lower line) decreases 
as angle of attack increases. Meanwhile, measured at a point near the trailing edge, the pitching 
moment coefficient increases as angle of attack increases. 


Figure 2.5.7 shows that the two straight 70-2 
lines measuring the pitching moment 

about both the leading edge and a point #04 
near the trailing edge intersect. This 
intersection point is special because it 


oO 


occurs at zero airfoil lift. This special 
value of the pitching moment coefficient (it 
is not zero if the airfoil is cambered) is the 
pitching moment coefficient measured 
about the aerodynamic center position 
defined in Eqn. 2.5.9. If the airfoil is “Oe 
symmetrical, the pitching moment at zero 


when the airfoil angle of attack is zero. “O-3k, 
- 16* 


Angie of attack, a 


< 
ot 
"Moment coefficient about 


f Angle of zero lift 


oieg 


Moment coefficient, Cy 


Returning to Figure 2.5.6 we see that the Figure 2.5.7 — Airfoil pitching moment coefficient about 


aerodynamic center position for the three reference points (Kermode). 
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example airfoil is near the airfoil 4 chord position and that, at this position, the pitching moment is 
slightly negative, but constant until near stall. Jf the airfoil is very thin and the flow is 
incompressible, aerodynamic theory places the aerodynamic center exactly at the % chord position. 
In addition, for symmetrical airfoils, the constant pitching moment is identically zero. 


The major difference between Eqn. 2.5.7, which uses the total lift and pitching moment about the 
leading edge, and Eqn. 2.5.9, which uses aerodynamic coefficient derivatives (rates of change 
with @) with respect to the angle of attack, a, is that the reference position changes when we use the 
center of pressure, even though we eliminate the aerodynamic moment from static equilibrium 
calculations. Aeroelastic problems are concerned with the effects of changes in aerodynamic loads 
due to lifting surface deflections. As a result, the early aeroelasticity researchers found using the 
aerodynamic center as a reference point for aerodynamic calculation more useful than the center of 
pressure. We will also use this point as our reference in the models developed in the next section. 


2.6 Lifting surface flexibility and lift generation-the typical section model 
In this section we introduce the typical 
section model used by aeroelasticians 
for at least eight decades. The model angle 


lift=constant 
structural force 


= structural moament 


used in this book differs slightly from 5 

this historic model in that it has a finite 

span, as will be explained in this 0 center of twist 
: structural 4 t 

section. twist Shear Genter 


low speed flight 


To begin, consider Figures 2.6.1 and 
2.6.3. which show two cross-sections. 
Figure 2.6.1 shows a symmetrical wing 
segment with the lift acting at the 
aerodynamic center, but with no pitching moment about that point. Figure 2.6.2 shows the same 
segement but with camber so that the pitching moment about the aerodynamic center is not zero. 
Both segments have a length dimension (span) that goes into the paper; these segments have a 
planform area, denoted as S. The wings 
section is mounted in the wind tunnel reduced 
and securely restrained by a structure pitch - 
that resists upward movement and angle Mo aie 

ae _ structural moment 
twisting. | ; = 


Figure 2.6.1 - Wing section twist in low speed flight 


litt=constant 


: ; ; ; structural center of twist 
Consider the symmetrical section first. twist 6 shear center 
As indicated in the free body diagram (Constant) ; 
shown in Figure 2.6.2, the pitching high speed flight 
moment about the aerodynamic center is _ Figure 2.6.2- Wing section twist in high speed flight makes 
zero since the section is symmetrical, up most of the required angle of attack. 


but the wing will twist because the 

resultant lift vector at the aerodynamic center does not pass through the shear center. The total wing 
segment incidence is the sum of the elastic twist angle @ and the angle of attack o, input by the wind 
tunnel operator. The sum of these two angles is small, only a few degrees. 
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If the wing lift is equal to the airplane weight then the lift is constant when the airplane weight is 
constant. As a result, the wing lift produces a twisting moment about the wing cross-section shear 
center that is also constant, no matter what the airspeed. If we go faster, the dynamic pressure 
increases and lift will increase if we do not reduce the wing angle of attack, O,. 


Let’s now add wing camber and consider the idealized unsymmetrical wing model with planform 
area, S, in Figure 2.6.3. The structural 
stiffness is represented by a spring, Ky), Lift. L 


that resists wing upward (bending) 


translation and a torsion spring, K,, that Co” ~ Kr. ae (plunge at 


nis =~ shear center) 
resists wing twist. Historically, the = << 


bending deflection has been called plunge 
(denoted as / in Figure 2.6.3 and measured 
at the shear center) with the positive 
direction is downward at the shear center. Figure 2.6.3 — Unsymmetrical wing with pitch and 


The shear center thus becomes our plunge degrees of freedom and pitch and plunge 
springs, showing offset distance, e, between lift and 
shear center and pitching moment at the aerodynamic 
center 


structural reference point. 
This structural model is similar to the 


model used previously in Figure 2.4.1. Here, we model the action of the two springs in Figure 2.4.1 
as a single spring to resist torsion and a single spring to resist plunge. As a result, 


K, = K,+K, (plunge spring) @6:1) 


_ K,+K,(a+by 
K,+K, 


K, (pitch spring) (2.6.2) 


When the operator gives the wing an angle of attack, @, this model generates aerodynamic lift ZL and 
a pitching moment M (about the shear center). L and M are functions of both the known quantity @ 
and the unknown twist, @ 


L=qSC,, (@) +9 (2.6.3) 
Myco= Myc + Le (2.6.4) 

or 
Myce = qeSCyac + GSeC,, (Gy + 9) (2.6.5) 


The two external loads in Eqns. 2.6.3 and 2.6.5 are written in matrix form as: 
> SC, oe + gSC. - + gScC, : 
= a c 
Mc Qo, 0 ella qo, a, , g2CCuac 1 (2.6.6) 
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The static equilibrium equation that defines the relationships between the aerodynamic force and 
moment (L and M) and the shear center deflection, h, and the twist rotation, @, is: 


0 x lla" as 


The loads expression in Eqn. 2.6.6 is substituted into the static equilibrium equation, Eqn.2.6.7, to 


give: 
fo Oh 0 -I|{h =i 0 
ar cee = qSC,, O26 + qSC,, a, : + gSeC ye 1 (2.6.8) 
T 


Gathering all of the deflection dependent terms on the left hand side of Eqn. 2.6.8, we have 


ca Ura Coie | gd a ome ves etree come te (2.6.9) 
Oe Re lele = a0) eel = eee he QOC/MAC } Ys 


In Eqn. 2.6.9 the right hand side contains the inputs and the left hand side has the unknown output, 
@. The 2x2 matrix proportional to aerodynamic pressure, g, and the section displacements is called 
the aerodynamic stiffness matrix. While the structural stiffness matrix is symmetrical and h and @ 
are decoupled, the aerodynamic stiffness matrix is unsymmetrical and h and @ are coupled. 


To solve for h and @, first combine terms and then divide each term in Eqn. 2.6.9 by K,. 


Ve one h) SC, @, [-1 0 
r r |S by t | + a | (2.6.10) 


0 1 eG, 6 K, e K,; 1 
K; 


The wing plunge displacement, h, and twist angle, 0, are found by inverting the matrix on the left 
side of Eqn. 2.6.10 and multiplying both sides by this inverse. 


-9SC;, -gSC;, / | 
1 K, Hy, K, 
K Sec K SeC 
h| _ q8C,,& tee fe “Te gScCyac ee he ° 
= T + T 
0 K; 


1 e K,; 1 


1 PeG, 1D, 
K; K, 


or 
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qSeC, —qSC, 
ay, eee oe ee ce 
K, gSeC gSeC 
h| _ qSC1,% 1— be , ScCyac | 1- ie (2.6.1 1a) 
6 K, 


= 7 eee 
i gSeC,, es gSeC,, 
K, K, 


This relationship is also written as: 


—qSC, 
Se ey 
1 sala A af“ 
_— ‘a qd e 
(et = K; id pee Cyie | l= ye 


2.6.11b 
gseC,, [Ke ea 


—— 1-7 
1" ‘ye K; 
£ 


The plunge displacement, h, and the wing twist angle, @, are composed of two parts, one due to 
section angle of attack and the other due to wing camber: 


ey Oy, 
K gScC K 
h= +_ la, wae o 2.6.12 
1, Pe, : K, 1, BG, ( ) 
K, K, 
qSeC,, 
K ScC 1 
O= r qee™ mac (2.6.13) 


———_1—__|a,,+ 
1 Be, . K, 1, BG, 
K, K, 


Substituting the expression for elastic twist, Eqn. 2.6.13, into the expression for lift, Eqn. 2.6.3, we 
have the expression for wing lift, again composed of two parts from two sources. 


qSC, gSC,, [Scene 


= air 
1 e"n, ° 1 79°", K 
K, K, 


Note the presence of a nondimensional “aeroelastic parameter” in Eqn. 2.6.13 and Eqn. 2.6.14. 
Both Eqn. 2.6.13 and Eqn. 2.6.14 have denominators that tend to zero as the nondimensional 


(2.6.14) 


T 
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SeC 
aeroelastic parameter cone K. tends to 1. The value of dynamic pressure at which this occurs, 
T 


Ip =", Cc? is called the divergence dynamic pressure and is a subject of discussion and 
ec, 


derivation in ensuing sections. With low torsional stiffness K,, the term in the denominator of Eqn. 
2.6.14 can be small, creating large deflections, even for small angles of attack. If this is so then the 
wing response (plunge and twist) will be very large for small loads. 


When the wing cross-section is symmetrical (C),,- =0) the lift expression simplifies to a single 


term. 


= gSC,, 


=| ———_—___*—__ |q@ 
1 DeCi, . 
K, 


The required angle of attack to generate lift equal to weight is 


a, = oul = aie (2.6.16) 
gSC,, Ke 


In this case the wing twist is constant, as discussed earlier: 


(2.6.15) 


gSeC, gSeC,, 


K. K. L qSeC 
re nc Se) es ee Pf lar 
1 9SeC,, 1 95C, qgSC,, a aa 
K, K, 


K 
Notice that, from Eqn. 2.6.16) the input angle of attack is zero when q=q, mers In this 
e L, 
special case all the wing lift is generated by structural twist, @ The reasons for this strange result 
will revisited when we discuss wing divergence in a later section. 


A comment on the typical section model used in this text 


At the beginning of this section we mentioned that the typical section used in this text is slightly 
different than the typical section used by early aeroelasticians. If the wing is nonuniform, with a 
tapered chord for instance, then aeroelasticians would go to a point located at the wing %4 span 
position. They would then measure or estimate aerodynamic data for that section, in particular, the 
section chord, the section lift curve slope, section aerodynamic center position. The section would 
have a unit span so its planform area is simply S=c*1. 


The difference between the historical model and the model we have discussed is that the structural 
coefficients were stiffnesses per unit span. In our case then these stiffnesses would be 
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Ky and aA where the dimension b is the semi-span. Placing these terms into our formula 


then we have gp) = Ve Cc, = Tee Cs The result is the same if the wing lift curve slope and 


the %4 span lift curve slope is the same and the wing has no taper. From an educational standpoint 
the author prefers to see a wing slab moving in the air and supported by concentrated as opposed to 
distributed springs. This subject will come up again in Chapter 4 when we consider the dynamics of 


the wing. 


A sidebar — why the %semi-span position? 


Why was the geometry and aerodynamic properties of the 34 span position thought of as “typical?” 


Tupolev 154 wing example 


Fuselage 
side 


Rrneness 


Fuselage : 
centerline Inboard flap 


Figure 2.6.4 - Tu-154 wing structural 
idealization with two flaps 


0.4 0.3 0.5 0. 


Figure 2.6.5 — Tu-154 wing twist with respect to 
the freestream for three loading conditions, 
nondimensionalized with respect to the tip twist in 
each case . 


To determine what is special about the %4 chord 
position a group of Russian researchers examined 
the static aeroelastic behavior of a_ tapered 
Tupolev 154 swept wing with two flaps shown in 
Figure 2.6.4. Three separate load conditions 
were used to compute three different wing twist 
distributions. The first condition was a one-g 
loading, the second was deflection of the inboard 
flap while the third condition was deflection of 
the outboard flap. The twist distributions for the 
load three cases are shown in Figure 2.6.5, 
nondimensionalized so that the maximum value 
in each case is 1. There is nothing remarkable 
about the twist values at the %4 span position. They 
are all different, as we would expect. 


When the wing twist in each case is divided by the 
integral of the twist distribution (different for each 
case) a distribution defined as 


515). 42) 
6, (z ) maT: 

[, 4 (2a 
2.6.6. The integrals are a measure of the total twist 
in each of the three cases and are different in all 


three cases. In this case the curves intersect near 
the % span point. 


(i=1,2,3) is plotted in Figure 


Figure 2.6.6 shows similar plots and results for a low wing span supersonic transport for five 
different loading conditions. In this case the curves intersect near the 0.82 semi-span position. 
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Figure 2.6.6 — Tu-154 streamwise twist 
for three load conditions plotted as 
twist divided by the integral of twist. 


Figure 2.6.7 — Russian SST 
streamwise twist for five different 
subsonic loading conditions, 
nondimensionalized with respect to 
integral of twist. 


2.7 Example problem — The effect of elastic axis offset on wing lift effectiveness and 
divergence 


The idealized wing whose cross-section is shown in Figure 2.7.1 will be tested in the wind tunnel at 
several different airspeeds at sea level 
conditions and at such low speeds that the 
flow field is incompressible. The position of 
the wing box indicated in Figure 2.7.1 can 
be moved to change the wing shear center 
position and the offset distance between the 
shear center and the aerodynamic center. 
For instance, when the shear center is at the 


mid-chord, when the shear center is at the shear center range 

wing mid-chord the offset distance e is equal Figure 2.7.1 - Wind tunnel wing segment model 
to c/4. Wind tunnel tests show that the wing cross-section 

divergence dynamic pressure is 


100 lb/ cg at this wing box position. 


The aerodynamicist predicts that the ratio of C,,,. to C — is -0.075 and the wing lift curve slope is 


4.0; the reference planform area is S =10 fae The section angle of attack @, for all tests will be 


pre-set to 3°. 
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Problem 
a) Develop the expression for the lift force when the shear center is located at the mid-chord. Find 
the twist angle @ as a function of wind tunnel airspeed. Plot lift and twist angle vs. airspeed. 


b) Consider the same wing segment with different wing shear center locations. The offset distance 
€, = C/4 when the shear center is located at the mid-chord in part (a) is used as the reference. We use 
e; to denote the shear center/aerodynamic center offset for other shear center positions. Develop the 


equations for wing lift L and twist @ as functions of the ratio vA . The divergence speed will 
oO 


change with changes in “/, . Plot lift and twist angle vs. airspeed for ve values of 0, 0.25, 0.50 
oO oO 


and 0.75. 

Solution 

(a) Begin by solving for the wing divergence speed. We know (by combining Eqn. 2.6.18 with 
given data) that the divergence speed when the shear center is at the mid-chord is: 


d po = = 100 / 2 
t. 
Cte fe (2.7.1) 


1 
From the definition of dynamic pressure, g = 3 pv’ 


Dae 200 ee 
Vise [7 -{ =290.1 
PO Vp ¥0.0023769 Tree 


(2.7.2) 
Using Eqns. 2.6.13 and 2.6.18 the twist angle formula becomes: 
a, + [ ou \ouse 2 
oO e, 
6, -( q 2 -( ud Oa (2.7.3) 
Qe 2 [ Y } Yoo) | _{¥ 
Ido Vi, 
where the effective initial angle of attack @, is: 
GC 
a, =a, +| — || 4° |=-0.2476 (2.7.4) 
€ C,, 


As a result the twist angle will always be negative. Substituting Eqns. 2.7.2 and 2.7.4 into Eqn. 
2.7.3, the twist angle expression is 
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v \| 0.2476 
= 5 (2.7.5) 
290.1 V 
| 
Vio 
Vv? 
or 6 =-0.29421x10> —— a (2.7.6) 
a) 
Voo 


Lift is composed of two different terms, the lift due to the initial angle of attack and the twist angle 
due to both angle of attack and camber. The lift expression computed from Eqn. (2.6.14) is 


L, = or. a, + q C Cuac 
ieee 7a A eNO Oe QD) 


Substituting our data we find the lift when the shear center is at the mid-chord to be: 


SC : 
| eee osme-09( (2.7.8) 


“ileal ae be 
Vo0 


Using the expansion qSC, = 4 (4,s¢,, ) = 4 (100*10*4) we have 


oO oO 


2 
ok ok 
t-( 4 suis z aosass-o3{ (2.7.9) 
po 1 (¥ Vo0 


The final result is 


2) 
i] y 
L, =| —“ matron 2) (2.7.10) 


Do 
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(b) Examine the effect of different shear center/aerodynamic center offsets on the amount of twist 
and lift generated. When the offset between the aerodynamic center and the shear center is e; the 


divergence dynamic pressure is: 


Ky 
Ini = 
Bers (2.7.11) 
The reference divergence dynamic pressure is given as: 
K. 
Im =~ Pe = 100 psf. 
Co\Lor (2.7.12) 
The ratio of these two divergence dynamic pressures is 
Ini _ &o _ Voi 
gt 
Ino & Voo (2.7.13) 


Using Eqns. 2.7.3 and 2.7.13, the twist angle for an arbitrary shear center position is: 


G 
a: ee (2.7.14) 
Vi Aes 


The lift on the wing is 


L, = SC, @, + qgSC;, 


or (2.7.15) 


Plots of twist angle and lift are shown in Figures 2.7.2 and 2.7.3. Note that even though the 
structural twist is always negative, at low speeds wing lift increases as airspeed increases. At higher 
speeds lift decreases dramatically. Note however that this example does not allow the wing angle of 
attack to change with airspeed. 
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airfoil twist angle (radians) 


airspeed (ft/sec.) 


Figure 2.7.2 - Elastic twist angle vs. airspeed for 4 different ratios of elastic axis offset 


0 25 50 75 100 125 150 
airspeed (ft/sec.) 


Figure 2.7.3 - Lift vs. airspeed for 3 different ratios of elastic axis offset 


2.8 Using the typical section model to understand early monoplane aircraft 
structural failures 


ScC 
When e=0, the wing twist is 0 “ee . The wing lift is L=gSC,, [«, Eu | . Notice 
T 4 
that the wing lift depends on @, but @ does not depend on the lift because there is no offset between 
the aerodynamic center and the shear center. This is a special case, not often encountered in aircraft 
designs. The Langley Aerodrome was one of these exceptions. 


The first monoplane wing designs had low torsional stiffness to permit wing warping control. Their 
wing segment sections also had substantial camber, creating a large nose down (negative) pitching 
moment. The effect of the nose down pitching moment on wing twist is seen in the second term in 
Eqns. 2.6.14 and 2.7.15. When the angle of attack is zero, this initial twist induced angle of attack is 
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_ SCCyac 


ie (2.8.1) 


M 


The term Cy,,- becomes more negative as camber increases. The twist due to camber becomes 


more negative as airspeed increases; its effect on lift is amplified by aeroelastic interaction. 


Garrick and Reed" believe that this situation may explain why one of Langley's Aerodrome’s 
crashed on take-off from its houseboat launch platform. Since the Aerodrome wing camber was 
large, a large nose-down pitching moment was placed on the torsionally flexible wings. Even 
though the shear center and aerodynamic center were close together so that the offset distance e was 
nearly zero, the nose-down twist placed a downward load on the wing instead of the upward load the 
Langley expected. Thus structural flexibility, but not aeroelastic interaction, gave the first flight 
prize to the Wright Brothers. 


The first term in Eqn. 2.6.14 is the lift due to an input angle of attack. The effect of the lift curve 
slope C i is amplified by the aeroelastic term 


1 


1 Pe, 
K, 


The amplification term increases with increasing airspeed when the aerodynamic center and the 


(2.8.2) 


shear center are offset and the AC is ahead of the shear center, the usual case for subsonic aircraft. 


The wing angle of attack @ necessary to generate lift to support one-half the airplane weight W in 
level flight is 


L= + = qSC,,@,+qSC, Oy (2.8.3) 
where 
é Cy 
o> ues (2.8.4) 
K, 
Solving for the wing angle of attack @% we have 
1(W \| l-q@ 
Ly 
SeC ScC 
with q= re fe and, as before, Q, = RMAC 
K, Rr 


Equation 2.8.5 indicates that, for a wing with a symmetrical wing segment section (@,, = 0), the 
pilot will have to reduce the wing attitude as dynamic pressure increases, but the angle must be 
increased if Q < 0, the case for an unsymmetrical wing cross-section. 
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2.9 Load factor sensitivity to aeroelastic effects-fatal changes in response 

In Chapter 1 we mentioned that, in 1912, the British War Office banned Royal Flying Corps pilots 
from flying monoplanes after a series of accidents involving the Bristol Prier monoplane. Thurston 
(Reference 2) asserts that the cause of these accidents was overloading the brace wires. Even if we 
accept this conclusion, the cause of the overload may be aeroelastic. Some witnesses to these 
incidents reported in-flight explosions. Our simple model explains these incidents in which 
airplanes with small torsional stiffness could easily overload the wings when flying at high speed. 


A slight change in wing angle of attack will increase the wing load substantially when we fly at 
“high” dynamic pressure. Our typical section model defines “high” in terms of the parameter 


SeC 
qg= ace fe . When qg is near g =1 then the dynamic pressure is high. 
T 


To show this, differentiate the lift expression, Eqn. 2.8.3, with respect to the angle of attack @ to 
obtain 


AOL _ qSC,, 
0a, 1-q@ 


(2.9.1) 


L 
The aircraft load factor, n, is defined as n = W . When n=] the airplane is in straight, level flight. 


From Eqn. 2.9.1, the change in aircraft load factor with wing angle of attack a is: 


aly, = an _ qC,, 
aa, da, (WK)-q) a 


Equation 2.9.2 shows that the increased load factor due to a change in aircraft attitude from a pilot 
input will increase with increasing aeroelastic interaction (increasing values of q ). The load on an 


—  @gseC 
airplane wing flying close to the divergence speed 7p = Ke will be extremely sensitive 


to small changes in angle of attack. 


2.10 A-simple single degree of freedom (torsion only) aeroelastic model 


In the last three sections we have developed a two-degree-of-freedom analytical model to 
demonstrate how aerodynamics and structures interact on an unswept wing. Although the typical 
section model has bending (plunge) freedom as well as torsional freedom, only the twisting 
deformation interacts with the aerodynamic loads. This interaction creates bending but the bending 
does no create further aerodynamic loads. We also found that there is a nondimensional aeroelastic 
parameter that controls the interaction process. When this parameter has a value close to unity we 
find that the wing deflection becomes extremely large. We called this special point “divergence.” 


In this section and others to follow we will examine divergence and find that it is a static aeroelastic 


instability. To solve for this special point we need to understand more about the mathematical 
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features of a _ general divergence Lift 
problem and also understand the origins 
of aeroelastic stiffness related problems. — V— torsion spring 


To focus on static aeroelastic stability = = — -------- ™iti\Qgp-->~------ 
issues we will use the modified 
idealized wing, with planform, S, 


mounted on a torsional spring, as shown : : F . : 
Figure 2.10.1 - Pin mounted wing segment with torsion 


in Figure 2.10.1. The translation spring sueine rosirannt 


K), and the plunge degree of freedom 


are removed and replaced by a frictionless pin is inserted at the shear center. A torsion spring is 
inserted at the shear center to resist the aerodynamic moment. The wing has an initial angle of 


attack, Q and twist 6. 
As before, the aerodynamic moment about the shear center in Figure 2.10.1 isM,, 


M, = Le +Mac 
with 
Mac = 9cSCyac 


The wing lift, L, is 
L=qSC, (@+ 0) 


The elastic restoring torque of the torsion spring is Mg, defined as M,=K,@ in the 


counterclockwise direction. The unknown elastic twist angle, @ is again found from the equation of 
torque static equilibrium written about the shear center. Summing moments at the pinned shear 
center gives 

2X Msc = M, -M, = 0 


or > M spear =K,6 gSeC, ( Q, 7 6) ~ gScCyac =0 (2.10.1) 
positive center 
clockwise 


qSeC, a, +{£) Cue 
id e 


Solving for @, we find@= Divide both the numerator and 


K, —qSeC,, 


denominator by Ky to obtain 
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ae 
K, 


d= 2.10.2 
1 qSeC, 
K, 
where 
c\C 
a, =a,+ (=| See) (2.10.3) 
eC, 


: : —_ qSeC,, 2 : = 
Once again, the size of the aeroelastic term ¢ = K, 18 very important because, as q 


approaches unity, @ will approach infinity if the angle of attack is not reduced. Aeroelastic 
divergence occurs when the dynamic pressure, 7 , equals g, =1, where 


K 
Ip = Sec, (2.10.4) 


The wing lift is L=4qSC, 


or 
gSC, 
fae «(£5 Cue (2.10.5) 
1-g Wp )\e C,, 


Equation 2.10.5 is identical to Eqn. 2.6.14. No matter how small the angle @ is, the theoretical twist 
angle @ tends to infinity as g approaches qp. 


Stiffness is defined as the rate of change in a moment or force (in a given direction) per unit change 
in displacement. Displacements include translation such as plunge and angles such as @ The 
aerodynamic torsional stiffness is computed differentiating Eqn. 2.10.1 with respect to @ From this 


relationship we notice that there is an aerodynamic torsional stiffness term ae = —qSeC L, » and a 


OM ett ; 
structural stiffness term 6 “=K,. Aerodynamic stiffness resembles a negative torsion spring. 


Structural stiffness is always greater than zero, meaning that it resists twisting. When eé is positive, 
the aerodynamic stiffness is negative and decreases (becomes a larger negative number) with 
increasing flight dynamic pressure. 
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The sum of the two terms (K,; —qSeC,, ) 


torsion spring stiffness, K; 


<x 
K=K,|1-¥Y |=K,(1-9@))i D 
(or T T q)) is the 0“) 
Ip ® 
. : : : c 
wing aeroelastic torsional stiffness. As = 
, ’ 7 - 
G— 4p the aeroelastic torsional stiffness aa he, ~.. dynamic 
declines to zero, as shown in Figure 2.10.2. c Laas pressure, q 
At dynamic pressures larger than the 3 Ke. “as 
je) 


divergence dynamic pressure, the 
aeroelastic torsional stiffness is negative; 
any slight increase in angle of attack of the 
wing will cause it to twist without limit. 


Figure 2.10.2 - Aeroelastic torsional stiffness vs. q 


Unless the aeroelastic torsional stiffness is positive, there is no solution to the static equilibrium 
equation for this wing. A negative stiffness gives a negative twist angle for a positive angle of attack 
input, clearly an impossible answer. 


2.10.1 Twist angle and load aeroelastic amplification, the feedback process 
From 2002-2006 the Defense Advanced 


Research Projects Agency (DARPA) 2% \ a | * i 
funded a project known as the Morphing 6° a - 
Aircraft Structures Program. The ' . 


purpose of this project was to design 
multi-role aircraft whose configurations 
would drastically change in flight.'" This 
required complex linkages and flexible 


skins. One of the concerns was 
aeroelastic behavior, both static and 
dynamic. Two companies built large Figure 2.10.3 - NEXTGEN reconfigurable aircraft 


prototype. This swept wing aircraft had linkages that 


scale, half span model demonstrators that 
P moved a flexible skin into five different positions. 


were tested successfully in the Transonic 
Wind Tunnel (TDT) at NASA Langley 
Research Center at simulated altitudes up to 50,000 feet and speeds up to Mach 0.9 


16,17 
3. 


The wind tunnel tests were intended to test 
and to demonstrate the structural integrity of 
these morphing wings. This demonstration 
included the ability to predict stresses in a 
reconfigurable structure with numerous 
linkages and joints and to predict flutter 
speed. One of the two models demonstrated 
stress predictions to within 5% of stresses 


observed in the wind tunnel while the other y 

model had larger errors of the order of 10%. Figure 2.10.4 - Typical CFD pressure distribution 
detail showing color coded pressure distribution over 

What was surprising about this 10% a transport aircraft wing 
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difference was that the mathematical model used computational fluid dynamics (CFD) modeling to 
predict the external aerodynamic loads; these highly accurate loads were then incorporated into the 
structures finite element model. The model with the 5% errors used a less complicated aerodynamic 
paneling method, but demonstrated better correlation between the observed and predicted stresses. 


The CFD accuracy problem was traced to the aeroelastic modeling or rather the lack of aeroelastic 
modeling for the static loads. The feedback between the wing surface deflections and the 
aerodynamic loads was nonexistent. The CFD calculations were run for the rigid wing surface, then 
applied to the rigid structure; the process ended there. More specifically, this procedure, involving 
exquisitely accurate aerodynamic loads calculation, proceeded as if there were no aerodynamic 
stiffness matrix like that displayed in Eqn. 2.6.9. 


BG PND ges Os tate ge eae 2 (2.6.9) 
0: Kelle 0 sedhel Le GOCUMAC } We 


Of course the large finite element model had thousands of degrees of freedom, but the equation 
structure was the same. CFD aerodynamic models require exact knowledge of the wing surface 
geometry in Eqn. 2.6.9 this data is represented by the loads on the right hand side that are functions 
of the angle of attack and camber. CFD models are not efficient when used to compute aerodynamic 
influence coefficients that go into the aerodynamic stiffness matrix. This deficiency can be 
remedied by iteration, that is, computing deflections using a version of Eqn. 2.6.9 without the 
aerodynamic stiffness matrix: 


a eh =qSC, a : +qScC 
Cc 
0 K, N16 Qo, &, 4 QOC¥ vac 1 


Where the terms /1, and @, refer to the first estimates of deflections. We can then iterate to find the 
corrections to the loads by solving the equation: 


[ella eee a) teem 


Eventually the process should converge and the deflected shape will be the sum of the series of 
terms generated by the iterative process. On the other hand, there is a possibility that the sum of the 
terms will diverge. We need to examine this possibility. 


Let's return to our single degree of freedom model to examine the importance of including 
aeroelastic interaction effects. Equation 2.10.1 for the elastic twist can be written as 


a 
wg 90 (2.10.6) 


G= 
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1 
where the amplification factor, f, is j= _ which can be expanded as an infinite series: 
—q 
1 rn ae ee ‘ sn 
ere ae al aa +..=14+ 09 (2.10.7) 
—q n=1 
This infinite series convergences if g < 1, but diverges if q bar = 0.1 


q 21. This explains why the deflections are infinite 
above the value g=1. At g=1 we have a critical 
condition. 


Let's compare the size of the first 5 terms in the series in 
Eqn. 2.10.7 for different values of g to see how quickly 
the series converges and to see when aeroelastic effects 
are likely to be important. The first term in the series 
will always be equal to 1. 


Figure 2.10.5 shows that terms beyond g may be 
neglected with little inaccuracy when gq is small. 
Aeroelastic interaction adds about 10% to the 
nonaeroelastic twist angle. However, in Figure 2.10.6, 
we see that when g is equal to 0.50, terms in the series 
out to g * are important to the computation of the twist 
angle. Finally, as shown in Figure 2.10.7, the 
correction terms when gq = 0.90 are strong contributors 


to the infinite series well past the g * term in the series. 


Equations 2.10.6 and 2.10.7 are combined to get a 
series expression for elastic twist. 


SeC, a eats 
=P 45477 +...) 2.10.8) 
T 
The first term in this series is 
SeC, a 
eae (2.10.9) 


a 


The angle @, is twist angle when no aeroelastic interaction 
or feedback is included. The second term in the series is 


= Sec. 
6. = 90, =e 


2.10.10 
0 K ( ) 


Tr 


The term @, is the wing twist in response to @, (not a). 
The third term in the series is 
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Figure 2.10.5 - Series terms for twist under . 
load when q =0.1 


0.0 


series term size 


series expansion terms 
Figure 2.10.6 - Correction terms for twist 
when gq =0.5 


series term size 


1 2 3 4 5 
series expansion terms 


Figure 2.10.7 - Correction terms for twist 
when g = 0.90 
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0,=9 0,=98, (2.10.11) 


This term is the system twist in response to @. Each term in the infinite series for twist is a 
corrector for the added moment due to the twist added by the previous term. As a result, the series 
for the twist including aeroelastic effects is written as 


0=0,+>.6, (2.10.12) 
n=l 


where 9, = q'0, =q6,_, so that when @ is near 1, the additional twist added to each term is nearly 
as large as the original input twist. 

From these brief results we can see that neglecting the interaction between structural distortion and 
the aerodynamic loads on a wing surface can lead to underestimation of deformations and resulting 
stresses. This underestimation can be slight, moderate or severe depending on how close we are 
operating to the wing divergence speed. In Chapter 3 we will see that the divergence dynamic 
pressure is negative, however, it still will furnish us with a benchmark that is not to be ignored. 


2.10.2 Static stability — defining the critical dynamic pressure point 


For constant airplane weight, W, the wing supports 2 the weight so the wing angle of attack is: 


aL, -1() 1-q = gScCuac (2.10.13) 
2\ S/\ gC, K; 


The wing twist angle at constant weight is 


1 {ew 
@= rales + 48eCunc| (2.10.14) 


At the theoretical divergence condition, the wing input angle of attack in Eqn. 2.10.13 is 


Ge -( ae ) = -(<) Cuac with the wing twist angle from Eqn. 2.10.14 
Mi e Ly 


6, = a7 (5) ces . The wing twist angle is not infinite, as would be the case if we did not 
T e Es 


change the angle of attack. 


When the lift is controlled, as it would be in flight, we cannot learn anything about stability by 
looking for an airspeed or dynamic pressure at which the deformation becomes unbounded (goes to 
infinity). Only when the wing is given a fixed angle of attack will this happen. The wing is in 
equilibrium, as indicated by our equations above, but it is unstable if something moves it off this 
equilibrium point. 


Stability is the ability of a physical system to return to its original equilibrium state when disturbed 
or “perturbed.” The stability of any structural system can be “tested” by applying “Euler’s stability 
test” or “Euler’s criterion.” This is also called the energy test. Simply stated, “a system in static 
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equilibrium is statically stable if it returns to its original equilibrium position after being disturbed 
by some extraneous action.” This boundary between stable and unstable occurs at a transition point, 
in our case, the divergence dynamic pressure. 


The most familiar type of static structural instability is column buckling in which an axial force 
creates side deflection at a critical load called the buckling load. At buckling we encounter 
“bifurcation” in which there are multiple static equilibrium states in addition to the original 
equilibrium state. These states will have different structural deformation shapes. Buckling involves 
solving an eigenvalue problem in which the 
eigenvalues are related to the critical loads (in our 
case, critical dynamic pressure, gp) and the 
eigenvectors (deformed shapes) are related to the AQ>*-- 
structural shape at the critical load. 


lift + perturbation lift 


Ms= K;(6+A8) 
aniston eerne 


Consider flying at a dynamic pressure so that the 
wing angle of attack is constant, given by Eqn. 
2.10.13 and the wing twist @ is given by Eqn. 
2.10.14. We disturb (in stability jargon, we Figure 2.10.8 - Flexible wing segment 
"perturb") the wing twist slightly from its static perturbations in angle and load 
equilibrium value @ (given in Eqn. 2.10.14) to a 

new position 9 + A@, using a disturbance moment created by a force of very short duration to move 
the wing away from its equilibrium position. 


As shown in Figure 2.10.8, the aerodynamic and structural moments about the shear center change. 
These changes are computed to be: 


AM , = qSeC, (A@) 
and 


AM, = K,(A8) 


These two increments or “deltas” are usually not equal. The net restoring moment (the difference 
between the perturbation moment exerted by the spring and the perturbation moment provided by the 
wing aerodynamics) due to the change in twist is 


AM, = AM, - AM, =(K, —qSeC,, )A@ (2.10.15) 


At all values of g, except g =1, Eqn. 2.10.15 is not zero and the perturbation in the twist angle will 
create angular motion; the system will not be in static equilibrium. 


Static structural stability analysis uses Eqn. 2.10.15 as a starting point and asks whether the 
perturbed wing tends to come back to its original equilibrium twist value or whether it moves away. 
If the wing tends to move back towards its original static equilibrium position the original static 
equilibrium position is stable and unique. If the wing tends to move further away from the perturbed 
position, then the static equilibrium position is unstable. 
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There is a third possibility. If Eqn. 2.10.15 is zero when the perturbation A@ is very small, but non- 
zero, then the system also will be in static equilibrium in the perturbed state. This condition is called 
neutral stability and we are at a transition point, also called a bifurcation point (two equilibrium 
positions) or critical point. The arbitrary displacement angle A@can be factored out of Eqn. 2.10.15, 
indicating that the test for static stability does not involve the size of A@. 


In Eqn. 2.10.15, if AM, >AM_, at the wing position 9+A@ the spring torsional moment 


overpowers the aerodynamic moment and causes the wing to return to its original equilibrium 
position, @. (Actually the wing will oscillate about the equilibrium point, as we will see in another 


section in the text). On the other hand, if AM, < AM, , the aerodynamic moment overpowers the 
structural moments and the wing will move away or “diverge” from the static equilibrium point. 


If Ky is greater than gSeC, then static aeroelastic stability is assured. If we fly at a condition 
where Ky is less than gSeC Lo, We Will have torsional instability in which a small disturbance causes 


the wing to move away, uncontrollably, from its static equilibrium position. 


The crossover point between stability and instability occurs at the divergence dynamic pressure gp. 
This is our zero stiffness condition occurring when 


K, = qpSeC,, 
so that 
ee (2.10.16) 
oP Bee, - 


This result is identical to our previous results in which we found that flying above this dynamic 
pressure created infinite deflection on uncontrolled wings. 


2.11 Illustrative examples — doing the math 


Let's consider two examples to illustrate how to solve basic static stability problems. 


2.11.1 A non-aeroelastic example — static stability 
of a magneto-elastic device bending spring 


The idealized configuration shown in Figure 2.11.1 is a Kr (in-lb./rad.) 
“magneto-elastic” device. This involves a magnetic 
armature located at the end of a flexible shaft and suspended 
between two magnets. The magnetic field creates attraction 
forces on the armature; these forces depend on the distance 
between the magnets and the armature. This model replaces 


the flexible shaft with a rigid rod free to rotate about a pivot 


armature 


whose rotation @ is restrained by a bending spring that resists 


rotation @, the “bending angle” shown in the figure. 


If the armature is displaced sideways by a small angle @ the . . . 
; : Figure 2.11.1 — Magnetic device model 
magnetic forces change. The forces are given by the 
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expressions P = —— and P, = ——— 
(d+y) (d—y) 
distance when @ is zero and y=bsin@ is the lateral displacement of the end of the rigid. Only 


where @ is a field strength parameter, d is the gap 


when @is zero are the forces equal. 


Problem 
Find the value of the magnetic field parameter a that makes this configuration neutrally statically 


; d 
stable. For your calculations, restrict deflections to be small so that sind =@ and 0 << a Note 
2 


that (d- y)'(d+ yY =(42-y?). 


Solution 
Perturb the device an amount @ and then sum moments about the pivot point. Use the counter- 


clockwise direction as positive 
OM , =P; b-P, b-K, 0=0 (2.11.1) 


Substitute expressions for P,, P) (and note that y=b@) into Eqn. 2.11.1 to get following equation. 


K,0 


ym, = ab as ab a 
(d—b@) (d+bé@) 


b(4db) (2.11.2) 
a 
>”, | oe. = K, |e 
Since 0 << yy we can neglect higher order terms in @and simplify Eqn. 2.11.2. 
4ab? 
>M, = 7 —~K;, |0 (2.11.3) 


Equation 2.11.3 is the negative of the system stiffness. There are two ways to solve this problem. 
The first is by computing the value of a at which the system stiffness becomes zero. This leads us to 
the neutrally statically stable condition, written as 


aM _ 4 


—= (2.11.4) 
00 
From Egn. 2.11.3 we get 
OM 4ab* 
—=0= —K (2.11.5) 
00 
Solving for a 
K,{d° 
Aerit = f(a (2. 1 1.6) 


Copyright Terrence A. Weisshaar 
Chapter 2 


page 68 - Static aeroelasticity-structural loads and performance 


The second method of obtaining the solution is to set Eqn. 2.11.3 equal to zero. 


Aab? 
vm, -(*8- x, Jao om ule 


Ke 
The result is the same, a,,;, = Be 


2.11.2 Example —a flexible aircraft in steady turning flight 


An idealized aircraft is shown in Figure 
2.11.2 flying at an angle of attack @, with 
its idealized, uncambered "wings" twisted 
an amount @ The wing lift for each of its 
two, idealized, wing segments is 
L, =qSC,, (a, +0) where 


Q. = wing angle of attack 
O= wing elastic twist 
S = wing area (per wing segment) 


During a banked, constant altitude turn 
this airplane develops a load factor n, is 
defined as 


_ Tota 1 Lift Figure 2.11.2 - Aircraft and wing free body diagram 
Total Weight 


where Total lift = 2gSC,, (a, + 9). 


The total aircraft weight isW=(2m+M)g, where m=wing mass (per wing segment) and 
M=fuselage mass (plus everything else). With these definitions, the load factor is: 


2qSC,,(a@, + 9) 
n=— 


2.11.7 
(2m + M)g ( ) 


We will solve for the aircraft angle of attack @ and the wing segment twist angle @ as functions of 
n. We will also solve for the divergence dynamic pressure. 


Equilibrium requires that 
L=nW=2L, (2.11.8) 


In this problem, 7 is the independent variable with &, determined from Eqn. 2.11.8. 
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The restoring torque, Ms, on each wing segment is Ms = Kr@. When we sum moments (with 
L,, =qSC, (@,+@)) about either wing shear center we get the following static equilibrium 


equation: 
YM ging = 0= Ky - Lye —nmgd (2.11.9) 


K,0=L,e+nmgd (2.11.10) 


From Egns. 2.11.8 and 2.11.10 we get 


Ky8=Ly (e+ m8") (2.11.11) 
Ww 
Defining the new parameter {Z = a we have 
K,0 = Ly (e+) =(qSC,, (a, +9))(e+u) (2.11.12) 
Solving for the twist angle: 
SC, (e+ 
BCA te) (2.11.13) 


7 K,; — gSC, (e+) 


From Egn. 2.11.13 it appears that the twist angle @ becomes infinite when the denominator 

approaches zero. This implies that aircraft weight ratio (and dimension d) determines the 

divergence dynamic pressure. However, Q, is a function of dynamic pressure so we need to explore 

Egn. 2.11.13 further. 

= K, 
SeC,, 


The divergence dynamic pressure for the restrained wing (with o, held constant) is ID inet 


4: [1+ Ja, 
qd e 
6 = het (2.11.14) 


jot [1 " 4 
Os, € 


The lift force acting on each wing segment is Ly, 


Rewrite Eqn. 2.11.13 as 


(2.11.15) 


Rewriting Eqn. 2.11.15 with a common denominator, we get 
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(2.11.16) 


The angle of attack of the two wing segments, 0, is found by setting each wing segment lift equal 


to % the aircraft weight, W, so that L, = ae . Solving for @%, using Eqn. 2.11.16, we get 


a, = (2.11.17) 


Now substitute Eqn. 2.11.17 into Eqn. 2.11.14. This gives an expression for the wing segment twist 
angle @in terms of n. 


(2.11.18) 


qSC,, |1-—4 [+4] 
° Wes e 


-( i (“4 (2.11.19) 
2 )\ Ky e 


Although the twist angle @ is a function of the load factor, it is not a function of airspeed once we 
take into account that the wing angle of attack must change with airspeed. It is not possible to 
determine anything about wing static instability by simply examining the denominator in Egn. 
2.11.19. 


Equation 2.11.18 reduces to 


The reason that the static equilibrium value of @ does not become infinite as airspeed increases is 
that we are constantly reducing the angle of attack to keep the lift on the wing segments equal to the 
aircraft weight. As we fly faster, we get enough lift with reduced @ and constant @ A static 
stability analysis using perturbations from the static equilibrium state is necessary to find the 
divergence dynamic pressure. 


To examine the effects of increased dynamic pressure on torsional stiffness, we re-write the torsional 
equilibrium equation as 


K,0=qSeC,,(&% + 0)+nmgd (2.11.20) 


Combine @ terms to get 
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(K, —qSeC,, )() = qSeC,, @, +nmgd (2.11.21) 
If we perturb that wing we transform Eqn. 2.11.1 into 
(K, —qSeC,, )(0+A) = qSeC, a, +nmgd (2.11.22) 
This leaves us with two separate parts 
(K, -qSeC,, )(@)+(K, —qSeC,, )(A@) = qSeC, a, + nmgd (2.11.23) 
From Eqn. 2.11.23, we see that for two solutions for static equilibrium to exist we must have 
(K,—qSeC,,)=0 or K, =qSeC,, (2.11.24) 


The load factor does not appear in Eqn. 2.11.24 so that the divergence q for this example is the same 
as found from our fixed wing stability analysis. We conclude that inertia loads (or gravity) don’t 
change the divergence speed. 


2.12 Linear analysis vs. nonlinear analysis-divergence with a nonlinear aerodynamic 
load 


The onset of static structural instability is predicted by locating the bifurcation point at which two or 
more structural deformation forms are possible for the same external, applied loading. In this 
section we will show that nonlinear theory predicts the appearance of multiple equilibrium states, but 
with finite, not infinite, deformations. The purpose 
of this section is to show how linear system stability 
concepts are related to nonlinear system behavior. 
The problem we will solve is one for which the 
structural restoring moment is linear but the 
aerodynamic load is nonlinear. 


Our previous static equilibrium analyses assumed 
that there was a linear relationship between the 
structural and aerodynamic loads and the twist 
deflection. These analyses predicted that, at the Figure 2.12.1 - Single degree of freedom airfoil 
divergence speed, the wing twist and the with nonlinear aerodynamic load and linear 

aerodynamic load become infinite no matter how structural spring 

small the external applied load. Linear analysis also 

predicts that the aeroelastic stiffness becomes negative once the system becomes statically unstable. 


Consider the single degree of freedom uncambered wing with the cross-section shown in Figure 
2.12.1. A concentrated load P is applied aft of the shear center (represented as a pin in Figure 
2.12.1) to create a moment Pd about the shear center to create twist 6 The wing has no initial 
incidence angle. Twist,@, is resisted by a linear torsion spring, but the aerodynamic lift is a 
nonlinear function of @(measured in radians) with lift coefficient C7, given as: 
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CaG,. [ota G 121) 


The effects of @ on the moment arms e and d are ignored (in fact, these moment arms are ecos@ 
and dcos@, respectively). The torsional equilibrium equation is: 


Y Myc =0 = Le— K,0+Pd (2.12.2) 


Substituting the lift expression in Eqn. 2.12.1, the torsional equilibrium equation in terms of twist 
angle @ is: 


16 
(K, —qSeC,, )O+—-qSeC,, 6 = Pd (2.12.3) 
_ gSeC, 
Dividing all terms in Eqn. 2.12.3 by Ky (with the notation g =————“) we obtain the 
T 
nondimensional torsional equilibrium equation: 
= = Pd — 
(1-7)0+~7e =— =P (2.12.4) 
3 Kr 


en 
The parameter P = — is the elastic twist that would occur if no aerodynamic load acted on the 
T 
section. 


The effective torsional stiffness K, for this single degree of freedom system is defined as: 


OPd 
K,=— 2.12.5 
ay ( ) 
Flexibility is the inverse of stiffness: 
1 00 
exibility = — = 2.12.6 
a rage (2.12.6) 


From Eqns. 2.12.3 and 2.12.4, the effective torsional stiffness of the nonlinear system is 


K, =(K, —qSeC,, )0+16qSeC, 0” (2.12.7) 


e 


Divide the effective stiffness in Eqn. 2.12.7 by K, to obtain the nondimensional effective torsional 


Stiffness: 
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K =e =F = (9+ 1696) (2.12.8) 


The effective torsional stiffness is a 
nonlinear function of twist 0. A 
stable system will resist the applied 
torque; the effective stiffness of a 
stable system is positive. The 
effective stiffness expression in Eqn. 
2.12.8 has hardening features. As 
twist angles increase the system 
becomes more resistant to twist. 


Agar divergence 


relative effective 
torsional stiffness 


A plot Eqn. 2.12.8 at four different 
values of g: 0, 0.5, 1 and 1.5, is 
shown in Figure 2.12.2. The 
expression for effective stiffness is 
symmetrical with respect to 9. When 
@ is zero the linear and nonlinear Figure 2.12.2 - Effective torsional stiffness vs. twist angle ® 
expressions for effective stiffness are (Eqn. 2.12.8) at four different dynamic pressures. 
identical. The figure indicates that 

effective stiffness has a value of unity at zero dynamic pressure and then declines, becoming zero at 
divergence and negative at dynamic pressures above divergence. 


00 O01 O02 O3 04 05 
twist angle 


When g = 0.5 (50% of the divergence dynamic pressure parameter), the system is stable because 
the effective stiffness is positive; the effective torsional stiffness increases as the section twists and 
the effective stiffness increases (hardens). 


When g =1 (the curve marked divergence) the effective stiffness at 9= 0 is zero. The linear system 


model predicts neutral stability, divergence, at this point. However, as load increases the wing twists; 
the effective torsional stiffness becomes positive and thus stable, but not with zero deflection. This 
indicates that the section is stable when 6 # 0 and will resist further twist, even when g =1. This 


is a major difference between the nonlinear model and the linear model 


When g =1.5 (the curve marked 150% divergence) we have a negative effective torsional stiffness 


at @=(; the system is statically unstable. Figure 2.12.2 shows that the effective stiffness of the 
nonlinear system does not become positive until the section twists about 0.14 radians, a fairly large 
displacement. The linear and the nonlinear systems are unstable (the effective torsional stiffness is 
negative) for a small range of twist angles when g = 1.5, but the nonlinear system will restabilize if 
it rotates more than 0.14 radians. 


Figure 2.12.3 shows the plot of the static equilibrium twist angle as a function of applied load 


P Ye given in Eqn. 2.12.4. When there is no aerodynamic load, the relationship between load 
T 


parameter and twist angle is linear, as indicated by the line marked linear system. At dynamic 
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pressures below the divergence speed, 
for a fixed value of Pe Ye , there is 
T 


only one unique value of @ When 
operating at the divergence dynamic 
pressure, Figure 2.12.3 indicates that, 
while the solution for twist angle is 
unique, the region near 0 = 0 is very 
sensitive to small changes in applied 


torque; a small change in P Ye 
T 


nondimensional applied torque 


creates large changes in @. 


40.50 V ' 
sesaengg a, ALL TTT 


is zero indicating that the stiffness is ae = 04 03 02 01 00 01 02 03 04 O05 
zero at this point (the flexibility is Twist angle 6, radians 
infinite). As a result, large twist 
angles result from small changes in 
applied torque Pd. 


Figure 2.12.3 - A plot of static equilibrium twist angle Ovs. 


load parameter P (Eqn. 2.12.4) showing multiple 
equilibrium states when operating above divergence where 
the wing is statically unstable. 
When the system is _ neutrally 


statically stable, as predicted by linear theory, there are an infinite number of nearby, but 
undetermined, equilibrium states, Euler’s so called proximate states, at the same value of load. For 
the nonlinear system, a slight load change produces more deflection but there are no infinite 
deflections and the load deformation relationship is unique. 


From Figure 2.12.3, when g =1.5, linear system analysis predicts that the effective torsional 
stiffness of the section is always negative and, as a result, statically unstable. A positive applied 
torque produces in a negative twist - a physically unacceptable result. In fact, the nonlinear analysis 
predicts that, when a small positive torque is applied, three possible equilibrium states exist, one of 
them zero but two relatively far away from zero twist. 


The two nonzero equilibrium states shown in Figure 2.12.3 when g =1.5 are stable because the 
slopes of the curves, the effective torsional stiffness, have positive slopes. The third static 
equilibrium angle at @= 0 is unstable because the effective stiffness is negative. The conclusion is 
that, above the neutral stability point identified by linear analysis, there are multiple equilibrium 
twist angles; the equilibrium state depends not only on the equilibrium equations but on which 
direction the system is perturbed after the load is applied. While the wing will not undergo infinite 
deformation, the nonlinear analysis predicts a finite, but large deflection that is most likely 
unacceptable. 


Homework problem-nonlinear torsional spring stiffness 
The single degree of freedom wing model shown previously in Figure 2.12.1 has a linear 


aerodynamic applied load given as L=qSC, @ but it resists applied torque by developing a 


nonlinear structural restoring moment given by the relationship: 
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M, = K,0+ Ky, 


Ky, can be either positive or negative. When K7, is positive we have a so-called "hardening 


spring.’ 


' When K7, is negative the spring is called a "softening spring." When K7, is zero, we 


have the usual linear spring relationship. 


a) 


When 7! Kc Q).2 derive the nonlinear equations of static equilibrium when an initial 
T 


moment M, = qSeC,,,,@, is applied at a flight dynamic pressure q. Plot —F vs. @ for 
T 


La~o 


four different values of dynamic pressure parameter: g = 4 0, 1, 0.5, 1.0, 1.1, with 
Ip 


K 
Ip =~ for 0.5<0<0.5 


CCrg 


Answer - EO (1-g)0+0.26° 


ve 


b) The effective torsional stiffness is defined as the change in M, with respect to 0 


c) 


oM, 
Keg = 76 The normalized effective torsional stiffness is defined as 
Keg 1 OM , 


so that K, =(1-@) +0.6*@° Using the same values of dynamic 


pressure parameter as in part (a) plot normalized effective stiffness as a function of @for 


VIeGe0s. 


Remove the external load and verify that the static equilibrium equation becomes 
K. =. K. 

oe +7 (1-9) 8 =0 . This equation has three solutions @=0 or 9@=+ me , 
yj ij 


When q <qp then there is only one equilibrium point 8, =O because the term under the 


am | (note that this is the 


radical is negative. Plot these solutions as a function of O< 
T 
inverse of the term in the twist solution above and is a measure of the spring nonlinearity. 


2.13 The effect of Mach number on divergence-the match point 

Flow compressibility introduces added complexity to aerodynamic calculations. The importance of 
flow compressibility is judged by computing the airstream Mach number. Mach number is defined 
as the ratio of the airspeed divided by the local speed of sound. Flow compressibility effects are 
subdivided into four different regions: 1) the subsonic regime 0 < M < 0.8; 2) the transonic regime 
0.8 <M <1.2; 3) the supersonic regime 1.2 <M <5; and 4) the hypersonic regime M > 5. In this 
section we will briefly examine how flow compressibility affects the divergence speed in the 
subsonic flow regime. We will define and identify the importance of the “match point.” 
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Airstream compressibility begins to become important above M = 0.3. As discussed in Section 2.5, 
wing lift (and the aerodynamic coefficients) depends on the local flow velocities at points above and 
below the wing surface. The Mach number at 

some points on the wing itself may be ee 
considerably larger than the airstream ahead of /7~ maximum velocity 


the wing. 
_ Negative 

2 Positive pressure 
Consider the airfoil shown in Figure 2.13.1. ' pressure 
The static air pressure, p, at points on the airfoil 
section decreases on the upper surface so that it 7 
is lower than the freestream static pressure. The ‘Maximum pressure—minimum velocity 
static pressure measured on the lower surface Figure 2.13.1 — Typical pressure distribution on 
may also decrease, but not as much as it does on an airfoil surface. (Reference 18) 


the upper surface. The difference between the 
two pressures defines the lift coefficient. 


The static pressure, p, on the upper or lower wing surface is: 
p=p,+Ap (2.13.1) 


where p, is the airstream static pressure ahead of the wing and Ap is defined as: 


| ee: 
Oe Pa Gs =qC, (2.13.2) 
Cc » is the local pressure coefficient on the upper or lower wing surface, written as 
y2 
C, -1-(¥) (2.13.3) 


V is the local flow velocity at any point on the wing surface and V, is the wing airspeed. This 


oO 
expression is only valid for incompressible flow. At a stagnation point where the flow speed is zero, 
(for instance, on the nose of the airfoil), the C, is equal to 1, its maximum value in incompressible 


flow. At wing surface positions where the flow speeds increases V >V, ; the pressure coefficient is 


negative. 


The lift coefficient is computed by integrating the differences between the pressure distributions 
over the upper and lower wing surface 


1 trailing 
= | edge = 
} oe (om © oper Jax (2. 13.4) 


Cc, = 


When the flow is compressible, it is shown that the maximum value of C, is greater than 1; the 
values of C, are larger those computed for incompressible flow. Prandtl and Glauert* independently 


* “Science and aeronautics have suffered a severe loss through the fatal accident to Mr. Hermann Glauert on 


August 4. Mr. Glauert was walking with his brother and his three children, and stopped to watch the blowing- 
up of a tree-stump; a large piece of wood, projected nearly 100 yards, struck him on the temple and killed him 
instantly. Born in Sheffield on October 2, 1892, Mr. Glauert was educated at King Edward VII School and 
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suggested a pressure coefficient compressibility correction factor to account for flow 
compressibility (now known as the Prandtl-Glauert transformation): 


C= ae (2.13.5) 


This approximation is accurate up to about M= 0.8. Since wing aerodynamic coefficients are 
integrals of the pressure coefficient, they will be modified by the same transformation given in Eqn. 
2.13.5. The same is true of the lift curve slope. For instance, the lift curve slope for a two- 


dimensional airfoil in incompressible flow is c, =2%. The lift curve slope, corrected for 


compressibility, is c, = ae 
p 1-2 


C 
We can now use the relationship C, = Fao 5 in the equation for the divergence dynamic 
“  J1-M 
pressure so that 
1 DD 1 a 1 2 2 K. K, 1-M? 
= VS D l= a My = —— 2.13.6 
a aoa ¥) a a oe ret 


The parameter “a” is the local speed of sound which changes with altitude. The divergence Mach 
number appears on both sides of Eqn. 2.13.6. The torsional stiffness coefficient, K;7, does not 
depend on Mach number or altitude unless there is significant aerodynamic heating. 


1 
The equation on the left hand side of Eqn. 2.13.6 g= a pa’M” is called the atmosphere line. The 


airspeed computed from the left side of Eqn. 2.13.6 must match that computed for the right hand 
side. The divergence dynamic pressure on the right hand side of Eqn. 2.13.6 is written as the 
product of an incompressible term times a Mach number term. 


_K,v1-M* 


SeC, 


‘a0 


K, 
SeC,, 


=Qp,V1-M”* with gp, = 2.13.7) 


Ip 


The intersection point between the curve defined in Eqn. 2.13.7 and the atmosphere line is called the 
match point. 


The match point changes with altitude. An example of this change is plotted in Figure 2.13.2 for a 
wing with a computed divergence dynamic pressure of 250 psf in incompressible flow. Three 
different atmosphere lines are plotted in Figure 2.13.2 for three different altitudes. The intersection 


Trinity College, Cambridge.” A similar notice was given in Flight on August 9, 1934. “The British aviation 
world was shocked to read in the news of the untimely accident which caused the death of Mr. H. Glauert near 
Aldershot last Saturday. Hermann Glauert was Principal Scientific papers Officer at the Royal Aircraft 
Establishment, Farnborough, Mr. Glauert, his wife and three children were watching the Royal Engineers 
blowing up tree stumps, and a piece of wood hit Mr. Glauert on the temple. He was a well-known authority on 
aerodynamics. Among his works the best known is, perhaps, "The Elements of Aerofoil and Airscrew 
Theory." Mr. Glauert was 42 years of age.” 
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point moves to the right as altitude increases so 250 
divergence occurs at a higher Mach number at i) 
altitude. 2 = 200 pore . q 
SF 150/20, | 
Write the equation for the match point as: a a 
— 
E 2 1001 
1 c 
| 2 22 2 
Ip =Ipov1-Mp =e Mp =4repM p (2.13.8) ne 2 Se ee 
ith 4... =~ pa’ 0 : 
Le ae a 0.00 0.25 0.50 075 1.00 
Mach number 
The reference dynamic pressure g,. is a known ; ; 
function of altitude. Now, square both sides of Eqn. Figure 2.13.2 - Match point Mach number for 


di 
2.13.10 to get a quadratic expression for M te the ll aia 


divergence Mach number 
2 2 
M4 +| 422 | 2 -| 420] =0 (2.3.9) 
Vref Vref 


so that 


altitude (ft.) 


2 4 2 
{2=| ae ta val 
Dref Vref Treg (2 3 10) 


M} = ; 


0 200 400 600 800 1000 1200 
divergence velocity (ft/sec) 


Figure 2.13.3 - Altitude vs. divergence Mach number 
(da. = 250 psf.) 


Figure 2.13.3 plots divergence speed as a 
function of altitude for the example wing. Also 
shown is the result of the calculation of 
divergence speed when the effect of Mach 
number is ignored. The computed 
divergence speed is smaller if 
compressibility effects are included. 


2.14 Divergence of multi-degree of 
freedom systems 

The objective of this section is to expand 
our study of aeroelastic stability to a two 
degree of freedom system in which both 
degrees of freedom affect the result. We 
will also show that this leads to an 
eigenvalue problem. We will use a two- 
degree-of-freedom model shown in Figure 
2.14.1, with semi-span, b, and wing chord, 
c, to develop the static stability solution. Figure 2.14.1 - Two degree of freedom wing showing 
This model consists of two wing sections definition of twist angles (degrees of freedom). 
mounted on a rigid shaft that excludes 
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wing plunge but permits torsional rotation. 
Each section is uncambered and has the same 
initial angle of attack, @%. Two torsion springs 
model the wing structure resistance between 
the inboard and outboard wing sections and 
between the wind tunnel wall and the inner 
section. 


View A-A in Figure 2.14.1 shows the two 
wing twist degrees of freedom measured and 
the initial angle of attack O. The lift on each Figure 2.14.2 - Free body diagram of the wing 
wing segment acts at its aerodynamic center segments. 

and is calculated to be 


L, = SC, (@, +4) 
a (2.14.1 a,b) 
L, =qSC,, (@, +) 
Figure 2.14.2 shows a free-body diagram of each of the two wing sections, together with the external 
and internal forces and moments for a deformed position of the system in which the outer section is 


rotated more than the inboard section. The shear forces in the rod are not shown. Summing twisting 
moments about the rod support for each section in the positive, nose-up pitch, direction, we have: 


> M, =0= 2K,(@, — 6)-3K,6,+ Le (2.14.2) 


>_M, =0=-2K,( 6, -,)+ Le (2.14.3) 


Substituting for L; and L, in terms of @, 0, and @ we create the following two equations. 


5K,0, —2K,0, =qSeC,, (a, + 4) (2.14.4) 


—2K,0, + 2K,0, = qSeC,, (@, + ,) (2.14.5) 


Collecting like terms, these equations are written in matrix form as: 


Pi leaee le eased (2.14.6) 
Ty Selig oo eel oo eee ah 


Summing forces does not provide useful aeroelastic information because these equations only relate 
the lift forces to the shear reactions in the rod support. 


In Eqn. 2.14.6 we have taken the portion of the aerodynamic loads proportional to deformations 6, 
and @ from the right hand side and moved it to the left side of the equation with the structural 
stiffness matrix. As written in Eqn. 2.14.6, the matrix contribution from the aerodynamic loads is 
referred to as the aerodynamic stiffness matrix. Divide Eqn. 2.14.6 by Kr and define a 
nondimensional parameter g as 
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_ gSeC, 
= 2.14.7 
q K, ( ) 


(note that this nondimensional parameter is identical to that used for the single degree of freedom 
wing) The equations of static equilibrium become: 


ae PF ee 2.14.8 
5 3-q\\6|° 7? ot 


The 2x2 matrix in Eqn. 2.14.8 is the aeroelastic stiffness matrix. The solutions for the unknown 
static torsional deformations, 6; and @,, are 


=o = (2.14.9) 
6, A |7-q 


The term, A, is the determinant of the aeroelastic stiffness matrix. 


A=q -—7q+6 (2.14.10) 


The values for @ and @ (divided by @) are plotted as a function of the parameter g and shown in 
Figure 2.14.3. From Egn. 2.14.10 we 
see that A+0 as gq ->1 or g > 6 so outboard 


that both 6; and @ approach infinity as 
q-lorqg-6. 


The two dynamic pressure conditions, 
q=1 and g=6, are associated with 


i. ri 


i ee aes 

Fate 2 
aerodynamic) in Eqn. 2.14.10 is zero. (|| Re ea EE 
For the  single-degree-of-freedom EE eee 


example, when the single aeroelastic 0 { 2 3 4 5 6 7 
torsional stiffness term became zero, 
divergence occurred. It appears that 


large torsional deformations and loss 
of stiffness in this two wing segment 
system. These large values occur 
because the determinant of the total 
system stiffness matrix (structural plus 


panel twist/initial angle of attack 


dynamic pressure parameter, q 


the divergence condition for the two Figure 2.14.3 — Calculated twist angles 0, and @ vs. dynamic 
degree of freedom system is (we will pressure. 
prove this in Section 2.14.1): 

A=0 (2.14.11) 


Combining the panel lift force expressions in Eqns. 2.14.1 a,b with the solutions for the static 
torsional deformations, 0, and @, Eqn. 2.14.9, the lift on each segment is: 
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feo Et Te etaf S| en 


The total lift is L; +£,. Let’s say that this two degree of freedom system is required to produce a 


Ww 
constant lift so that L= L, +L, = = (this wing segment is one of two wings lifting an aircraft with 
weight W.) The wing angle of attack required to produce constant lift equal to W/2 is: 


a -(5) (7-1)(4-6) (2.14.13) 
° (2 J3qSC,, (4-9) 


Combining Eqns. 2.14.12 and 2.14.13 we find the panel section lift forces: 


(atats) om 


The lift on each section changes with dynamic pressure. As dynamic pressure increases the 
outboard panel carries a larger fraction of the total weight. The angle of attack becomes zero at 
q =1 and. On the other hand, the lift on each segment is finite at g =1 and g =6. It does become 


infinite at ¢=4, but as we shall see, this outside the operational envelope because the wing is 
statically unstable at this point. 


2.14.1 Static aeroelastic stability of multi-degree-of-freedom systems 


The two-degree-of-freedom aeroelastic wing static equilibrium equations are written in general as: 


[lks]-al4.I1]9 | 16 (2.14.15) 


0)... 8, 
where are input loads. The vector 
Q, A; 


| is the wing deformed system equilibrium state. 


A statically stable linear system has a unique equilibrium state. But, when a system is in a state of 
neutral stability, more than one equilibrium state exists. This is the so-called bifurcation point. Let's 
see what special set of conditions must be present to obtain not one, but two solutions to Eqn. 
2.14.15. 


{| _ ja | {60, 06,| . ee a | 
Define = + where is the vector of perturbation twist angles and iS 
0, C 00, 00, os? 


the vector of original equilibrium twist angles. Use the notation [Ki ] = [K 5 | = q[K A to denote the 


aeroelastic stiffness matrix. Then Eqn. 2.14.15 becomes 
[K, }{4°}+[X, }{68,}={0,} (2.14.16) 
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By definition {Q,}=[ K, |e} so we have the following requirement for the existence of the 


perturbed equilibrium state. 


(x, }f98,}=10,}-[&, lo} (0 aun 
The condition for the existence of the perturbed static equilibrium state is: 
[ K; |{69,} = {0} (2.14.18) 


Equation 2.14.18 is simply the wing static equilibrium equation, Eqn. 2.14.15, with all external, non- 
deformation dependent loads removed, set to zero. As a result, static stability analysis is simply a 
by-product of our usual static equilibrium analysis. 


When we drop the notation 66 in Eqn. 2.4.18 we can write this equation as: 


He beef oof ltr Tet easy 


with the understanding that the displacements @ are perturbations, not real displacements. In these 
equations we clearly identify the critical dynamic pressure as an eigenvalue and the displacement 


6 
vector : as an eigenvector. 
2 


For our example the eigenvalue problem is, from Eqn. 2.14.6: 
gSeC, \\-1 0 }JG]| | 5 -2)/4 
K, J|0 -1}|af |-2 2 |]6, 
Aft. O42, > 2/12 
q = (2.14.20) 
0 -1)|4 —2 2 || 
From Eqn. 2.14.8 this eigenvalue problem can also be written as: 


7 ed Raed | (2.14.21) 
le) oS he =a 


2.14.2 Divergence eigenvalues 
One obvious solution to Eqns. 2.14.18, 2.14.20, 2.14.21 is {60,}={0}, but this solution simply 


defines the original, undisturbed equilibrium state since it means that we have no perturbation of the 


or 


system; {66,}={0} is called the trivial solution. As mentioned previously, the existence of a non- 


zero perturbation solution is known as Euler's criteria for neutral static stability. The neutral 
stability condition involves finding a set of self-equilibrating deformations that do not require 
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external inputs such as @. At neutral stability the system is able to remain in a specially shaped 
deformed state without external loading. 


Linear algebra tells us that the perturbations 60, and 6@,can be nonzero only if the determinant of 


the aeroelastic stiffness matrix| K;, | is zero. 


=A=0 (2.14.22) 


Equation 2.14.22 is identical to the neutral stability condition we found in Eqn. 2.14.11; A is called 
the stability determinant. It provides us with the characteristic equation for stability, a polynomial 
function of dynamic pressure. 


Dynamic pressure is an eigenvalue of the aeroelastic stiffness matrix. When gq is equal to an 
eigenvalue, two static equilibrium states will be possible. This special value of dynamic pressure is 
called a bifurcation point; Eqn. 2.14.22 is 
a formal mathematical condition for 
neutral static stability of a linear, multi- 
degree-of-freedom aeroelastic system. 


In general, a structural system with n 
independent displacement degrees of 
freedom (the present example has n=2) 
has n eigenvalues. The eigenvalue that 
yields the lowest positive value of q 
usually has the most physical importance. 
(Note that there might be negative 
eigenvalues or even imaginary eigenvalues 


aeroelastic stiffness determinant 


which would not make sense physically). 0 2 4 6 8 

For our example, the eigenvalues g, and dynamic pressure parameter, q 

q for Eqn. 2.14.16 are g, =1 and g, =6. Figure 2.14.4- Aeroelastic stiffness matrix determinant 
Avs. q. 


The strain energy stored in the springs is 

always positive so the determinant of the aeroelastic stiffness matrix is positive when the dynamic 
pressure is zero. As shown in Figure 2.14.4, the determinant decreases as the dynamic pressure 
increases. 


At g =1 the total energy stored as the result of wing torsion is just equal to the work done for some 


specific combinations of panel deformations (eigenvectors). At this condition, we lose the ability to 
resist additional aerodynamic moments. Further increases in airspeed will cause the determinant to 
become negative, alerting us to the fact that static equilibrium is impossible and that any additional 
loads will cause motion. This motion is divergent motion that we will examine later. 


At values of dynamic pressure above g =6 the determinant is again positive, indicating that the 


system appears to be able to store energy. While this may be so theoretically, the result is 
meaningless because divergence has already occurred and we have moved well beyond the initial 
equilibrium position. 
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All solutions to Eqns. 2.14.11 and 2.14.22 are associated with neutral static equilibrium. On the 
other hand, from Figure 2.14.4, some solutions to Eqn. 2.14.11 are associated with entry into an 
instable region while others are associated with exit from the instability. The test of entry or exit 
from static stability is whether or not the determinant becomes positive or negative when the 
dynamic pressure increases. 


The plot of the stability determinant in Figure 2.14.4 shows that the divergence dynamic pressure is 
the lower of these two values. 


Ip = 4 = (2.14.23) 


Note that if we change the value of one of the two torsion spring stiffnesses, say from 3K, to 5K,., 
the eigenvalue in Eqn. 2.14.23 will change. 


2.14.3 Divergence eigenvectors 


Wing divergence eigenvectors are found by substituting, one at a time, the values of divergence 
dynamic pressure g, into the equilibrium equations in Eqn. 2.14.18 and then solving for 60, and 


60,. We will eliminate the deltas from the 60, and 60, notation in the discussion that follows to 
simplify notation. When qg=1, the two homogeneous equilibrium equations in Eqn. 2.14.18 
become 

46-26, =0 (2.14.24) 

-26, + 6,=0 (2.14.25) 


When g =1, we find from either Eqn. 2.14.24 or 2.14.25 that 


6, = of (2.14.26) 


Notice that any twist displacement vector with the ratio ws =2 satisfies Eqns. 2.14.24 and 2.14.24. 
1 


We arbitrarily set 8, =1 and write the expression for the divergence "mode shape" associated with 


lal luo} 
a, 1.0 (2.14.27) 


Equation 2.14.27 also can be expressed as 


g)° 0.5 
af ~V10 
2 (2.14.28) 


Since any arbitrary constant C; either positive or negative, can be used and still have the wing twist 
angles satisfy the perturbation equilibrium equations in Eqns. 2.14.24 and 2.14.25. Equation 2.14.27 
(or Eqn. 2.14.28) is called a divergence mode shape corresponding to the eigenvalue gq, =1. 


eigenvalue g, =1 as the vector 
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With g =6, the perturbed static equilibrium equations become 


-9,-28, =0 (2.14.29) 
-26, —46, =0 (2.14.30) 

The solution is 
6, =-286, (2.14.31) 


With @, = 1, the mode shape (eigenvector) for ¢g =6 is written as 


a)? J-2 9 
a, = : or C, ; (2.14.32) 


A 
Look again at Figure 2.14.3, where the static response Ve is plotted as a function of g, When 


6, 
Qa, 
q — 1, the ratio “Yh tends to 0.5; this is the first eigenvector. Similarly, as ¢ — 6, the ratio % 
2 2 
is negative and approaches -2.0, the second eigenvector. 


Let's compare the strain energy stored when the system is deformed in the divergence mode shapes. 
The strain energy, U, stored by torsion springs is: 


U == (3K,6? +2K; (0,-8)') (2.14.33) 


When deformed into the first divergence mode (with C; = 1), the strain energy is: 


3 1 1 a 
U,= 3 K, (+ +K, ) 0; = ; K,0; (2.14.34) 


When deformed in the second divergence mode (with C2 = 1), the strain energy is: 
3 
U;= Sx, (4)+K, ()| 0; =15K,0; (2.14.35) 


The strain energy stored in divergence mode 2 is 24 times larger than that in divergence mode | 
because it takes more work by aerodynamic forces to deform the system into mode 2; this requires a 
much larger divergence dynamic pressure to generate these forces. 


2.14.4 Summary - divergence of multi-degree of freedom systems 


To construct the eigenvalue problem for static aeroelastic stability we first assume that the flexible 
structure is in a deformed equilibrium state. We then develop static equilibrium equations in terms 
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of the system deformation. Then, all terms on the right hand side of the equation (for instance, those 


depending on @) are set to zero. This gives us the static equilibrium equations for the neutrally 
stable perturbation states in terms of perturbation displacements. 


We can either directly expand the determinant of the resulting matrix equation in terms of 
parameters such as g, Kr and then solve for gp or, if we know the numerical values of K; and other 
problem parameters, we can let the computer find numerical values of qp. 


In general, the eigenvalue problem for static stability of an n” order wing system is: 


[K {9} =q[K, KO} (2.14.36) 
This equation has the form 
[A]{x}= ALB] x} (2.14.37) 


The eigenvalues, 4; are related to values of dynamic pressure for which the system has neutral static 
stability. 


2.15 Example - Two degree of freedom wing divergence 


A side view of a two segment wing idealization is shown in Figure 2.15.1. These segments are 
connected by two springs, one with spring constant k /b/inch and the other with spring constant 2k. 
Each wing segment has planform area S. 


F, =2kb0,— 


Biggee Sead endo wine geeneey Figure 2.15.2 - Free body diagram of configuration 


Each spring develops internal forces in response to relative deflection between its ends. The lift on 
each of the wing segments is given by 


L, =qSC,,(9, + @,) 


2.15.1 
L, =qSC,_(8, + @,) ee 


where @, is an initial angle of attack common to both segments and @, and @, are additional angles 


of incidence caused by lift on each wing segment. The free-body diagram for this configuration is 
shown in Figure 2.15.2. 
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Moment equilibrium gives the following (assuming small angles) 


Sum clockwise about pin 1 


>i M, =0=Le-kb’ (6, +6,) (2.15.2) 
Sum clockwise about pin 2 
>) M, =0= Le—kb’ (6, + 0, )—2kb’6, (2.15.3) 


Substituting L, and L, from Egns. 2.15.1 and 2.15.2 into Eqns. 2.15.3 and 2.15.4, we get 


kb°@, + kb°@, — gSeC, 8, = qSeC, @, 


(2.15.4) 
kb*0, + 3kb*@, — qSeC,, 0, = qSeC, & 


Equation 2.15.4 is written in matrix form as 


eee | OO eee adit iss) 
e = e eLO. 
1 3/(6, aa ee 0, fae 


Using the nondimensional pressure parameter 


gSeC, 


a (2.15.6) 


q = 
Then Egn. 2.15.5 becomes 


ae ee ae ee (2.15.7) 
i Biel Or el mn 


Checking the stiffness matrix 
Before continuing, let’s use the strain energy method discussed in Section 2.3 to check to see if our 
structural stiffness matrix is correct. First we form the strain energy expression. 


U == kb" (6,+0,) + 2kb*(8,) 


We then differentiate this result, using Eqn. 2.3.7 as the guide. 


2 2 
S028 (8.+0,)()=K(8.48,) may? SU 
1 


and 
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OU 2 a? (6, 40,) (1) + 2b? (0,) = kG, + 3kb? (8,) 
a0, 2 2 


o-U o-U 
= kb” 7 = 3kb* 

00,00 00: 

Comparing these expressions with the structural stiffness matrix in Eqn. 2.15.5, we see that the 

stiffness matrix matches. 


Note also that we have developed the stiffness matrix from the equilibrium equations in Egns. 2.15.2 
and 2.15.3. In these equations the structural stiffness terms are negative. We chose to take these 
terms to the left side of the equations so that they became positive. No matter how we solve these 
equations we will get the correct result. However, it is always better to organize the equilibrium 
equations so the diagonal elements of the structural stiffness matrix are positive. Then the 
equilibrium equation results will match the energy method results. 


The divergence dynamic pressure 
The divergence dynamic pressure is found by writing the homogeneous static equilibrium equations 


using Eqn. 2.15.7 in the following eigenvalue form (note again that the displacements in Eqn. 2.15.8 
are perturbation displacements): 


gi At jt 14)4 (2.15.8) 
Talli alla, ~ 


Equation 2.15.9 can also be written as 


¢ - ea 7 tof (2.15.9) 


Calculating the determinant of the aeroelastic stiffness matrix in Eqn. 2.15.9 and setting it equal to 
zero, we have 
A=q -49+2=0 (2.15.10) 


Equation 2.15.10 has two roots. Gp, =0.5858 and = gp, =3.4142 


The stability determinant in Eqn. 2.15.9 is positive between g = 0 and gp = 0.5858 but becomes 
negative at larger values. This means that the divergence dynamic pressure is: 


kb* 


eu, 


Ip = 0.5858— (2.15.10) 


2.16 Example — the shear center, center of twist and divergence 

A wind tunnel test article consists of a thin, relatively rigid shell supported on a relatively flexible 
support mechanism consisting of three beams, as shown in Figures 2.16. As a result, the wing 
surface can only rotate with respect to the wind tunnel airstream and move upward or downward. 
The structural stiffness effects of these flexible beams are idealized to be three springs, K;, K2 and 
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support frame 


MK 


ta >! ‘ 
' la abe, 


rigid wing shell | Wind tunnel wall 
al 
planform view Figure 2.16.2 — Two-dimensional idealization 
Figure 2.16.1 — Example-wind tunnel model Fy 7 
with frame support ) ™\1=Px 
xo! 
x no oo 


K; as indicated in Figure 2.16.2. In our case, 
measurements indicate that the support beams — Kk A sexe 
are of equal size so that the three spring center 
constants are all equal to K _ (lb/in.) ( F,=K(h+aé) F.=K(h) Fy=K(h-b6) 
K,=K,=K,=K). RW MW MAA ANWW Wi 


a—— 


The upward or downward wing displacement 
is represented by the displacement, h, at the 
center beam (positive downward). The wing 
rotation is due to wing elastic twist, @ as 
indicated in Figure 2.16.2. 


Figure 2.16.3 — Free body diagram showing spring 
forces and shear center position. 


A free body diagram showing the deflected springs is shown in Figure 2.16.3. This free body 
diagram shows the deflections due to h and @ create restoring forces and a moment about the 
arbitrary reference point we have chosen at the center spring. Also shown in Figure 2.16.3 is the 
unknown position of the shear center. We cannot at present use this position for our reference point 
because we don’t know where it is. Finding this position is one of our objectives. 


The structural stiffness matrix of this model is developed using strain energy methods. From the 
deflections shown in Figure 2.16.3 we conclude that the strain energy due to deformation of the three 
springs is: 


U => K(h+a8) + 5K (hy + 5K (h-08) 


2 


Using our relationships for the elements of the structural stiffness matrix (differentiating with respect 
to h and @) we find 


Kix, Bee. 3K K(a-b) 
K,, Ky | |K(a-b) K(a’+b’) 


When a downward force F (downward forces are positive, upward forces are negative because of the 
fact that we chose, arbitrarily, the downward displacement h to be positive) and a nose-up (measured 
at the reference point at the center beam and positive counter-clockwise) torque T are applied to this 
idealized wing model, the relationship between F, T, h and @is represented as follows: 
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F Ki, Ky |{h 
= (2.16.1) 
‘i Ki Bex | 
Solving Egn. 2.16.1 for h and 4, 
aE _ 
Aer K(a’+b’) K(b-a) ("| (2.16.2) 
0) Al K(b-a) 3K - 


where A= |K | =2K°’ (a +ab +b") is the determinant of the stiffness matrix. Note that since there 


are no aerodynamic forces (yet) this determinant is always positive. 


2.16.1 Locating the model shear center 


The wing sectional shear center is defined as a point on the wing section where a concentrated force 
may be applied without creating rotation. Suppose that the shear center is located at a distance x to 
the left of the center spring as indicated in Figure 2.16.3 (if we find that the shear center is to the 
right then x will be a negative number). When we apply the downward force P at the position x to 
the left of the center spring in Figure 2.16.3 we create external loads F=P andT = Px . Equation 


2.16.2 becomes: 
h K. -K P 
=1| 42 “ (2.16.3) 
0 Al “Ka Px 


P 
From this equation we see that 0 = oe (-K gtk ao) and that this twist angle is a function of both P 


and x. The size of the load is known because we are free to choose it. The distance x is unknown. 
If the value of x is the shear center then 9@=0. Using Equation 2.16.3, we find the relationship 
between P, X and the stiffness matrix elements when P is placed at the shear center. 


P 
9 =F (Kn t Kuk) =0 (2.16.4) 
Solving for x we find 
=* Koy R@=KD ad 
x= = = (2.16.5) 
Ki 3K 3 


Note that the size of the load P does not affect the shear center position. Although the rotation is 
zero, the displacement, h, is not zero, but depends on the size of P. 


eae +ab+b) = 2.166) 
3\ A 3K 
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2.16.2 Locating the center of twist 


The center of twist is defined as the point about which the wing section appears to rotate when only 
a torque is applied. To find the center of twist, we apply only a torque 7. From Eqn. 2.16.2, the 
torque will cause both twist @ and deflection h: 


_ K(6-4), 


h 2.16.7 
A ( ) 
and 
g = 2K + (2.16.8) 
A 


The downward deflection of a point y a located a distance x (to the left of the center spring) is: 


yhtxb="(b-a+ RT 
A (2.16.9) 


We set y = 0, to find the point where the wing displacement is zero due to an applied torque: 


a-—b 


2.16.10 
; ( ) 


x= 


Comparing Eqn. 2.16.10 with Eqn. 2.16.5, we see that the shear center and center of twist are located 
at the same point. 


2.16.3 Wing divergence 
To find the divergence dynamic pressure for this model, we apply an aerodynamic force 
L=qSC, a@, + qSC,,@ at the aerodynamic center located at a distance d to the right of the middle 


spring. In this case, the force and moment are F = -L and T=Ld. The aerodynamic force is negative 
because we took the positive displacement direction downward and the lift force is directed upward. 
Similarly the lift causes a positive counter-clockwise torsional moment. The static equilibrium 


matrix equation is: 
-L K, «K h 
=| (2.16.11) 
Ld K,, Ky ||@ 


The static equilibrium equation becomes 


1 Ki, (Ki. +qSC, ) | {hp 
gSC, &, = 2 (2.16.12a) 
d Ki, (Ky —qSC,,d) 6 


or 


(2.16.12b) 
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The notation K ,, 1S applied to only the two stiffness terms with q in them. The deflections are: 


h SC, a|K —~K -1 
oP a) RO Big (2.16.13) 
0 A —-K,, K,, d 


where A= K,,K5) —K,.Kj>- 
To solve for divergence we set the determinant equal to zero. 

A=K,Ky, - Ky Ky =0 (2.16.14) 
Expanding the determinant, 


K,, Ky — GSC, K,,d- Ky 7 FSC, Kis =0 


(2.16.15) 
We find that the divergence dynamic pressure is 
2. 
je) 
eel ee, 
‘ SCia d+ Ky 
mee (2.16.16) 
Stee ae Ki, 
The shear center location is ¥ =——, so that 
11 
= Ky — KX 1 
Ip Ane SC (2.16.17) 


Notice that when d =—x (aerodynamic center coincident with the shear center), divergence will not 
occur. Substituting for x, K,, and K,,, we have: 


Ki Ko- Ky 3K'(a +b')=K*(@=b)" 


SC, = 
a me ae ee a 3Kd + K(a—b) 
Pea ap ab 
or G22 = (2.16.18) 
SC, [a+ 
a 3 
a-b 


Notice that the distance e=d + is the distance between the shear center and the aerodynamic 


center, so that Eqn. 2.16.17 is similar to the expression we found from the analysis of our typical 
section, single degree of freedom system. 
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2.17 Example - Divergence of a three-degree-of-freedom system 


A wind tunnel model shown in Figure 2.17.1 consists of three rigid wing sections mounted to a shaft 
so that they are free to rotate about the shaft. The gaps in the figure are accentuated for illustrative 
purposes. Three torsion springs restrain relative motion of the wing sections with respect to each 


support wind tunnel wall 


wer 


3K,0, 


qSeC_ .(0,+8; 1 1 


aerodynamic 
ee O,+3 
K7(63-85) 


A A views A-A Kr(@3-62) 
a, \ 3 


qSeC | ,(0..+63) 


qSeC 1(049+82) 


Figure 2.17.1 - Three degree of freedom wing 
planform and side view 


other and the wind tunnel wall. The spring 
connecting wing section 1 to the wind tunnel 
wall has spring constant 3K; Ib-in/radian. The 
other springs have spring constants 2K;y and Kz, 
respectively. Each wing section has equal 
aerodynamic and geometric parameters given as S, e and C Ry 


Figure 2.17.2 - Three degree of freedom free-body 
diagrams 


The wind tunnel is operating at a simulated altitude of 35,000 feet when the system is given an initial 
angle of attack @, at airspeed, V. The angles of incidence of the three wing sections are shown in 


Figure 2.17.2. 


Problem 
1) Find the divergence airspeed assuming incompressible flow; 2) find the divergence speed in 


compressible flow. We are given that x =1500 Ib/ fr. 


Solution 
First develop the equations of static equilibrium in matrix form and then find the polynomial 


characteristic equation to calculate the divergence dynamic pressure when = =150 Ibl fe , where 


C_,, 1s the incompressible flow lift curve slope for each of the three sections, uncorrected for Mach 


number. 
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Summing torsional moments for each of the three wing sections shown in Figure 2.17.2 produces the 
following equations. 

<M ,=0=Le-3K,0,+2K,(0, -9,) 

EM, =0=L,e-2K,(0, -0,)+ K,(0,-9,) 

EM ,=0=L,e-K,(0,-9,) 


These three moment equilibrium equation, written in matrix form, become 


(SK,—qSeC, ) —2K 0 6 I 
a, 1 
—2K, (3K,- qSeC;, ) -K, 0,-= qSeC, a,\1 (2.17.1) 
0 -K, (K,-qSeC,, )| | 
- q gSeC,, 
Arbitrarily define a nondimensional dynamic pressure 4 = = = Eqn. 2.17.1 
2K, 2K, 
SeC i 
becomes 
5-2q -2 0 6, 1 
-2 3-2g -1 |)6,}=2qa,41 (2.17.2a) 
0 = MSO g es 1 
Setting the external loads to zero the divergence eigenvalue problem is 
0, 7 Sele |e, 
q4\9,7= ms 2) 3) Se, (2.17.2b) 
0, Oo eT Ee, 


The divergence dynamic pressure for incompressible flow is found by forming the determinant of 
the aeroelastic stiffness matrix in Eqn. 2.17.2a and then setting it to zero to solve for gpo. This gives 
the following characteristic equation. 


A=-8q° +36q° —36q+6=0 (2.17.3a) 
Dividing by 2 

4q° -18q° +18g-3=0 (2.17.3b) 
The three roots of the characteristic equation, Eqn. 2.17b are 


g, =9.2079 gq, =1.147 = g, =3.145 The same results are found from a MATLAB eigenvalue 
calculation using Eqn. 2.17.2b.. Divergence occurs at the lowest value of 9 so the divergence 


dynamic pressure is 


Gp, = 9.2079 so that qy, =0.4158 


ec, 


Copyright Terrence A. Weisshaar 
Chapter 2 


Static aeroelasticity (4"" edition) - page 95 


When the wind tunnel is operating at a simulated altitude of 35,000 feet, the speed of sound and air 


lb—sec? 
(me oe 


t 
density are, respectively, @,, = 973-2 p=7.365x1 The incompressible flow 
s 


2 : 2 
divergence speed V, is Vv? =“ so that ve — af aie - :| ai =1129 ie ay, F 
Pp 7.365x10™ ) SeC, SeC, sec 


so that the airspeed predicted with incompressible flow data is 


K, 


/, 1 
V, 336 = 33.6(1500)”? =1301 ft/sec (2.17.4) 


CU, 


The divergence Mach number is 


V 
M,=— or M,, =1.34 
a 


This Mach number is supersonic; the assumption of incompressibility is not valid. 
With compressibility included, the divergence dynamic pressure g, is found by using the Prandtl- 


Glauert transformation C, = = ; and then solving Eqn. 2.13.8 


Vv1-M 


_| 2K 1 
Ip =I) =, : V1-Mp => pa.M> (2.17.5) 


Squaring both sides of Eqn. 2.17.5 gives 


ns (1-M2)=M4 (2.17.6) 
pa, Sec, 


This equation is rewritten as 
M3 +A°K’M;- A’K’ = 0 (2.17.7) 


Aq 
with A -( 4D, ) K -| a } The two solutions for M,, in Egn. 2.17.7 are 


SeC, 
~A?K? +VA'K* + 44°? 
2 (2.17.8) 
One of these values is negative. The positive value of M Fs is: 


2 _ —A’K? 4+ V AK 4 442K? 


My 
2 (2.17.9) 


M,, = 
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Taking the positive square root of Eqn. 2.17.9 (Mach number must be positive) we find the 
expression for the divergence Mach number at the match point. 


eae ies 
M, =AK ae 


2 (2.17.10) 


With the values of Ip, , P.and a, inserted into A and K we have: 


2 
K SeC 
M , =0.84x10°| —*— | -1+ beaanxo'[ 


catia T 
(2.17.11) 
The divergence velocity V, , with flow compressibility, is 
V.=4,.Mp (2.17.12) 
1 2 |2 
so that V, =0.817(1500)| -1+,/1+2.83x10° (a) (2.17.13) 


When we divide V. by V,, we get 


1 22: 
Ye. = 0,02439(1500)2 —1+,/-1+2.83x10° (a5) = 0.67 (2.17.14) 
V 1500 


oO 


so V, = (0.67) (1301) =871.7 fps. and M, = 90.898. This result is transonic and a bit larger than 


recommended for the validity of the Prandtl-Glauert transformation, but shows that divergence for 
this configuration is not supersonic. 


The divergence mode shapes are found using 2.17.2b with each of the eigenvalues, taken one at a 
time, to solve the following equation. 


3-24), —2 0 6, 0 
-2 3-29, -1 |,6,7=40 
0 -1 1-29, |(@J [0 


For the first eigenvalue, gp, = 0.208 , these equations become: 
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4.584 —2 0 A, 0 
—2 2.584 -l |,6,7= 0 (2.17.15) 
0 -1 0.584 || 8, 0 


Because the 3 x 3 matrix is singular (we made it so when we formed the characteristic equation) we 
have only two independent equations. If we arbitrarily let @,= 1, the divergence mode shape for the 


wing is 


6,|\ {0.2549 
6, $=40.5842 (2.17.16) 
6, 1 


Note that we could have arbitrarily set any of the three displacements equal to | or any other number 
we might choose. 


2.18 Example - Divergence of 2 DOF offset wing segment sections 


In this section we will work a problem three different ways to show that how we choose the 
coordinate system and develop equations may change parts of the analysis, but does not change the 
result. A swept flexible lifting surface is modeled as the two interconnected, uncambered wing 
sections shown in Figure 2.18.1. This model simulates wing sweep by positioning the outer wing 
segment at a distance tan A aft of the inner section (here the term / is a dimension, not the plunge 
displacement). The relative rotation of these sections with respect to each other and with respect to 
the wing root connection is restrained by two linear elastic torsion springs with equal spring 
constants K7. These springs develop a restoring torque proportional to the twist angles @, and @, as 
shown in Figure 2.18.1. 


The elastic axis of each section is located a distance e downstream of the section aerodynamic 
center (the 4 chord line). Each wing section has a planform area equal to S$ =hc. The lift curve 


| v 


slope for each section is equal to C Tee 


offset=h tanA 


A+ 
N 
torsion Y A-« 
springs ra « h > 
gaps exaggerated nee “<r. 
for illustration oo es 
2 
section A-A 


Figure 2.18.1 - Staggered wing section model 
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Problem 
(a) Develop the equations necessary to determine the divergence speed for this configuration. 


(b) Determine the divergence dynamic pressure gy when a =10 Plot gp vs. angle A. 


Solution 

Begin by finding the equations for the twist deflection of each section. Because divergence is a self- 
equilibrating condition, we examine the system with no initial angle of attack so that the twist angles 
are perturbations from system static equilibrium. The Free Body Diagram for this problem is shown 
in Figure 2.18.2. 


Figure 2.18.2 - Free body diagrams for inboard panel (segment 1) and 
outboard panel (segment 2) 


Static torque equilibrium and shear equilibrium are written at and about the torsion springs as 
follows. 


Section 2- Equilibrium at/about the outer rod (elastic axis) 


Force equilibrium R,=L, (2.18.1a) 
Torque equilibrium (+ c-clockwise) Le—-K,@, =0 (2.18.1b) 
Section 1- Equilibrium about/at the inner rod (elastic axis) 
Force equilibrium -KR,+L,+R, =0=-RK,+L4+L, (2.18.2a) 
Torque equilibrium 
(+ c-clockwise) Le-R, (htanA)+K,8, — K,@,=0 
Substituting Eqn. 2.18.1a we get Le-L,(htanA)+ K,6, — K,@, =0 (2.18.2b) 


Note that Eqns. 2.18.1b and 2.18.2b are elastically coupled and the stiffness matrix will be 
unsymmetrical. We could remove this coupling and asymmetry by using Eqn. 2.18.1b to change 
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Eqn. 2.18.2b to Le—L,(htanA)+L,e—K,6,=0. However, we will use Eqns. 2.18.1b and 
2.18.2b as our building blocks for this part of the example. 


Note that the shear force equilibrium equations for R; and R> are required because the shear forces 


enter into the torque equilibrium equations. Substituting for the wing segment lift in terms of twist 
angles @ and @, in Eqn. 2.18.2b becomes 


gSeC,, 6, — qShtan AC, (6, +0,)—K,0, + K,9 =0 


(2.18.3) 
or K,6—-K,@, —gSeC, 6, + gShtan AC, (6, + 6, ) =0 
While Eqn. 2.18.1b becomes 
gSeC, (0+0,)—K,0@, =0 
Bm 1 5 T 2 (2.18.4) 


or K,0,—qSeC, (0,+8,)=0 


ec. 
Define the nondimensional parameter g = “2 and divide Eqns. 2.18.3 and 2.18.4 by K, to 
T 
get 
yee) _h 
1—q+q-—tanA |@ —| -1+q—tanA |@, =0 (2.18.5) 
e e 
—G9, +(1-4)@ =0 (2.18.6) 


Equations 2.18.5 and 2.18.6 are written in matrix form as 


[1-7 +a4 tana) [-1+74 tana) (e 
e e 0 =0 
2 


(-7) (l-q) (2.18.7) 


h 
When — = 10 the divergence dynamic pressure is found by forming the determinant of aeroelastic 
e 


stiffness matrix in Eqn. 2.18.7 and setting it to zero to get the characteristic equation for divergence. 
Gp + Jp (10 tan A—3)+1=0 (2.18.8) 
The divergence dynamic pressure from this polynomial is 


(3-10tan A+ /100tan? A —60tan A+5) 


Gp = , (2.18.9) 
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The plot of g, versus A is shown in Figure 
2.18.3. Notice that there are 2 solutions to 
Eqn. 2.18.9 and that these two roots bound 
an instability region for which the 
characteristic equation is negative, 


Gp + Jp (10tan A—3)+1<0. 


dynamic pressure parameter 


When the term under the radical in Egn. 
2.18.9 is negative, the divergence dynamic 
pressure is a complex number; there will be 
no real divergence speed. These values 


occur to the right of the value A=5.71° ; ; _ . 
(tanA=1/10). Figure 2.18.3 - Divergence parameter q as a function of 


offset angle (degrees) 


angle A 


When A=5.71° the divergence dynamic 

pressure parameters found from Eqn. 2.18.9 are equal; above A=5.71°, no real divergence speed 
exists. At A=26.57° the term under the radical is again zero. When A>26.57° Eqn. 2.18.9 gives 
real, but negative, divergence dynamic pressure parameter values that are real, but negative. The 
plot of Eqn. 2.18.9 is entirely in the lower, negative region for these sweep angles. 


The physical explanation for these “divergence impossible” values of sweep is that, when A>5.71° , 
the twist from the outer wing segment twists the inner wing segment nose-down so much that the 
configuration will not diverge. 


2.18.1 Solution using a different coordinate system 

Changing the wing segment coordinate system changes the equilibrium equations but cannot change 
the divergence speed. Let's illustrate this by considering the same example with a slightly different 
wing segment twist definition. For the outer wing segment (wing segment 2), the angle of attack is 
re-defined as shown in Figure 2.18.4 so that both wing segment twist angles are measured as 
rotations from a common line parallel to the airstream. This changes the equation for L, and the 
equation for the torsional reaction between the two wing segments. 


| V _ offset=h tana 


a i / \«— h —+ — 


gaps exaggerated oe 
for illustration -“@ 


section A-A 


Figure 2.18.4 - Re-definition of wing segment elastic twist angles 
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The torsional moment equilibrium equations for this configuration are developed using the free body 
diagram shown in Figure 2.18.5. 


R 


htanA 


Figure 2.18.5 - Free body diagram with new twist angle coordinates 


These equations are 


—K,6,+ Le+ K,(6,-6,)—LhtanA =0 (2.18.10) 
-K,(6,-6,)+Le =0 (2.18.11) 


Lift L, is now written as 


L, = gSC,,, (2.18.12) 


Substituting the lift expressions into Eqns. 2.18.10 and 2.18.11 results in the following two equations 
for torsional equilibrium 


h 
K,@, — qSeC,, A, + K,9, = K,0, +48C,,7 tan A| 0; = 0 (2.18.13) 


K,(0, - 6,)— qSC_,, e0, = 0 (2.18.14) 


Equations 2.18.13 and 2.18.14 are written in matrix form as 
2K, —Ky |} 4 1 anh A, 
— qSeC, e 
-Ky Ky \@ 0 1 A; 


The stiffness matrix is symmetrical. With g = 


(rh 2.18.15 
0 (2.18.15) 


2 we write Eqn. 2.18.15 as 


T 


h 
= G. —l+q,—tanA 
(2-q) oo ye -(a (2.18.16) 
- q, 


6, 
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The divergence dynamic pressure is computed by forming the determinant of the aeroelastic stiffness 
matrix in Eqn. 2.18.16, setting it to zero and solving for gp. The determinant is 


= = _/h 
0-a\-a,)+(-1+a, tan) =0 (2.18.17) 
e 
h 
With — = 10 the equations became 
e 


h 
2-39) +9p—1+Gp—10tanA = 0 (2.18.18) 
e 


or 
Gp - B-10tan Ag, +1=0 (2.18.19) 


Equation 2.18.19 is identical to Eqn. 2.18.8. We have shown, by example, that the divergence 
dynamic pressure is independent of coordinate system definition. 


2.18.2 Equation development using the Principle of Virtual Work 


Analytical mechanics provides an alternative way of developing static and dynamic equilibrium 
equations. Also referred to as Hamilton’s Principle and the Principle of Virtual Work, these 
methods use work and energy concepts, together with variational calculus to methodically develop 
the equations necessary to find divergence. Some readers will be familiar with energy method 
formulations while others will not. The purpose here is not to give a tutorial, but to outline the 
solution in a way that is at least understandable at a low level. 


A statement of the Principle of Virtual Work for a static system such as tose we have been 
considering is: 


“Tf a body is in static equilibrium under the action of prescribed external forces (in our case the two 
lift forces acting on the two sections) the virtual work done by these forces due to small compatible 
displacements (called virtual displacements) is equal to the change in the strain energy of the 
system. ” 


There are two new terms, virtual displacement and virtual work in this definition. If we have an 
infinitesimal displacement, say d@, then the inboard section actually rotates an infinitesimal amount 


so that some force must be applied to move it. A virtual displacement 69, is an infinitesimally 
small rotation, but this rotation is hypothetical or imaginary. The notation 6 ( ) is an operator used 
for variational calculus, much like the operator d ( ) used for calculus so that, for instance, 


5(6°) =20,66,. 


Virtual work is the product of a force times a virtual displacement. For instance, when the inboard 
section rotates a small amount, the upward deflection of the aerodynamic center is d, = e@,. If the 


lift were held constant during this actual displacement the work done would be W, = Ld, = Le@.. 
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During a virtual displacement the lift is held constant so that the virtual work is 
OW, = L.6d, = L,ed0, = qSC,, 0,e60, = qSeC,, 6,60, . (Notice that the lift is a function of the twist 
angle but we require that the lift remain fixed even as the twist angle changes. The virtual work 


OW, + Ld, = SqSeC, €6; = 2qSeC,, 6,08, .) 


Return to the original coordinate system shown in Figure 2.18.1. When we give the wing segments 
virtual displacements 60, and 6@,, the lift vectors move upward as indicated in Figure 2.18.6 


htanA 


displacement =(htan A)606, 
Figure 2.18.6 - Lift vector virtual displacements for wing segment sections 


The virtual work OW done by the two aerodynamic loads is: 
OW = Lcd, + L, 6d, (2.18.20) 


The virtual displacements dd, and dd, are the virtual displacements of the aerodynamic centers of 
the wing segments due to virtual torsional displacements 66, and 68, 


od, =e00, (2.18.21) 
and 
6d, = (e—htanA)66, + e680, (2.18.22) 


Equation 2.18.20 becomes 


OW = Led, + L,(e—htanA)66, + L,ed8, (2.18.23) 


or 


ow =| asec, 6,+qSeC, (6,+6,)—gSeC, (4, +0) ean) 
‘a ‘a ‘a e 


(2.18.24) 
+qSeC, (6,+0,) 60, 
The potential energy of the two springs is: 
1 2 1 2 
U=5Kr (6, ) +5 Kr(4) (2.18.25) 
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The virtual work done by the springs is equal to the variation dU, the result of the variational 
operator 0 ( 2 
OU = K, 0,06, + K, 0,00, (2.18.26) 


The virtual work of the lift forces (Eqn. 2.18.24) must equal the variation of the internal forces in the 
springs. 


OW = dU 


asec, 6,+qSeC, (6,+9,)-—qSeC, (8, +0, eran} 54 
‘a ‘a ‘a e 


or (2.18.27) 


+qSeC, (6,+0,) 60, = K,0,60, + K,0,60, 


Gather like terms in Eqns. 2.18.27 and write these two equations in matrix form. 
Kr 0 J4 [2- tan) [1—Atana) 0, 
= gSeC 2.18.28 
0 Fe Dra oF lhe. heen 


_  gSeC, 
With g =——— Ean. 2.18.28 becomes 
K, 
[12g +74 tana) [-a+a4tana) 0, 0 
Do igee eee 13 = (| (2.18.29) 
=¢ (l-q) : 


Calculating the determinant of the aeroelastic stiffness matrix in Eqn. 2.18.29 and setting it to zero 
gives the characteristic equation: 


Gp t+ Gp(10tanA —3)+1=0 (2.18.30) 


This characteristic equation for divergence is identical to Eqns. 2.18.8 and 2.18.19 developed 
previously. 


As mentioned previously, the matrix relationship in Eqn. 2.18.29 is different than Eqn. 2.18.7 even 
though the same twist coordinates are used. Note also that when g = 0, Egn. 2.18.29 is symmetric, 
while Eqn. 2.18.7 is not. A stiffness matrix found using the Principle of Virtual Work will always 
yield a symmetrical stiffness matrix. 
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2.19 Torsional divergence of an unswept flexible wing - differential equation model 
So far we have limited our static aeroelastic models to one, two or three degrees of freedom in which 
the aerodynamic surfaces and structures were 
lumped into wing segments with 
interconnected torsional _ springs. The 
interaction between aerodynamic loads and 
twisting deformation of slender wings can be 
idealized so that equilibrium conditions are 
represented with differential equations. In this 
case the structure is uniformly distributed 
along the wing so that any deformation will 
be continuously distributed along the wing. 
In this section we will use this classic model 
(discussed routinely in all aeroelasticity texts) 
to discover the effects of aeroelasticity on a higher fidelity wing model. In particular, the model 
allows us to predict how the 

airloads redistribute themselves | dx 


because of aeroelasticity. T [| T+dT/dx V 


| 
line of aerodynamic centers 


reference, elastic axis 


line of centers of mass 
Figure 2.19.1 - Unswept wing planform 


ip 


The continuous flexible wing root 
structure, shown in a planform ae 
view in Figure 2.19.1, has only 
twisting freedom, despite the fact 


that the wing will also bend boundary element boundary element 
under the distributed 
aerodynamic loads. The Figure 2.19.2. Unswept wing torsion idealization showing a 


: typical torsion element 
structure extends continuously 


from the root to the tip and is embedded within a wing with constant chord. The model has a 
reference line of aerodynamic centers and a line of shear centers (called the elastic axis) so that 
there is a constant distance, e, between the two. In addition the wing mass is accounted for by 
distributing a constant mass per unit length along the wing at a distance d behind the elastic axis. 


Figure 2.19.2 shows the continuous wing structure subdivided into infinitesimally small sections 
with lengths dx together with a typical free body diagram. Lift forces /(x) are distributed, directed 
upward from the page, along the line of aerodynamic centers. Lift forces produce a torque about the 
wing elastic axis. The lift per unit length, in terms of the initial angle of attack and the distributed 
twist, Ox), is: 

I(x)= qe c,, @ = qca,(a, + 8) (2.19.1) 


The notation a, is used as shorthand for the local lift curve slope c, . In reality the local lift curve 


slope is not constant but changes with wing location, but this model ignores this feature. This 
aerodynamic representation is often referred to as strip theory in which the aerodynamic interaction 
between elements is ignored. 


In addition, an aerodynamic pitching moment is distributed along the span at the aerodynamic 
center. 
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i =Gqe (2.19.2) 


mac 


where c,. is the sectional moment coefficient (also assumed constant along the wing) about the 


mac 


aerodynamic centerline. Finally, the distributed wing weight w creates a torque, f,,. 


t, = nwd (2.19.3) 


where vis the load factor normal to the wing surface. Combining these loads, the distributed torque, 
shown in Figure 2.19.2, and positive as shown, is: 


t(x) =qcea,(@,+0)+qc'c,,,. +nwd 


rie (2.19.4) 
(note a, =C, ) 


The relationship between the twist deformation along the length of this structure and an internal 
cross-sectional resultant torque, T(x), shown in Figure 2.19.2, is given by the relationship: 


T(x)=GJ ae (2.19.5) 
dx 


GJ is the effective torsional stiffness of the wing cross-section. 


Equilibrium of internal and external torques about the x-axis gives: 


EM, =0=T-rde—[T+ Zax) =(1-2 an (2.19.6) 
xX 
or 


dT 
—— =-+4(x) (2.19.7) 
dx 


By using Eqn. 2.19.7, we can derive a relation between the torsional deformation @ and the applied 


external torque f(x). This is 
dT d dé 
= [as = —t(x) (2.19.8) 


Substituting the expressions for the distributed torque, Eqn. 2.19.8, becomes: 


(as ao + qcea,0 = — (qcea,a, +qc’c,, 


+nwd 2.19.10 
7 = : ) ( ) 


a 


K =(qcea,a, + qc°c +nwd)/GJ. With GJ constant along the wing, Eqn. 2.19.10 is 


mac 


do 
caer (aed @=—-K . Define an aeroelastic parameter, Rh = gee ; Eqn. 2.19.10 is: 
dx GJ GJ 


6° +NO=-K (2.19.11) 
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The solution to Eqn. 2.19.11 is: 
(x) = Asin Ax+ BcosAx—K/7 (2.19.12) 


where A and B are unknown constants. At the wing root, x=0, we have a boundary condition. 


K 
(0) =0 = B~ 75 


K 
that gives B= 2 . The expression for O(x) becomes: 


K 
A(x) = Asin Ax—— (1—cos x) (2.19.13) 
At the wing tip, x=L, the torque T(L) is zero because it is a free end. From Eqn. 2.19.5, 


T (L) =GJ0'(L) =0=GJ [Arcos eee aL | 
A (2.19.14) 


Solving for A and substituting it into Eqn. 2.19.13, we have: 


O(x)= [1-¢0s Ax —(tan AL) (sin Ax) | 
(2.19.15a,b) 


a(s)=-[a.+] =] 4] te ie eavAn= (tin At) (ain 22) 


qcea, 


Equations 2.19.15(a) and (b) describe the twist distribution along the wing from root to tip. We 
should not forget that the wing also has bending 
deflection, but bending deflection does not affect 
unswept wing airloads. 


Equation 2.19.15(a) shows that the tan AL term 
has a big effect on wing twist. A plot of tan AL 


tan(AL) 


5 
vs. AL is shown in Figure 2.19.3. At AL= me 


N 
tan AL becomes infinite. If aati small 


initial wing loads produce (theoretically) infinite 
torsional deformations. Physically, this cannot 
happen because structural and/or aerodynamic AL 

nonlinear effects would make our model Figure 2.19.3- Twist amplification coefficient tan AL 
assumptions invalid. 


0 
0 T/10 20/10 30/10 40/10 m/2 
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We suspect (and will prove in the next section) that torsional divergence of this wing model occurs 


when 
2 
yh a -(2) (2.19.16) 
GJ 2L 
or 
2 2 

-() GJ -(4) [| 1 idan 
gr cea, 2 L }\ Sec, own 


K. 
The typical section idealization predicted the divergence dynamic pressure to be gp = Ks ae 4 
La 


the dynamic pressure at divergence is directly 
proportional to the stiffness parameter Gy 
and inversely proportional to wing area S$; the 
greater the wing semi-span, the lower the 
divergence gq. Notice also that the sectional lift 
curve slope c,, =a, has replaced the wing lift 


curve slope in __ the 


_K;, 
Io ~ /SeC,, ° 


Figure 2.19.4 plots Eqn. 2.19.15(b) when K=1 
(the plot is @K) and q equal to 10% of the 
divergence dynamic pressure. This twist 
distribution is very close to that found when 
aeroelastic effects are ignored. 


equation _ for 


twist angle per wing root 


tion at 10% 


angle of attack (radians) 


nondimensional distance from wing root 


Figure 2.19.4 - Twist angle (radians) vs. distance 
from root at 10% divergence q 


2.19.1 The perturbation solution for differential equation divergence dynamic 


pressure 


Stability is determined by considering perturbations away from a static equilibrium position. We 
develop the perturbation equilibrium equation by taking Eqn. 2.19.10 and setting the right hand side 


to zero. 
fe {Metal g 9 (2.19.18) 
In this case, @ represents the perturbation twist distribution along the wing. The is 
Ox) = Asin Ax + Boos Ax (2.19.19) 
with the constants A and B undetermined. The boundary conditions are 
(0) =0 (2.19.20) 
and T(L)= GJ@’(L)= 0 (2.19.21) 
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Substituting Eqn. 2.19.20 into Eqn. 2.19.19, we have 
6(0) = Asind*0+ BcosA*0=0 (2.19.22) 
or A(0)+ BU) =0 (2.19.23) 


The wing free end torque condition, Eqn. 2.19.21, becomes 


O'= AAcosAL-— BAsin AL (2.19.24) 


Equations 2.19.23 and 2.19.25 are written in matrix form as 


hee AL) oe aha a {of (2.19.26) 


This is an eigenvalue problem. The determinant of the 2 x 2 matrix in Eqn. 2.19.27 must be zero if 
A and B are to be nontrivial. This leads to the condition for neutral stability at divergence. 


A=AcosAL=0 (2.19.27) 


This equation has the following solutions 


wae “Ggecy 1 
GJ 


VL = 


(2.19.28) 


where n is an odd number so that 


_a 3a (Qnti)a 


AL ; saree (2.19.29) 
2D 2 


The lowest value in Equation 2.19.28 (n=/) defines the divergence dynamic pressure. This leads to 
the divergence condition when n = J. 


(2.19.30) 
4 GJ 


Equations. 2.19.30 and 2.19.16 are identical. 


When the eigenvalue in Eqn. 2.19.28 and the boundary conditions are inserted into the solution for 
the eigenfunction, Eqn. 2.19.19 we get the divergence mode shape. 
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O(x)=AsinAx= asin( 22) 
2L 


This is the mode shape for divergence of this 


continuous wing model. It has an arbitrary =f 
amplitude, A, and is a half sine wave. & 

2 
Figure 2.19.5 shows the effect of increasing the = 
value of the nondimensional dynamic pressure to 
75% of the divergence dynamic pressure. The 
same plot for the 10% divergence dynamic "90 02 04 06 08 41.0 
pressure was previously presented in Figure 2.19.4 nondimensional distance from wing root 
and is barely visible in Figure 2.19.5. An increase Figure 2.19.5- Twist angle parameter, @K vs. 
in dynamic pressure from 10% to 50% and 75% of distance from root for three different dynamic 


the divergence dynamic pressure creates a ls 


substantial increase in twist (and lift) if the angle of attack is not changed. 


2.19.2 Example - Divergence of a wing with flexible wing/fuselage attachment 

The wing root attachment to an aircraft fuselage is fairly rigid. However, to investigate the effect of 
a flexible wing root/fuselage connection we will use the differential equation model. The planform 
of our uniform cross-section, uncambered, unswept wing is shown in Figure 2.19.6 (the wing semi- 
span is the dimension b instead of L and the independent coordinate is now y, not x). The wing- 
fuselage attachment joint is torsionally flexible; a torsional spring with spring stiffness K, (in- 
Ib./rad.) is attached between the wing root and the fuselage so that the attachment restoring moment 
is K,@(0) when the wing root twists an amount O(0). In the previous example, the root was 


restrained so that 0 (0) =0. 


V line of 
aerodynamic 
centers 


tow 


\ elastic axis 


flexible attachment 


fuselage 
Figure 2.19.6 - Continuous wing with flexible wing/fuselage attachment 


The purpose of this example is to investigate how wing attachment flexibility changes the 
divergence speed of this wing model. To do this we will: 
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e develop the boundary conditions for wing divergence and show that the value of the 
dimensionless stiffness ratio B= OY b determines the wing divergence speed. 
T 


e develop the characteristic equation for wing divergence and solve this equation for several 
values of f. 


The equilibrium equation for the continuous wing does not change. From Eqn. 2.19.18, the 
homogenous torsional equilibrium equation for the wing perturbed twist distribution, A(y), at 


qcec, 


divergence, is 9"+ 4°79 =0 with A? = 
GJ 


Figure 2.19.7 - Free body diagram of wing/fuselage junction 


The internal torques present at the wing root are shown in Figure 2.19.7. The attachment spring 
torque, M,, and the internal wing torque, 7’, at the wing root are defined as 


M, =K,0(0) (2.19.32a) 
and 

T (0) =GJ@'(0) (2.19.32b) 
Torque equilibrium at y = 0 is a boundary condition. As indicated in Figure 2.19.7, the torque from 
the torsion spring, M,, and the torque at the wing root must be in equilibrium so that 
T(0O) — M, =0, or 

K,@(0) =GJ@'(0) (2.19.33) 


K; , 
=! 9(0) = 0’(0 2.19.34 
or a (0) = @’(0) ( ) 


At the wing tip, y = b, the torque is zero. 
Since GJ is nonzero; we divide by GJ to write the second boundary condition as: 


0’(b)=0 (2.19.36) 
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Equations. 2.19.18, 2.19.34 and 2.19.36 form an eigenvalue problem. The general solution to the 
differential equation is O(y) = Asin Ay + Bcos/y so that the derivative of the twist along the wing 
is 0’(y) = AAcosAy—BAsinAy. The boundary condition at the wing/fuselage junction (Eqn. 
2.19.34) becomes: 


Bt gj) BEE waar 


GJ GJ 
pice = AAb 
GJ 
B=ABAb 
Finally ABAb- B=0 with B= = (2.19.37) 


T 
The torque free condition at the wing tip (Eqn. 2.19.36) is written as: 
0’(b) = AA cos Ab— BAsin Ab = 0 (2.19.38) 


Written in matrix form, Eqns. 2.19.37 and 2.19.38 become: 


Bab —1 A 0) 
. = (2.19.39) 
AcosAb —AsinAb||B 0 


The determinant of this matrix equation is the characteristic equation for the divergence problem. 
A= A(-BAbsin Ab + cosAb)= 0 (2.19.40) 


When the attachment is rigid, K; is infinite, B = 0, and Eqn. 2.19.40 reduces to Eqn. 2.19.27. Since 
the airspeed is greater than zero, 2 #0, so we can divide Eqn. 2.19.40 by 4 to find the 
characteristic equation for divergence of this wing model in terms of the eigenvalue Ab . 


PAbtandAb=1 (2.19.41) 


or tan Ab = an (2.19.42) 


1 
The solution for Ab is the intersection between the two functions, tan Ab and ——. When f= 0 


BAb- 
(Kr is infinite) the root is rigid and the eigenvalue for Eqn. 2.19.42 is Ab = Mf, Figure 2.19.8 
shows the solution to Eqn. 2.19.42 for three different values of B: B= vA p= Vi. and B= Vas 


Ab 
The abscissa in Figure 2.19.8 is — where the clamped wing root reference value is + ; 


%, 
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Figure 2.19.8 - Eigenvalue solutions for three different values of B 


For example, with = 0.25, Figure 2.19.8 shows that 


Ab qcec, 
—=0.8051 or Ab= “b=1.2646 (2.19.43) 
a, GJ 


Stiffer (larger) values of the wing/fuselage attachment spring stiffness correspond to small values of 
B, so, as the fuselage attachment becomes stiffer, the eigenvalue associated with wing divergence is 


closer to that associated with a torsionally clamped wing root. 


Using Eqn. 2.19.43 and the divergence dynamic pressure for B = 0.25 is: 


ge a 0 -(2*) 
GJ b 
2 2 
. a -() Gye [o.so5% Se (2.19.44) 
b cec, cb°ec, 2) cb’ec, 


since the airspeed involves a square root of the answer in Eqn. 2.19.44 we conclude that with 
B= 0.25 the flexibility of the wing root attachment reduces the wing divergence airspeed by about 
20% when compared to a rigid attachment. 


Figure 2.19.9 shows a plot of the wing divergence speed relative to the clamped wing divergence 
speed, as a function of £. 
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Figure 2.19.9 - The effect of the support stiffness conditions on divergence 


From Eqn. 2.19.39 we find that AfAb = B so that the eigenfunction for divergence in this case is 
O(y)=AsinAy+ Boos Ay= A(sinAy+ BAbcos Ay). In this case we have both sine and cosine 
participation in the mode shape. 


2.19.3 — Stiffness criteria — an energy approach 


The importance of wing torsional stiffness was recognized in the early years of aviation, both for 
load and stress determination and flutter. Because high fidelity analysis was not available, much of 
the certification for safety was left to flight test. While this ensured a safe airplane, if a problem was 
discovered during flight testing, the modifications were likely to be costly. During the early 
development, instead of extensive analysis, designers were likely to use stiffness criteria, requiring 
that the wing torsional stiffness exceed a certain level. 


Flax'’ developed one such criterion for divergence. His analysis began with the continuous model 
shown in Figure 2.19.1, although this model was not restricted to a uniform planform. Assume that 
we are operating at the divergence dynamic pressure and we give the system the usual perturbation, 
in this case a perturbation from static equilibrium. From Conservation of Energy we know that the 
wing structure strain energy and the work done by the external aerodynamic loads must be equal. 
This gives the following equation: 


lps (doy 1p 
U =o) GJ (<2) dx=W.,.. =| (asec, (x))(e0(x)) dx (2.19.48) 


The factor 2 appears in the aerodynamic work expression because the load is assumed to have been 
very slowly applied from a zero value. In Eqn. 2.19.48 the torsional stiffness is a function of the 
spanwise coordinate, x. Similarly, the lift distribution expression need not be so simple as it was in 
our previous analysis. We recognize that the solution for x) involves considerable effort and a 
closed-form solution may not even exist. 
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The next step is to choose a reference torsional stiffness, GJ,, at the wing root for instance, and 
similar reference values for the wing chord, c,, and the elastic axis offset, e,. This produces the 
following relationship. 


2 
b| GJ \( dé b 
GJ, F (2 (4 AX = GpC,e, I, [=< (£4 ) dx (2.19.49) 


Now define the torsional deflection and lift distribution as 0 =@, f(x) and ¢,O=c, Of (x) 


then we have: 
2 
>o{ GJ \( df _ 2fo{ Cc @ 29 
GJ 0: (ZI) dx= nC, €,9; I, ot a dx (2.19.50) 


br c \l e 
Dp o£ 061, i Ere 
Co €, 
and GJ, = 3 (2.19.51) 
b| GJ \{ d 
Nea eek 
0\ GJ, }\ dy 
Equation 2.19.51 is a stiffness criterion for divergence at the specified divergence dynamic pressure. 
Although we do not know the exact deflected shape at divergence (the eigenfunction) we can 


approximate the shape f(x) and estimate the required torsional stiffness. If the wing has constant 
chord the expression reduces to 


b| C e 2 _ 
10 0801, [, [\\u )dx _ I.E 1, I, (f*) dx 
2 of df 2 
\ (4) o 


With an approximation for the torsional divergence deflected shape an approximation for the wing 
root torsional stiffness is generated. When the exact mode shape is known, Eqn. 2.19.52 provides 
the exact required value for torsional stiffness. (Use the mode shape sin Ax to find the required GJ 
and compare it to the exact value computed in Section 2.19.1.) 


(2.19.52) 


Summary 


Aeroelastic models must include enough details to allow us to predict what happens to the lifting 
surface structure when we fly at different speeds and altitudes. It does not need to have exquisite 
detail like a stress analysis model. The models considered in this chapter were simple and assumed 
that both the structure and aerodynamic behavior is linear - output is proportional to the input. We 
introduced the typical section model with pitching (structural twisting) and “bending” (called 
plunging) displacement degrees of freedom. This model illustrated that aeroelasticity is a natural 
feedback process between structural deflections and aerodynamic forces and moments. 
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From a structural mechanics standpoint, aeroelastic problems are statically indeterminate problems; 
the structural deformation creates loads that enter into the equilibrium equations. This chapter has 
presented an overview of the matrix method structural analysis approach and introduced the shear 
center and elastic axis concepts as structural reference points. It also illustrated the relationship 
between the strain energy stored in a structure and the structural stiffness matrix. We reviewed 
fundamental low speed aerodynamics and introduced aerodynamic coefficients and aerodynamic 
reference points — the center of pressure and the aerodynamic center. 


Aeroelastic and load amplification and static stability is a stiffness problem; aerodynamic loads 
create aerodynamic stiffness (or “de-stiffness”) just as structural geometry creates elastic stiffness. 
The effect of aerodynamic stiffness is to reduce unswept wing torsional stiffness; this effect leads to 
aerodynamic load amplification and ultimately static instability known as divergence. 


If the total wing lift is held constant (the wing angle of attack can change) twist deformation is small, 
even at the divergence speed. It is only when the angle of attack is fixed and the load is allowed to 
become large that twist deformation becomes unbounded at divergence and extremely large even at 
speeds considerably below the divergence speed. As a result, we place the wing in static 
equilibrium with a small (or even zero) angle of attack, perturb the wing and then observe what 
happens to the response. 


With a linear analysis, at divergence we encounter a “bifurcation” situation in which there are 
multiple static equilibrium states in addition to the original equilibrium state. This is Euler’s criteria 
for static stability; the boundary between stable and unstable states occurs at a transition point 
defined by the appearance of multiple static equilibrium states. We showed that a nonlinear model 
also predicts the appearance of multiple equilibrium states, but with finite, not infinite, deformations. 


We also examined the effects of flow compressibility on divergence; flow compressibility introduces 
a compressibility correction factor, the Prandtl-Glauert transformation. As a result, the Mach 
number is introduced into the problem. We found a “match point” for which the Mach number at a 
given altitude matches the Mach number at divergence. 


Divergence is a theoretical problem that requires that the wing support is securely fixed to an 
immovable aircraft. This “immovable aircraft” condition cannot occur in flight. Inertia effects such 
as weight distribution cannot appear in the results of a static stability analysis since they cannot 
produce deflection dependent loads without acceleration. Still, the theoretical wing divergence 
dynamic pressure allows us to understand when aeroelasticity is likely to be an important problem. 


We also considered multi-degree of freedom problems and determined that the divergence problem 
is an eigenvalue problem involving the search for self-equilibrating conditions, combinations of 
dynamic pressures (eigenvalues) and mode shapes (eigenvectors or eigenfunctions). The multi- 
degree of freedom eigenvalue problem can be solved in simple cases by forming the determinant of 
the aeroelastic stiffness matrix and setting it to zero. This is the so-called characteristic equation, a 
polynomial involving the divergence dynamic pressure. 
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We also introduced the Virtual Work approach to developing the required homogeneous equilibrium 
equations for divergence analysis. In addition we developed a differential equation model, as well as 
presenting a torsional stiffness criteria, based on an energy approach, to define the minimum 
required torsional stiffness for specified divergence dynamic pressure. 


In the next sections we will continue aeroelastic studies by examining the effects of static 
aeroelasticity on control effectiveness, particularly roll effectiveness, and then examine how we 
might link together control systems and the wing to change aeroelastic stability. 


2.20 Control surface effectiveness 


A wing with planform area S is shown in cross-sectional view in Figure 2.20.1. A control surface 
(an aileron in the case of a wing, a rudder in the case of a vertical tail or an elevator in the case of a 
horizontal tail) is connected to the wing section by a rigid hinge. The wing may be cambered 
(although the cross-section in Figure 2.20.1 is 


not); control surface rotation 6, (positive 


oO 


Lift 
Mac ee spring K; 
| shear center 


clockwise) creates additional camber, lift and a 
pitching moment that twist the wing an amount 0. 
As we will show in this section, the twisting Oy+ 8 
deflection creates negative lift, with the result that —\V pest 
the amount of lift generated by a flexible wing is ! 
less than that for a rigid wing. 


I 
Figure 2.20.1 - Torsionally flexible wing segment 


Control effectiveness is defined as the ability of a with moveable trailing edge control surface. 


control surface such as an elevator, aileron or 

rudder to produce lift or a moment to change the airplane pitching moment, rolling moment, or 
yawing moment, respectively. Horizontal tail, wing or vertical stabilizer torsional flexibility reduces 
control effectiveness of unswept lifting surfaces. 


When the control surface is deflected, the flexible wing lift is 

L=qSC, (a, + 0) + qSC,,6, (2.20.1) 
If there were no twisting, the lift due to the control deflection is L, = qSC, . é, . The twisting 
moment about the wing aerodynamic center is: 

M gc = WSCOyyc + gScCyac, é, (2.20.2) 


The terms Cy, andC,,,-, are defined as 


aC 
Cis = 35 


OCuac, 
0d 


(2.30.3) 


and Cyac, = 


Since positive 6, causes a negative C,,. value, this term results in nose-down wing twist that 
Oo 


decreases the wing angle of attack. Summing moments (clockwise positive) about the shear center, 


Le+M,-=M,=K,0 (2.20.4) 
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Substituting Eqns. 2.20.1 and 2.20.2 into 2.20.4 and solving for 0, 


gSeC, a.+(2)Sefrax{e +(£] Cue Je 
e) €, ae 3 a 


K, —qSeC, 


(2.20.5) 


The effects of ,, and 6, on @ are additive. We confine our attention solely to the second term of 
Eqn. 2.20.5, written as 


1+ [<) | ae 
Sec e 
a eee ‘8 (2.20.6) 
é, K, _ Ee a 
Ip 
From Eqns. 2.20.1 and 2.20.5 the lift due to O, is 
L=as[ 6, +6 Ja 
or 
Sc{ C 
nn (Bs |(-ne 
T Ls 
L=qsC, 6, ° (2.20.7) 
‘Oo 1-4 
Ip 


The second term in the numerator in Eqn. 2.20.7 contains a reference dynamic pressure, called the 
reversal dynamic pressure, qr: 


EK, |G, 
r= E 2 (2.20.8) 
SC aoe (Gs. 


The term —Grac, is a positive number so that the reversal dynamic pressure is positive. The ratio 
CL, OC, 0a oa 

= x — 
C, 05 JC, o& 


(2.20.9) 


a ratio that represents the change in wing angle of attack due to a unit control deflection, 5,. 
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Two-dimensional airfoil control surface aerodynamic coefficients 
In the late 1920's Herman Glauert developed a theory and the equations to predict coefficients C/, 


and Cus for a two dimensional airfoil (as opposed to the wing model we have here) as functions of 


& 
the ratio of a “‘flap-to-chord dimension,” defined as EF = — » The flap is what we are calling the 
c 


control surface and its chord dimension is cs. These two relationships are 


aoe —(cos” (1- 2£)+2JE(1-£)) (2.20.10) 


Cie 


ents 2 (ir) (02 BE (2.20.11) 


6; a 
a 


Note that, for the two-dimensional airfoil in incompressible flow, the lift curve slope is c, =27. 


Before continuing with the aeroelastic analysis, let’s examine how the flap-to-chord ratio affects the 
location of the airfoil center of pressure location due to the control deflection on a two-dimensional 
airfoil. 


Effect of control surface size on center of pressure location 


The control surface deflection creates lift and a nose-down pitching moment. The additional lift 
increases twist while the nose-down pitching moment decreases twist. Which effect dominates? To 
answer this question, let’s return to the three-dimensional wing model. 


Figures 2.20.2 (a) and (b) show two different ways of representing the lift and pitching moment due 
to 6,. Figure 2.20.2(a) shows the lift acting at the aerodynamic center, together with the pitching 
moment due to the aileron deflection. Figure 2.20.2(b) shows the lift due to the same aileron 
deflection acting at the center of pressure location so that there is no moment. The center of pressure 
is located a distance d downstream (aft) of the airfoil quarter chord. 


aileron lift = SC | 555 


aileron lift = qSC ; 55 
MAC = QScCync 59 


aerodynamic center 


/ at es 
mp gy, do So 
: / 
(a) aerodynamic center (b) flap aerodynamic center 
coordinates coordinates 


Figure 2.20.2 (a,b) - Force and moment equivalence for aileron deflection showing center of pressure 
coordinate 


Copyright Terrence A. Weisshaar 
Chapter 2 


page 120 - Static aeroelasticity-structural loads and performance 


The two representations must produce the same moment about the % chord point. This is called 
moment equivalence and is written as 


—Ld = qSdC,,.6, = M ac, = qScCyac, 9 
—C,,.d =CCyac, 


d Cuac 5 
: e (2.20.12) 
O 
For the two dimensional airfoil, this gives 
d Ce 1-E).{1-E)E 
=-—= ( ) ( ) (2.20.13) 
Cc 


Gs (cos (1-2) +2,/E(1-E)} 


The airfoil ratio ee in Eqn. 2.20.13 is plotted in Figure 2.20.3 as a function of aileron flap-to-chord 


ratio, E. This figure shows that when the 
aileron chord is extremely small, the 


center of pressure is near the midchord 0.25 


(d/e = 0.25) but then moves forward 5 
towards the 4 chord as E increases. 8 -_ 0.20 
o6 
This plot explains why wings with some go 
control surfaces will reverse. The dashed 5 = 0.10 
red line represents a possible location of ss 
the shear center, in this case it is 35% of 8% 
the chord aft of the leading edge. The 2 
roughly triangular area above the dashed = 
line represents a design region where the 
aileron center of pressure is behind (aft) aileron flap / airfoil chord 
of the shear center and configurations that (E=cr/c) 
will twist the wing nose down, even as the 
control increases lift according to the Figure 2.20.3 - Aileron center of pressure location with 


respect to the airfoil 4 chord vs. to chord ratio 
relationship L, = qSC, . é, . As a result, . fe flap 


the sum of the fundamental lift 

L. = qSC i é, and the lift due to twist is smaller than the fundamental lift. Only when the control 
surface has a large chord will the center of pressure lie forward of the elastic axis. In this case the 
flap to chord ratio would have to be about 0.5, a very large control surface. 


The reversal dynamic pressure and its relationship to divergence 


T 


SeC, 


‘a 


We already know the expression for divergence dynamic pressure for this wing is qd, = 


Where is the reversal dynamic pressure in comparison to divergence? Substituting Eqn. 2.20.13 into 
Eqn. 2.20.8 the reversal dynamic pressure expression becomes 
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= C = c 
Gp = K, | |. ee ae = K, [£)=a0(<) (2.20.14) 


ScCyac, \ Ci, 


Equation 2.20.7 can be written as: 


4 
L=L,| —% (2.20.15) 

2 

Ip 


where L,=qS Crs 6, is the lift generated by control surface deflection 5, ignoring aeroelastic 


effects. Control effectiveness, when the primary purpose of the control deflection is to produce 
lift, is defined as: 


ie: 

aS Ir 
be |y-4 (2.20.16) 

Ip 


Assuming that the wing dynamic pressure is less than the divergence dynamic pressure, then, as 
dynamic pressure increases and approaches qe, the ratio L/ L, decreases and becomes zero when 
qd =p. At this point no net lift is produced by the control deflection, although the wing segment 
will twist and thus cancel the effects of control deflection lift (see Eqn. 2.20.6). 


Two different plots are necessary to examine lift effectiveness since qr can be either greater than or 


less than 4p, depending on where the aileron center of pressure is (which in turn depends on how 
large the aileron is). The divergence dynamic pressure appears in Eqn. 2.20.16 so we define a 
parameter, R, as 


rR —K, C., 
or WseCur, C (2.20.17) 
RaL Curc, _¢ d 7 d 
e C., ec e 


(= 2: (7. 

L -(2] Ir | _ Ir 

L. \ej/|d_q ‘) q 

i pie a (ee eect | Meee 2.20.1 
€ AR (‘ Apr ed 
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In Figure 2.20.4, the ratio L/ L, is plotted 
when R>1. Note that there are two 
limiting conditions in Figure 2.20.4. 
When gr=qp, then R=1 and rigid 
surface and flexible control effectiveness 
are identical. When R=1 then 
Y=, =-Cy,1C,,. An additional 


limiting case in Figure 2.20.4 occurs when 
R is very large (q p >> dp) So that 


flexible lift to rigid lift ratio 


200 02 04 06 08 1.0 lim Lg 


— =]|-— (2.20.19) 
/Oreversal (Adiv > Ireversal) Rok, Ir 
Figure 2.20.4 - Control effectiveness as a function of 
dynamic pressure ratio, q / dp (R>1). Notice that Eqn. 2.20.19 plots as a straight 
line in Figure 2.20.4. 


Figure 2.20.5 plots control effectiveness 
vs. dynamic pressure when R<1. In 
this case, control effectiveness always 
increases when dynamic pressure 
increases. For instance, when R=0.5 
the aileron control effectiveness is 
always greater than | in the range of 


practical values of ¢/qp. 


Flexible to rigid lift ratio 


The reversal dynamic pressure defined 
in Eqn. 2.20.8 contains Mach number 
dependent aerodynamic derivatives. 
We can rewrite Eqn. 2.20.8 as 


0.0 0.5 1.0 1.5 2.0 
W/Greversal (Adiv < Ireversal) 
Figure 2.20.5 - Control (lift) effectiveness as a function of 


q/ dp (RSD. 
aK, Gr 
gS | |e (2.20.20a) 
SeCis (Ci. 
or dr = (pn, N1-M 3 (2.20.20b) 


where qz, is the reversal dynamic pressure computed assuming incompressible flow. As in the case 
of divergence, we must match the reversal dynamic pressure including Mach number with the 


1 
atmospheric flight value Gp = 5PM = Ir |l-M,. . 
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2.20.1 Aileron effectiveness including rolling motion 


In the last section the examination of control effectiveness was restricted to a wing/control surface 
combination fixed to a wind tunnel. In this case aileron effectiveness was measured as an ability to 
create lift which was in turn transmitted to a wind tunnel support so that the wing/control surface 
could not move. This situation approximates what happens with a horizontal tail/elevator 
combination or the rudder/vertical tail 

where motion is small or nonexistent. It v=pb/2 

does not simulate wing/aileron rolling 
motion about an aircraft roll axis. The 
purpose of this section is to introduce 
another measure of control surface 
effectiveness - aileron effectiveness — in 
which motion occurs. 


aileron 


The ability to create motion, in particular 
a large terminal or steady-state roll rate 
(defined as p in Figure 2.20.6), is the Figure 2.20.6 - Rolling airplane with ailerons deflected, 
measure of wing/aileron effectiveness. es yer roms ie Tear. 

Figure 2.20.6 shows a rear view of an 

airplane flying at an airspeed V-. 

Downward aileron deflection on the right wing (and upward deflection on the left accelerates the 
wing upward to produce a roll rate p(t) (radians/sec.). Airplane roll produces a velocity component 


v(y) along the wing, given by the relationship V (y ) = Py , as shown in Figure 2.20.6. If the aileron 


span is small the aileron rolling velocity v(y) at the wing/aileron surface is approximately 
A 


From the wing’s perspective, air appears to be moving downward with velocity v and rearward with 


airspeed V so that there is an apparent negative angle of attack, @, =— y, . This creates a lift term 


WY pb 
OG, (=) =«86, (2) called the aerodynamic damping force that acts downward, opposing 


upward motion. This force produces a rolling moment that retards motion and is called “damping in 
D: stim 
roll.” Eventually, as the roll rate increases and v= increases, the damping in roll becomes so 


large that the roll rate will cease to increase and p and v will become constant. The unsteady 
aerodynamic effects during this event are complicated and will not be discussed here. 


In figure 2.20.6 the lift on the left wing/aileron combination at the wing aerodynamic center has 
three input components: 


. 
L=4qSC,, [e-Z}+asc,,6 (2.20.21a) 


This lift creates a rolling moment M,,;,(t) that accelerates the wing upward so that (the factor of 4 on 
the right hand side is there because there is an aileron on the left wing doing the same thing). 
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b b b Td 
M,(t) -1(2) =gSC,, [Fl(e-2} +050, Ga oe (2.20.21b) 


Since 8 is unknown, the roll rate is an unknown. The steady-state condition in which yy i 0 is 


seldom reached in flight, but it represents a limiting condition. To find a reversal speed we first find 
the (unknown) steady state rolling velocity, p, as a function of the aileron deflection angle 6, by 


using a typical section idealization to simulate wing/aileron motion at a constant speed. 


To capture this situation with an idealization we imagine a very tall wind tunnel with a wing on 
frictionless bearings shown in Figure 2.20.7 that permit upward or downward frictionless movement. 


The velocity generated by the aileron is v= = : 


Figure 2.20.7 - Rolling effectiveness idealization showing wing/aileron combination free to move 


Because of the opposing aerodynamic damping force, the wing will accelerate until v becomes so 
large that the acceleration is zero. In this case Eqn. 2.20.21b becomes 


M, (t)= (3) = SC, [S)[e-2}+ase, (4)e- 0 


or (2.20.22) 


L=4qSC,, [e-2 gSC 5 =0 


Equation 2.20.22 has two unknown, the terminal value of v and the wing twist @ he problem 
resembles a statics problem. We use Eqn. 2.20.22 to solve for the terminal velocity. 


C, 
=04+—£6, (2.20.23) 


This analysis can be confusing because in the previous section we used the L=0 condition to define 
control effectiveness. Now we use this condition to define the terminal velocity. Next we solve for 
the wing twist. 
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Figure 2.20.8 — Free body diagram — wing/aileron combination free to move upward with constant speed, v. 


To find the twist angle we need to consider the wing/aileron free-body diagram shown in Figure 
2.20.8. A torsion spring with a spring constant Ky lb-in/radian restrains the wing twist, @. 
Although the lift is zero when the terminal condition is reached, the wing pitching moment, M,.#0 
and is resisted by the torque in the torsion spring. Moment equilibrium (with L=0) requires that: 


K,0=M je = gSCCyac, ©, (2.20.24a) 


Cy 
or 6=qSc—6, (2.20.24b) 


T 


1) 
From Eqn. 2.20.23, the steady state or terminal velocity, v= = , 1S: 


6 Cc, \| 
YD. is wt) ca) (2.20.25) 


The pitching moment coefficient is negative so, if the dynamic pressure is large enough, a 
downward aileron rotation will produce a negative terminal velocity. The reversal condition occurs 
when Eqn. 2.20.25 is zero. 


7-2 = inSe{ Se) = 6, 
as (2.20.26) 
dn =-| 2 “ts 

SiO ype EC, 


The reversal dynamic pressure in Eqn. 2.20.26 is identical to that found in the previous section. We 
now write v/V as 
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ee ae oem | een 
Vo OV C ds 
ae (2.20.27a,b) 
C 
p = 7 Vv ey J, 
b Cy AR 


b 
The term P ve in Eqn. 2.20.27a is a helix angle traced out in space by the wing/aileron located at 


distance b/2 from the roll axis, rolling at a constant rate p and moving forward at airspeed V . 


Figure 2.20.9 plots Eqn. 2.20.27b (divided by b/2 and oo) for this wing/aileron model when 
Ip =150 psf. As airspeed increases, steady-state roll rate reaches a maximum and then declines 


as airspeed increases, finally going to zero at the reversal airspeed. Compare Figure 2.20.9 with 
Figure 1.17 in Chapter 1. 


Aileron reversal and aileron effectiveness were important topics for researchers in many countries 
during the 1920’s and 1930’s. Since torsional stiffness plays an important part in this phenomenon, 
the wing skin thickness becomes important since it determines a major part of torsional stiffness. 
The Japanese Zero, the Mitsubishi 
A6M2, was one of the finest 


150 150 
fighters produced before World tier A 7 
; , a. 
War I. However, the designer’s = 100 4-Oub-o plan 100 > 
objective of producing an ultra- 8 ty = 
ii ; : 2 Is, D 
ight-weight maneuverable fighter ® o 50 50 % 
led to a thin wing skin. At 160 © s g RE a 
mph (260 km/h), the Zero had a = ro aye > ——_| 0 = 
roll rate of 56° per second. % | | | LL Sy 
Because of wing _ torsional FE -50 j oD 
flexibility, roll effectiveness 
dropped to zero at about 483 km/h -100 -100 
(300 mph). This fault led U.S. 0 50 100 150 200 250 300 350 400 


airspeed at sea level 


pilots to guard against low speed 


encounters and try to get the Figure 2.20.9 - Steady state terminal damping in roll airspeed vs. 
Japanese pilots to engage at higher forward airspeed, V 
speeds where they were less 


maneuverable. 
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2.21 Two degree of freedom wing aileron reversal 


In the previous section we considered a one degree of freedom wing model and defined aileron 
effectiveness two different ways: 1) 
the ability of the aileron to increase Vv 
lift on the wing, defined as 


l 
line of aerodynamic centers 
Latex 7 a i a SS 
L.. when @,=0 and <y- 
rigid elastic axis planform area S 
or 12 
= 0,; 2) the ability of the aileron to panel 1 eae 
planform area S 
provide a rolling moment. The 
reversal airspeed was the same in > 
both cases. In this section we will 
show that the aileron reversal 
airspeed depends on the definition of gaps shown are negligible aileron surface 
aileron effectiveness if the wing 
model has more than one degree of Figure 2.21.1 - Wing/aileron idealization 


freedom. 


The configuration shown in Figure 2.21.1, 
previously considered in Section 2.14, is a 

wing idealization with two elastically 
connected surfaces, each with planform area i 
S, with an aileron attached to the outer 

panel. For the present example the 

torsional spring constants have been 

changed so that they are both equal to Kr. 


As indicated in Figure 2.14.1, shown 
previously, the section twist angles are 
measured from the initial angle of attack 
position. The difference between this model 


and the model in Section 2.14 is that an moment due to 
aileron has been added to the outboard \ 
anel. 
p 0, 
Figure 2.21.2 - Wing/aileron combination free body 


Figure 2.21.2, shows a free body diagram; 


the internal torques due to these torsional diagram 


restraints are given as T; = K,;@ and T, = Kr (0, = 6,) The aerodynamic forces acting on each 
panel are, with a =0, 
L = gSC,, A (2.21.1) 
L, =qSC,, @ + gSC_, 6, (2.21.2) 


Summing forces on each panel gives two equations. 
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Section 2 (outer) 


LF, =0 

=O (2.21.3) 

Ry = Ly = qSC,, & + gSC,,6, (2.21.4) 
and (Section 1, inner) 

YF, =0 

[+ Rk, — R, =0 (2.21.5) 

R=L,+L, 


Equations 2.21.4 and 2.21.5 are not useful, except to compute the internal shear forces. On the other 
hand, summing moments about the elastic axis (positive nose-up) in Figure 2.21.1 gives two 
additional, useful, equations 


=M,=0=Le+K,(0,-6,)—K,6, (2.21.6) 
and 
iM, =0=Le-K,(8,-8,)+qScCy,56, (2.21.7) 


Substituting the expressions for L, and L, into Eqns. 2.21.6 and 2.21.7 and writing these two 
equations in matrix form gives 


[2K -qSeC, -K, 4) e fe a 
- = gSeC, | —*|- 2 2.21.8 
ik” ree K; -qseC,, |, | PC, S) pom ee 


oO 


Let's use the following data to illustrate some features of aileron reversal for this type of model. 


Cc 
Geog “Once S2ig 222! 254 (2.21.9) 
i e C,, 06 
_ gSeC, 
With the nondimensional aeroelastic parameter g = = Eqn. 2.21.8 becomes 
T 
[2-¢g -1 |fA 0 
=g 0. 2.21.10 
| -1 1-a/laf~*)-05/ % ae 


The wing segment twist angles are 
6 g\1-g 1 0 
i|_all-9 | 5, (2.21.11) 
af alt 2=9||-05 
where the determinant of the aeroelastic stiffness matrix is 
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A=qG -3g+1 C2L1) 
The characteristic equation for gp is 

Gp ~ 395 +1=0 (2.21.13) 
From this polynomial we find two eigenvalues 

dp = 9.3820 and dp = 2.618 


The smaller of these two eigenvalues is the nondimensional divergence dynamic pressure. It is not a 
function of aileron parameters. 


The total lift created by aileron input }, , (with a =0) is 
Lae =, + Ly = GSC, (O,+9,)+ GSC, 6, (2.21.14) 


From Eqn. 2.21.11, the section twist angles are 


U5 < 
0, =—— 70, (2.21.15) 
A 
and 
5q(2-G -4+0.59 
6, a et é, (2.21.16) 


Substituting these expressions for @ and @, into Eqn. 2.21.14, we find the lift on the flexible 


wing. 
[-0.5 7-0.59° 
ets gq —0.5q 
Lpex = GSC, rae ee ara 6, [98% (2.21.17) 
SCAM acs, Ee pia, 1s 
ia aa | 094-940-954 a A\8, (2.21.18) 
Le 
Simplify Eqn. 2.21.18 to obtain 
SC Bt ip oe age HOE Sst cs 
L flex -P ls ~0.59-4 +0.5q¢ += (q° - 3941) |6, (2.21.19) 
A Ce 
gO? Ns. tho ot Ce | 
or Lex = *10.9¢q —2.7¢+— |6, (2.21.20) 
ral ka 
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The total lift due to aileron displacement when the wing is rigid is: 


Lyigid = qSCy, 8, (221,21) 


ep : 
Since —* =0.4 , the ratio —“& is written as 


La rigid 
, qSCy 
oO 
—ft = A — {09° -2.79 +0.4} (2.21.22) 
Liseta qSC_ 5, 
This simplifies to 
L 9q° — 279 +4 D 
Stes 20 TAT 02.21.23) 2 
Lag Mg” —3¢ +1) - 
2 
L 3 
A plot of wa 4 (Versus q_ is . 1 
rigi oa © 
presented in Figure 2.21.3, together with = o 
the plots of the numerator and - 
denominator of Eqn. 2.21.23. ® 
© 


If the purpose of the aileron is to generate 
lift, at reversal the lift is zero, or 


nondimensional dynamic pressure 


La = 0 (2.21.24) _ 
Figure 2.21.3 - fer ys. 7 showing reversal and 
At reversal, the numerator of Eqn. rgd - 
2.21.23 is zero. This condition is written divergence 1 values. 
as 
9q, —274, +4=0 (2.21.95) 
This equation is quadratic in g, so there are two values of qp . 
Ir, = 0.1563 (2.21.26) 
Ir, = 2.8437 (2.21.27) 


Reversal occurs at the lower of these two values. 
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Divergence occurs when the denominator in Eqn. 2.21.23 becomes zero. Note that aileron 


L 
effectiveness, defined as when “** becomes zero, occurs before the divergence gq is reached. 
rigid 


The traditional measure of aileron effectiveness - generating airplane rolling moment 
The function of an aircraft aileron is to create a rolling moment, not just lift. The rolling moment for 
our two-segment idealization is computed by assuming that the lift on each wing panel is centered at 
the respective panel mid-span. This results in the expression for flexible wing rolling moment 
created by the aileron input. 


b 3b 
M hex = (L—)+(, =) (2.21.28) 
4 4 
where the wing semi-span dimension is b. We begin with Eqns. 2.21.1 and 2.21.2 and substitute the 


previously derived expressions for the twist of each panel to get the panel lift due to aileron rotation. 
Equation 2.21.28 is then written as 


b (1 3 3 
Mptex = ? QC, Ne Gt @,) oF aS.) (2.21.29) 


2 


Cc 
Substituting expressions for @ and @, from Eqns. 2.21.15 and 2.21.16 (with “fe =0.4) 
La 


we have 


SC, (b 

M pex = laa (2) (2.7¢° — 7.19 +1.2)6, (2.21.30) 
A \4 

or 

_ gSC,,b1 2.79? — 7.19 +1.2 | 


Ma, = re) 221,31 
a aed ere ae ae 
The rolling moment created by the aileron on the rigid wing is 
3b 
Migia =F C1, % (2.21.32) 
M 
e ratio M vie is 

Mex _2.5(2.7q —7.1q +1.2 
J flex _ 2.52.7q —T1q +12) (2.21.33) 


Myigia 3(q° - 39+) 


In this case, control reversal is defined as M fice = 0; we find the reversal point by setting the 


numerator of Eqn. 2.21.33 to zero. 


2.7a¢ —7.1Gp +1.2=0 (2.21.34) 
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Equation 2.21.34 yields two values: 


ar =0.1815 (2.21.35a) 
and 


Gr = 2.448 (2.21.35b) 
Aileron reversal occurs at the lower of these two values. 
ar =0.1815 (2.21.36) 


Because the definition of effectiveness is different in this case, this reversal dynamic pressure is 
different than that given in Eqn. 2.21.27. Although we have not included the rolling motion in our 
analysis, this result is the same as we would find from that analysis. Aileron rolling moment 
effectiveness, defined as the ability to create a rolling moment (Eqn. 2.21.33), is plotted in Figure 
2.21.4. 


Nf ot 
LS 


Outboard aileron 
flex/rigid rolling moment 


‘ Z A M flex 
Figure 2.21.4 - Rolling moment effectiveness vs. 7 
rigid 


The effect of moving the outboard aileron to the inboard panel 


When we move the aileron inboard the free body diagram for the configuration becomes as shown in 
Figure 2.21.5. The lift on each panel is 


L, = qSG,,4,+9SC,,6, (221,37) 
L,= gSC,, 0, (2.21.38) 
Summing the twisting moments on each panel, we have 


<M, = 0= Le+ K;(@ — 6)— KO, + gScCy, 0, (2.21.39) 
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and mM, =0= Le- K;(0, -9,) (2.21.40) 


Ty 
7, =K,(q-a) Ts 


» 


iy 


Figure 2.21.5 - Free-body diagram for inboard aileron/panel 


These equations in matrix form are 


2K —gSeC -K 6 [c Cc... 
[2K a i i '}s qSeC, |— (<) pens | (2.21.41) 
L “a om qseC,,, | 8, oe ev C,, | 0 
With the parameter values defined previously, we get: 


[2-7 -11]/6) _{-05 
| 1 ig af=4] 0 5, (2.21.42) 


The left hand side of Eqn. 2.21.42 is identical to Eqn. 2.21.10, but the aileron input has been moved 
from the lower element to the upper element in the vector on the right hand side of Eqn. 2.21.42. 
The divergence dynamic pressure is unchanged. 


To find the lift reversal dynamic pressure for this configuration, the lift on the flexible wing is 
computed. The panel twist angles {9} are: 


@\ gfl-F 1 ]{-05 
lea 1 2-g]| 0 0, (2.21.43) 


so that 
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SGd-q) .  -0.59+0.57” 
jy aE) oe 


= . (2.21.44) 
0, = 0.54 é, 
A 
The lift created by deflecting the inboard control surface is: 
Lex = GSC, (A, + &) + GSC, 6, (2.21.45) 
Substituting the expressions in Eqn. 2.21.44 into Eqn. 2.21.45, we have 
[ 0.59? -0.59 - 0.59 | 
Lpex = SC_, ee +qSC;, 6, (2.21.46) 
Re-writing Eqn. 2.21.46 we have: 
gSC_ ~2 = 
Lnex =—Y [0.99° - 2.29 + 0.46, (2.21.47) 


With aileron reversal defined as the inability to create lift (L tex = 0) we have the following 


equation for reversal. 


0.99. —2.2G, +0.4=0 (2.21.48) 
The two solutions for @p are: 

q k, = 0.1978 (2.21.49) 

dp, = 2.2466 (2.21.50) 


Reversal occurs at the lowest value of gp. 
dp =0.1978 (2.21.51) 


Comparing the value of gp in Eqn. 2.21.51 to that in Eqn. 2.21.27 we see that the value of the 
reversal parameter g in Eqn. 2.21.51 is 1.27 times that presented in Eqn. 2.21.27. The wing with 


the inboard aileron has a higher reversal speed. However, it is not necessarily more effective in 
generating lift at some airspeeds. 


Rolling moment with the inboard aileron 


To compute the rolling moment due to the inboard aileron, we begin with the expression for rolling 
moment. 
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sat 4 ( 2) ( a 
Mae = |Get} t+ BCL, 815 (2:21.52) 
or 
Sa kal Lo 2] 


ata | 2 |g? 29 42 2.21.53 
HOS Ne AL oe ge 


Since rolling moment generation is the objective, then, from Eqn. 2.21.53, aileron reversal requires 
that the numerator in that expression be zero. This means that 


1 2 
5 aR ~2Gn +5 =0 (2.21.54) 


The solution to Eqn. 2.21.54 yields two values. 


IR, =().2112 (2.21.55a) 

IR, = 3.789 (2.21.55b) 
Reversal occurs when 

Gr =9.2112 (2.21.56) 


Comparing Eqn. 2.21.56 with Eqn. 2.21.36, we see that the roll moment reversal dynamic pressure 
with the inboard aileron is 1.163 times that of the outboard aileron so moving the aileron inboard 
increases the reversal speed. On the other hand, before it reverses, the rolling moment due to the 
outboard aileron is larger than that for the inboard aileron because it has a larger moment arm. 


With real operational aircraft, different surfaces, located at different positions on the wing, are used 
to generate roll moments at high airspeeds that at low speeds. The low speed ailerons are located 
farther outboard than the low speed ailerons, but are "locked out" after a high speed is reached so 
they are prevented from operating. To understand why this is so, let's look at the relative 
effectiveness of each of the two surfaces in this problem. 


First, compute the ratio of My, for the outboard aileron to M¥,,, for the inboard aileron, to get the 


following expression 


2 =, 
Mon = M hex joutboard = 2.79 = ge 1g +12 (2.21.57) 
ne M secrinboard 0. Sq = 2q +0.4 


or 
5.497 -14.29 +2.4 


M.. = 
ratio a —49+0.8 


(2.21.58) 


If, at a chosen dynamic pressure, M is greater than unity, then it is better to use the outboard 


ratio 


aileron for roll control. Figure 2.21.6 plots the ratio in Eqn. 2.21.58. 
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3.0 Za 

2.0 fies’ ; 
es 

Ant 

wa. tt 


0 0.05 0.10 0.15 0.20 
Dynamic pressure parameter 


Rolling moment ratio: 
outboard/inboard 


Figure 2.21.6 - Relative effectiveness of ailerons 
to create a rolling moment (Eqn. 2.21.58) showing relative effectiveness of the two ailerons 


The "cross-over" dynamic pressure at which the inboard aileron becomes more effective occurs 


when Mj; 18 equal to 1. The cross-over point is found by setting the numerator in Eqn. 2.21.58 


equal the denominator of Eqn. 2.21.58 so that, at the cross-over point 


5.49° -14.29+2.4=g7 -49+0.8 (2.21.59) 


or 


4.4g7 -10.29 +1.6=0 (2.21.60) 
The solutions to Eqn. 2.21.60 are 


g, — 0.1692 (2.21.61a) 
and 
Gy = 2.149 (2.21.61b) 


From Eqn. 2.21.61a and Figure 2.21.6 we see that the outboard aileron is more effective for moment 
generation as long as q < 0.1692. The outboard aileron configuration reverses at g, =0.1815, as 


indicated in Figure 2.21.6. At larger values of dynamic pressure the inboard surface is more 
effective but will reverse at dynamic pressures above g =0.2112. 


2.22 Artificial stabilization - wing divergence feedback control 


Since aerodynamic forces and moments are naturally fed back to create displacement dependent 
airloads that may create instabilities, we can use an “unnatural” feedback process to control loads 
and create an artificial situation to change the divergence speed. 


2.22.1 — A two dimensional example 


The purpose of this section is to show how the natural aeroelastic feedback process is changed by an 
actively controlled aileron. For our example we will use the configuration previously considered in 
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Section 2.14, Figure 2.14.1. The addition of an aileron to the outboard panel, labeled as panel 2, is 
indicated in Figure 2.22.1. Both angle of attack and control deflection are inputs. 


aero 
centers 


shear 
centers 


Figure 2.22.1 - Two degree of freedom wing idealization with aileron. 
When the aileron is deflected an amount 6, the lift on panel 2 is: 


L, = qSC,, (@, + )+ GSC, 6, (2.22.1) 


Control surface deflection also causes a pitching moment about the panel 2 axis of rotation, written 
as: 


M,=L,e+ gScCy,, }, (2.22.2) 
Expanding Eqn. 2.22.2, we have: 
C,, c Cy, 
M, =qSeC, (@,+6,)+qSeC, +| = é, (2.22.3) 
i "\C, Ves, 


We include this modified equation in our original model in Section 2.14 to find the matrix equation 
for static equilibrium. 


5-2], [-1 O7/f4l_—) [Naf Cue, eC \f0] 
go) lo allel toi c, ec, Sls (2.22.4) 


Equation 2.22.4 shows that the aeroelastic divergence problem in Section 2.14 is unchanged; @, 


and 6, are inputs, independent of 0, and @,. Let us suppose that we include a feedback control 
system so that the aileron deflection, é,. responds to a sensor that measures the twist deformation, 
@,. The sensor, perhaps a strain gage, senses and sends signals to an aileron actuator to rotate the 
control surface an amount 6, given by the relationship 


r) 
6,=|G 0| | A (2.22.5) 
2 
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The row matrix |G 0| is called the control gain matrix. The constant G is a “gain" and can be 
chosen by the control system designer arbitrarily; it can be either a positive or negative number. 


The second term on the right hand side of Eqn. 2.22.4 is now written as 


C, c\ Cu 0 0, 
oy “ G 0|\ . 
G,. <) on (itl Na ee 


; Gace b oat olla} 


with 


SI 


(2.22.7) 


C Cc 
E26 ae+(5] "5 (2.22.8) 
CL, e CL 


Eqn. 2.22.7 is a function of @, and @,; it belongs on left hand side of Eqn. 2.22.4. The divergence 
eigenvalue problem for the wing system with this feedback control loop is: 


(S-q) (2) |J@|_ Jo 
(-2-gk) (2-q)||4, 0 (2.22.9) 


The characteristic equation for static stability is again found by taking the determinant of the 
aeroelastic stiffness matrix in Eqn. 2.22.9 and setting it to zero. 


A=q —7q9+6-2gk =0 (2.22.10) 


Figure 2.22.2 plots Eqn. 2.22.10 for five 
different values of k, including the open 
loop case, k=0 (labeled “baseline.”’). 


When k>0 the divergence dynamic pressure 
decreases as the eigenvalue (the crossing 
points) move away from each other. This is 
because when@, >0 the aileron will rotate 
downward and increase lift on the outboard 
panel; as a result, the two segment system is 
destabilized. However, if k<Q when 
6,>0) the aileron will rotate upward to 
"dump" lift from the outboard panel. On 0 2 4 6 8 10 
The two-panel system is stabilized. Notice dynamic pressure parameter 
that, as kK becomes more negative, the value - an : 7 

of the lowest eigenvalue (crossing points) Figure 2.22.2 - Stability determinant vs. q for 5 
increases while the second eigenvalue different values of control constant, k 


aeroetastic stiffness determinant 
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(originally at g = 6) declines so that they approach each other. At some value k <-1 there are no 
real solutions since the determinant will not cross the horizontal axis. The solutions to the 
characteristic equation are complex numbers. 


To take a look at this from a different perspective, we separate the terms in Eqn. 2.22.10 by defining 
the term A, as: 


A,=q -7q¢+6 = (2.22.11) 


The term A,, is the characteristic equation for the "open-loop” system. The divergence condition in 
Egn. 2.22.10 is written in terms of A, as 


A=A,—-2qk=0 or A, =2qk (2.22.12) 


Equation 2.22.10 is again plotted, but this time the range of dynamic pressure values is smaller than 
in Figure 2.22.2. Five lines corresponding to k = 1, 0.5, 0, -0.5 and -lare shown together with the 
function f (¢) = 2qk plotted as a dashed line for the values k=1, k= -1 and k= -1.05. 


Intersection points (circled in Figure 
2.22.3) between the A, curve (labeled 
k=0) and the dashed lines correspond to 
eigenvalue solutions to the closed-loop 
problem and are indicated for the three 
values shown. These circled values 
give the same values for the crossing 
points of the full characteristic equation 
with the controller included. 


aeroelastic stiffness determinant, A 


When k = -1.05, the 2gk line is 
tangent to the A, curve at g, =2.45. 


8 
60 05 10 15 20 .25 
dynamic pressure parameter 


Figure 2.22.3 - Effect of control gain k on divergence 
point at this dynamic pressure where stability determinant crossing points 
the two eigenvalues merge. If k < - 


1.05 then the line 2qk never intersects the baseline determinant and divergence is impossible. 


As a result there is a single intersection 


2.22.2 Example - Active control of typical section divergence speed 


This section examines the typical section wing shown in Figure 2.22.4 with an actively controlled 
aileron used to increase the divergence speed. The controller operates all the time so it not only 
changes the divergence airspeed but it will also change the wing torsional equilibrium position. We 
want to: 1) find the control constant k (defined below) necessary to double the divergence dynamic 
pressure compared to the uncontrolled case; 2) find the value of k for which divergence is totally 
eliminated; and, 3) find how much the aileron must rotate in normal flight. 
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As indicated in Figure 2.22.4 
the lift acting at the 
aerodynamic center is a function 
of @,, @ and 6( the initial 
angle of attack, elastic twist, and 
aileron angle, respectively). A 
aileron hinge restoring moment K,@ is 
(rigid) furnished by the spring. 


Figure 2.22.4 - Typical section with control surface and free body 
diagram The aileron angle 6 is related to 


the elastic twist 9 by the linear 
relationship 6=kO@ where k is 
a control constant that may be 
positive, negative, or zero. The 
effect of the control relationship 
O=k@ is to torsionally stiffen 
(or de-stiffen) the wing to 
increase (or decrease) the 
divergence dynamic pressure. 


Figure 2.22.5 - Free body diagram showing forces and moments Figure 2.22.5 shows the free 


body diagram. The moment equilibrium relation is found by summing moments in the clockwise 


direction: 


+=M_,,=Le+M,.-K,0 (2.22.12) 


pin 


For divergence calculations, @, is zero and the twist @ is a perturbation. The lift and pitching 


oO 


moment about the aerodynamic center are 


Cc im 
L=qSC, |0+—26 (2.22.13) 
a CL, 
M ye =QSCCyc 6 (2.22.14) 
6 


Substituting the relations for lift and pitching moment into Eqn. 2.22.12 with 6 = k@ gives the 
condition for neutral stability. 


C 
K, —qSeC,, —qSeC,, Ga k-—gScCyyc,k |9=0 (2.22.15) 


Ly 


The term contained in the outer brackets of Eqn. 2.22.15 is the aeroelastic torsional stiffness. The 
divergence dynamic pressure is found by setting this stiffness to zero. 
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Kr 
SeC, 
a (2.22.16) 


If k=O, Egn. 2.22.16 reduces to the "open-loop" (control off) divergence dynamic pressure found 
in previous sections. 


Define the parameter q, = “Ys ec, (the value of qp with k=(0). We want to increase the open 
La 


loop divergence dynamic pressure by a factor, n, then 4p = n 
Jo 
in x ! 
x 1+k| <8 4(£) Cures 
C,, a a iy 
2.22.17 a,b 
Ya ( ) 
k= n 
C,, Ce Ls 
-l 
When n=2 k= (2.22.18) 


If the denominator in Eqn. 2.22.16 is less than or equal to zero, then g, <0 and divergence won't 


Cc C 
1+k Gl 14(£) ea (2.22.19) 


occur. 


or 


(2.22.20) 


This is identical to the result we get when we let n — coin Eqn. 2.22.17b. 
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Cc 
When @, #0 the wing lift is L=qSC,_ [« +O+ a 2 ; 


L 


‘a 


The equilibrium twist angle @ is: 


When cies n this equation becomes, after substitution of Eqn. 2.22.17b, 


qo 
qa, 


1-4/ 
n 
SeC, 
where g = tafe . The aileron deflection is: 
7 


Lax 
0=k0= qa, WZ ) 


When n > © we have 


to Galed 


Ores es = k..0 = k..qa, = qa, 


The wing liftis L = gSC,,, (a, +6)+ gSC,,6 so that 


£ 


ale 


Ca /4(2)| Sue 
C,, e)| C 


Cc 


At the divergence free condition n > 0 


L q 


SC, @, -4/ C 
q La 1 vi 4(£) Ls | 


one ea 
gol, a i+(£) Cr 
C 
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Combining the two terms in Eqn. 2.22.23 over a common denominator, 


L=qSC,, a, (2.22.24) 


If the aircraft weighs W pounds and this wing supports half of this weight, then the angle of attack is 


14(£) me 
W c Cuac, 


a, = (2.22.25) 
2gSC aor 
14(£) oe eg 
c Cuac, 
or 
14(£) a's 
c MAC 
a, =~ : (2.22.26) 
2qK_ e C, 
14(£) o 4g 
c Cuac, 


We can now solve for the aileron deflection. Since 6,;,_ free = K-»Q@, we find the expression for 


the aileron deflection required for the actively controlled divergence free wing. 


le 
We C)| Cuac is 


Osi free mre (2.22.27a) 


or 


(Sa 
1% | Ne} Coney 


div—free ~ 
2 FrCuacy 1 e 
+ —_ 
Cc 


(2,22,27b) 


eee, 
Fa 
o) 
SO) 
a le 
nH 
oa sy 
+ 
S| 
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Changing the divergence speed requires a steady, non-zero, value of 6 at all times, not just when 
the wing is perturbed. Equation 2.22.27a is plotted in Figure 2.22.6 as a function of g . From this 
plot we can see that the aileron will be required to deflect to large angle unless its chord is a 
substantial fraction of the chord of the wing or unless the elastic axis offset e is small. 


Figure 2.22.7 plots the results from Eqn. 21.21.27 at three different values of dynamic pressure and 
shows that the best size flap-to-chord ratio will depend upon the dynamic pressure parameter g at 
which the wing operates. 


q 0.0 0.1 0.2 0.3 
flap-to-chord ratio, E 


Figure 2.22.7 — The effect of flap-to-chord ratio on 
divergence suppression requirements. 


Figure 2.22.6 - Aileron deflection requirements for 
divergence suppression 


Control effectiveness - summary 

Unlike divergence, control reversal is not an aeroelastic instability. With divergence a small input, 
@,, produces a large output, @ For reversal a large input 6, produces no output. Above the 
reversal speed, a control input 6, produces lift or moment opposite to that intended (hence the term 


reversal). At pressures larger than the divergence speed, the wing is statically unstable and will 
diverge or twist off. 


Two examples of control reversal were considered. In the first, control effectiveness was defined as 
the ability of a wing/control surface (an elevator or rudder) to create lift which then produced a 
moment about an unspecified flight mechanics axis. In the second example, a wing/aileron 
combination produced rolling motion which was allowed to attain a terminal, steady-state condition. 
The procedure to find the terminal roll rate was illustrated. For this second example, the dynamic 
pressure at the reversal condition was found to be identical to the first case, even though the 
boundary conditions were different. 


In the final example, an aileron was used to increase the divergence dynamic pressure by creating a 
control law relating aileron displacement to wing twist. While not an effective method to produce 
divergence changes, this example does illustrate how feedback control changes aeroelastic 
eigenvalues. In our final section we will illustrate the relationship between static instabilities and 
dynamic response. 
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2.23 The relationship between divergence and dynamic instability-typical section 
torsional vibration with quasi-steady airloads 

In Chapter 4 we will consider flutter, a dynamic instability. To complete the discussion of 

divergence we consider what might happen if a system that diverges is perturbed dynamically 

instead of statically so that it vibrates in the airstream. For this illustration we woll use the simplest 

model, a wing model, shown in side-view in Figure 2.23.1 mounted on a pin is identical to that used 

in Section 2.10 and shown in Figure 2.10.1. 


This wing has span, /, chord, c, and total mass M=ml. The plunge degree of freedom of this typical 
section is zero(h=0Q); twist 6 is measured as a perturbation away from the static equilibrium 
position. The twist angle @ is created by an initial disturbance with initial conditions (twist 
displacement and twist velocity). We assume that the aerodynamic loads on the vibrating system are 
time dependent but identical in form to the static system. This is an approximation and not strictly 
not unless the oscillations occur very slowly. In Chapter 4 we will improve on this formulation. 


Figure 2.23.1 - Typical section model for vibration analysis 


Summing moments about the pin (positive clockwise) and equating them to / 50 where I, is the 


mass moment of inertia computed about the fixed pin, the equation of motion is 
1,0 =-K,0+qSeC,, 6 or I,0+K,0-qSeC, 0=0 (2.23.1) 
With Ky = K; —gSeC,, Eqn. 2.23.1 is written as 


1,0 +K,@ = 0 (2.23.2) 
To solve Egn. 2.23.2 assume that 

At)= QAe~ (2.23.3) 
where s is an unknown, f¢ represents the time parameter and @, is the vibration amplitude. Equation 
2.23.4 becomes 


s'1,4¢° +K,Qe" =(s'1,+ K, ae" =0 (2.23.4) 


Since the vibration amplitude is not zero and e” is not zero, Eqn. 2.23.4 requires that 
(°° Ig+ K. r= 0. This means that the solution for the parameter s is 
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K. K. 
s°=-—" or s=,J-—* (2.23.5) 
I, I, 


Since the term K,. can be positive, negative or zero, there are three different solutions to Eqn. 
2.23.5: (1) if K, >0, implying operation at a dynamic pressure below the divergence g, then 


K. [K. 
S= le =v-l1 7 =+i@; (2) if we are operating above the divergence dynamic pressure 
0 0 


T 


then K; <0 and we have s=+ ; and, 3) when we are operating right at the divergence 


0 


dynamic pressure, Ky =0, s=0. Let’s examine all three cases. 


Case (1) Below the divergence dynamic pressure, we have @ > 0 and the solution 


Kr 


At)= 6,sin a +0,cosa@ (@= oA >0) (2.23.6) 
6 


The constants 6 and @, are amplitude constants determined by initial conditions. The wing, when 


perturbed, vibrates forever, but its amplitude does not change. The natural frequency is determined 
by the aeroelastic stiffness. 


Case (2) The aeroelastic stiffness is less than zero. 


a -K. 
At) = Ge" + Oye P’ (s=p= : r 0) (2.23.7) 
6 


Static analysis says that the displacement is infinite, but dynamic analysis says that the motion has 
two parts. For one part the motion decays but for the other it is exponentially increasing or 
divergent. 


Case (3) K,=0. When K, is exactly equal to zero, Eqn. 2.23.2 becomes Ig 6 = 0 so the 


solution is @= 6,+6,t The motion is neither oscillatory nor divergent. If the system has a small 
displacement it stays there. If it is given an impulse it travels away from equilibrium at constant 
velocity. It is neutrally stable. 
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Figure 2.23.2 shows a plot of nondimensional natural frequency (actual frequency divided by the 


[K. 
reference frequency @, = a vs. airspeed this example. At zero airspeed the natural frequency is 
0 


identical to the so-called in vacuo natural frequency, our reference frequency @, = = : 
0 


Torsional 


torsional natural frequency 
(Hertz) 


0 2 3 
0 50 100 150 200 250 300 350 400 450 500 
airspeed (ft/sec) 


Figure 2,23.2 - Nondimensional natural frequency vs. airspeed 


Divergence is a special case of dynamic instability. At the so-called statically unstable point, 
perturbed dynamic motion changes from simple harmonic oscillation to exponentially divergent 
motion. At this transition point the natural frequency is zero. Divergence has been called flutter at 
zero frequency. However, this term is not really correct since flutter, as we will show in Chapter 4, 
has a different instability mechanism. 


2.24 Summary - lift effectiveness, divergence and aileron effectiveness 

After over 100 pages of text and examples (with pages and pages of detailed algebra and equations), 
what should we take away from this chapter? First of all, if we had this material in 1920, when 
aviation was in its infancy, we would have been very popular with airplane designers. Our simple 
models could be used to explain the source of observed problems that had begun to be important and 
suggested ways to fix them. Despite the tremendous progress made in theoretical modeling over the 
last half century, this ability to use fundamental models to understand, avoid, diagnose and outline 
fixes to aeroelastic problems remains. 


The aeroelastic problems we considered in this chapter have been exhaustingly discussed by others. 
Several interesting articles cover the early years of structural development and aeroelastic problems 
and are both informative and entertaining.*” ** ** ** Two references, one by Fung (Reference 25), 
the other by Bisplinghoff, Ashley and Halfman (Reference 26) stand out, both for their clear 
discussions and the wealth of back-up references. 


Torsional stiffness is fundamental to all aeroelastic problems 


The fundamental importance of wing torsional stiffness to the sources and solutions to aeroelastic 
problems was confirmed in all sections of Chapter 2. The interactions between wing torsional 
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stiffness and the aerodynamic loads produce an aerodynamic wing de-stiffening effect that works its 
way into phenomena such as load amplification, ultimately culminating in divergence, and control 
effectiveness, leading to aileron reversal. 


We also introduced and defined terminology with definitions that included terms such as elastic axis, 
shear center, aerodynamic center, stiffness matrix and aerodynamic stiffness matrix. These terms 
helped us to characterize the important features of configurations that allow us to make preliminary 
assessments of problem areas. 


Boundary conditions are important 


In the area of analysis we also showed the importance of boundary conditions, restrictions and 
assumptions. For instance, when the wing angle of attack is held fixed, as it would be in a wind 
tunnel experiment, twist deformation increases without bound as airspeed increases. This 
unbounded distortion is caused by the approach of a static instability, divergence. However, when 
the angle is adjusted to keep the lift constant, as it would be in flight, this unboundedness is not 
present. 


These analytical features led us to consider what happens when we introduce small perturbations 
into the static equilibrium analysis. In this case we developed an eigenvalue problem from which we 
could solve for the divergence speed. We also found that the eigenvalue problem was simply the 
static equilibrium problem, but with all inputs set to zero and with the interpretation that deflections 
were perturbations rather than real deformations. We also introduced a simple vibration problem 
and showed that there is a relationship between the more realistic dynamic analysis and the static 
analysis. At divergence, one of the system natural frequencies goes to zero and the dynamic 
response is exponentially divergent in time. 


Divergence and control effectiveness are different but have the same origin 


Wing divergence happens when a small input, @&,, produces a large output, @ Control reversal is a 


performance problem related to wing torsional stiffness. With aileron reversal a large control input 
0, produces no output. Two control reversal situations were examined. In the first case, control 


effectiveness was defined as the ability of a wing/control surface (an elevator or rudder) to create lift 
which then produced a moment about an unspecified flight mechanics axis. In the second example, 
a wing/aileron combination produced rolling motion which was allowed to attain a terminal, steady- 
state condition. The value of the terminal roll rate was found. Reversal was defined as the inability 
to create a roll rate, not a roll moment. In both cases the reversal condition was found to be 
identical, even though the goal was different. 


As indicated in Figure 2.24.1 the offset between the wing section aerodynamic center and the shear 
center (or elastic axis) is a critical parameter for divergence. The shear center (elastic axis) depends 
on structural geometry, primarily the placement and size of major elements like wing spars. For 
aileron effectiveness the shear center is not important. Instead, aileron size defines the aileron center 
of pressure. The difference between the center of pressure and the aero center, shown as the 
dimension d, defines reversal. In both cases the torsional stiffness is important. 
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aero center 
Shear 
Divergence 
Aileron center of 
pressure 
Aileron 
reversal 


d 


Figure 2.24.1 — Two different geometric parameter define divergence and lift effectiveness. 


In Chapter 3 we will complete our study of the fundamentals of static aeroelasticity by considering 
aeroelastic effects on swept wings. For swept wings, aeroelastic interaction is more complicated 
because bending introduces additional loads. 
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2.25 Homework problems 


2.1 An uncambered wing is mounted in a wind tunnel on a 

shaft that is free to rotate. The mount is restrained by the aan a Spring Kr ee 
usual torsion spring Ky. In addition, a pendulum is attached to 

the wing mounting shaft outside the wind tunnel. If the wing 

c.g. is located at the pivot point, solve for the divergence g d 
dynamic pressure. | 


2.2 An all-movable, uncambered wing segment with planform etapa 


area S and lift curve slope Chey is used to control an air P) 


vehicle. The idealized model of this system is shown. It a 
consists of a rigid rod that can rotate about a forward pivot, 

but is restrained by a torsion spring, K;. This rotation is 

measured as the small angle @ shown. No aerodynamic 

forces act on this rod. 


~+_oa 


At the right end of the rod is an wing segment that is “all- antes 


movable” (as a result, CLs = Cr and Cyac = 0) about a Problem 2.1 

pivot located a distance b from the left end of the rigid rod. This wing is mounted on the rod with 
another torsion spring K2 that resists wing torsion, @ The torsion spring Ky, is unstretched when 
a = 0 and the torsion spring K, is unstretched when 0 =0. The wing segment is given a rotation 
input angle 6,. The total angle of attack of the 


torsion spring 


wing is then -& + 0 + 6. torsion spring 


2.2 Problem statement 

(a) Break the deformed system into two pieces. 
Draw the two free body diagrams of the deformed 
system. Derive the coupled static equilibrium 
equations and write them in matrix form as 


Ki, Kil _ JQ — Sey, 
oa ee =q0, where gq =———~* 
Ky, Ky |(@ Q, Ky 


(b) Solve for the divergence dynamic pressure at which this system is neutrally stable. Input the 


b K, 

ratio —=11 into this expression and find the value of — for which this configuration will not 
e 2 

diverge. 
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2.3 A rigid wing is attached to a wind- 
tunnel wall by a linear torsion spring. 
This wing also has a small wing attached 
to its tip. The idealized tip surface 
produces lift LL, according to the 


assumed simple relationship 


Ls = DB incr os as 


where Q@, is the streamwise angle of 


attack of this small surface; C Las 18 the 


tip surface lift curve slope, while S tip 1S 


the tip surface area. 


The tip device is connected to the wing 
tip by a torsion spring with stiffness k 
in-lb/radian. When the tip device rotates 
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torsional spring Bi 
) K. 


in shear center 


gap accentuated 
for illustration 


view A-A V 
enlarged view of 
cross-sectional 

geometry 


Problem 2.3 - Wing with flexible tip 


an angle B with respect to the wing tip a restoring torque kf is generated to oppose this rotation. 


The idealized wing liftis L=qSC,,@ with @=@,+ 0. 


2.3 Problem statement 


a) Write the equations of torsional equilibrium for this system when the entire system is given 
an initial angle of attack. Express these equations in terms of the following non-dimensional 


gSeC, Sipe 


parameters; g = > Spe 
T 


k : we : 
a . Derive the characteristic equation for 


T 


aeroelastic divergence g,, in terms of these parameters. 


b) Plot the divergence dynamic pressure parameter g, as a function of Sp when k, =10 and 


kp =l 


Problem 2.4 

Two identical uncambered 
wing segments are connected 
to each other by a torsional 
spring and to the wind tunnel 
walls by two other torsion 
springs, as indicated in the 
figure. The three torsion 
springs have the same spring 
constant Kr. The two 
segments are mounted on 


bearings on a spindle attached 
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to the tunnel walls. Solve for the divergence dynamic pressure. Find the mode shape at 
divergence. Describe this mode shape (how do the surfaces move relative to each other?). 


Problem 2.5 

Two wing sections are mounted on 
shafts and attached to each other and 
to a wind tunnel wall, as indicated. 
Note that the torsion spring 
stiffnesses are equal, but are offset an 
amount d. This configuration is 
similar to, but not identical to, the 
example in Section 2.18. When the 
airfoils are placed at an angle of 
attack Oto, the springs deform. Lift on 


V 
line of aerodynamic centers | 


\ X\ torsional spring 


torsional spring 


stiffness K (in-lb/tad) Stiffness K (in-Ib/rad) 


/92t Oo 


the two, identical, uncambered wing - eels 
sections is snacagMtedda 
A 


view A- 


L, =4qSC, AG, +6, ) Problem 2.5 


L, =qSC,.(@, + 4) 


Derive the matrix equations of torsional static equilibrium for this model when it is placed at 
Q, 


} Do not solve for 
2 


— 16 
an angle of attack ,. Express these equations in the form [K ‘ | 5, -| 
2 


the deflections, but identify the structural stiffness matrix and the aerodynamic stiffness matrix 
in these equations (“Identify” means put some words together and attach an arrow.) 


Solve for the divergence dynamic pressure; solve for the value of d that eliminates divergence. 


The shaft attaching the two segments can be re-positioned, but the outer wing segment always 


Place an aileron on the outer (right-hand) section. Develop the static equilibrium equations for 
the system when there is no initial angle of attack, but the control surface is deflected 


(a) 
(b) 
remains in-line with the inner segment. 
(c) 
downward an amount 56. 
(d) 


Solve for the rolling moment generated by the aileron. Solve the reversal speed in terms of 
general parameters such as C ” and C MAC, * 
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Problem 2.6 

A side view of two idealized wings is 
shown in the figure. This example is 
similar to that discussed in Section 
2.15, but has one spring removed. 
These wings are connected by a 
single spring, with spring constant k 
Ib/inch; each wing has planform area 
S. The spring is pinned to each wing 
and develops an internal force in 
response to relative deflection 
between its ends. This configuration 
differs from that discussed in Section 
2.15 because it can rotate without 


stretching or compressing the spring. 


This is called rigid body freedom. 
Problem 2.7 - Tandem wing geometry with single spring 


The lift on each wing is given by 


L, =qSC,, (9, + &) 
L, =4SC,4(8, + &,) 


where @, is an initial angle of attack common to both wings and @, and @, are the two twist angles. 


a) Solve for the characteristic equation A for static stability 

b) Ife; = e.=e, identify a dynamic pressure parameter g and plot the characteristic equation as 
a function of this parameter. 

c) Find the value(s) of the pressure 
parameter at which A= 0. Solve 
for the mode shape(s). Will the 


system diverge? Why or why not? metallic layer 
piezoelectric sub-layer 


Problem 2.7 


A wing test article consists of a low density, symmetrical ae z 
stiffness, wrapped around the two layer piezoelectric plate 2 
creates an electric field E shown in Figure (b) that causes € Problem 2.7(a) - Piezoelectric layered 1 
surface is formed. actuator 

3 2 


The wing idealization is shown in Figure (c); for 
analysis purposes we have only drawn the chord electrodes 
line. This model has two lift components; the first 
is due to an initial angle of attack and the wing 
twist. The second is due to the applied voltage that 


creates camber. 
voltage 


Problem 2.7(b) - Piezoelectric actuator 
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L, = qSC,, (@, + 0) 


ae 


f 
ey >< 
f \ 


We will represent the added lift and pitching 
moment combination created by camber as a 
single force L, acting on the wing at a 


distance e, aft of the shear center (this is the 
center of pressure). This location of the center i a eee 
of pressure for the additional force and we V 
assume that it does not move with applied torsion spring Kr 
voltage. 


The relationship between lift L, and the 


Problem 2.7(c) - Active airfoil with applied loads 
applied voltage V is: 


L, = (qSC,, )V 
V is voltage, not airspeed. Since the change in camber produces lift and pitching moment, the 
oC Oo camber 
aerodynamic coefficient C, = ——_4+— ——— tells us how much lift we get per volt across 


VY @camber OV 


the active material per unit wing area per unit dynamic pressure. 


Problem statement 
(a) When the airfoil initial angle of attack is zero (@, =0) solve for the lift on the wing as a 


function of the applied voltage V . 


(b) Find the wing divergence dynamic pressure gp as a function of the aerodynamic and 
structural parameters. (Voltage should not appear). 


(c) Let’s create a feedback system where we measure the twist @ and feed this back a to the 
voltage, so that V=k@ (k is an arbitrary number that we can choose). In this case, the voltage 
term should move to the left side of your equation. If the wing is placed at an initial angle of attack 


Q, , solve for the wing twist and solve for the divergence dynamic pressure. 


Solution hints 
The total airfoil lift is written as 


L=L +L, =qSC,,a+qSC,,V 


where a=a,+0 
In this problem @, =O. When V=k@, the equation of torsional static equilibrium is 
gSe,C,_ (a, + @)— qSesC,,, kO- K,0 =0 


Problem 2.8 
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The 1 dof wing shown in Figure (a) has leading edge and trailing edge control surfaces. There is no 
initial angle of attack when the two control surfaces are undeflected. Three lift components act on 
the idealized wing. Each force has its 


own location with respect to the shear 2 

center pin shown in the figure. Twisting : C 0 

is resisted by the torsional spring with V s 
torsional stiffness Kr. 0.255 4 


The entire assembly is mounted at the end 
of a long bar a distance r from the center shear center/pin 
of rotation. This assembly will move 

upward at a terminal speed v = pr when 

the surfaces are deflected. The leading 

edge surface and the trailing edge flap are 

geared together so that the three lift components are 
as follows: 


Problem 2.8(a) 


L, =qSC,,9+ qSC,, (-Y4,) (where v = pr) 


Ly = 4SC1,509 


Ly = 0.25gSC,,56 


I 
Problem 2.8(b) — Top view of wing assembly 


The rolling moment is M,,, = Lr = (L; + Ly + L3)r. showing yotbiondl fates pi 


Problem statement 
a) Solve for the twist 0. 
b) Find the expression for the steady state roll rate p. 
c) Find the reversal dynamic pressure 


Problem 2.9 

An uncambered wing with span, b, is attached to both walls of a wind tunnel and placed at 
an angle of attack @. The wing is idealized as a torsionally flexible, uniform property 
element like that discussed in Section 2.19. Use Figure 2.19.1 as your reference, but with 
span b instead of L and with an additional support on the right. 


The differential equation of torsional equilibrium used to model wing twist deformation is: 


d’@ . qcea, 9 = fee 
dx GJ GJ 


(a) Find the divergence dynamic pressure of this wing. 
(b) Torsion springs with stiffness, k in-lb./radian is placed at the left and right ends of 
the test article (at x=0). Solve for the expression characteristic equation for 
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divergence in terms of a parameter B = kb/GJ. Plot the divergence dynamic 
pressure as a function of B to generate a figure similar to Figure 2.19.9. 


Problem 2.10 

Classical wing divergence ignores the fact that the wing is attached to the fuselage with pitch and 
plunge freedom. Consider the idealized, uncambered wing attached to a torsion spring which is, in 
turn, attached to an idealized fuselage, as shown in Figure 1. The torsional spring resists torsion by 
developing a restoring moment M, equal to K,@, where @ is the relative rotation of the wing with 
respect to the fuselage. 


The fuselage is a freely flying article, but for our model it is attached to a frictionless pin to simulate 
the ability of the aircraft to rotate about its own center of gravity (cg). The wing has two degrees of 
rotational freedom @ and @, that generate lift. 


An all-movable tail on the fuselage provides a pitching moment about the aircraft c.g. by 
generating lift, L,,,,, given by the expression 


Lait = Praite (a +B) (1) 


L Gail 


In Eqn. 1, Z is the tail rotation with respect to the fuselage while @ represents the rotation 
of the fuselage about the airplane cg or model pin. 


wing elastic axis 


a offset, e tail lift 
torsional spring S L 
Pe Ky 


note: airplane cg is 
usually ahead of 
wing aerodynamic 
center 


Figure 2.10.1 - Wing fuselage configuration geometry 


Problem statement 
(a) Solve for expressions for @ and @ in response to the tail angle, # using the following 
definitions 
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: (2) 


As noted on Figure 2.10.1, the pin representing the c.g. is usually ahead of the wing 
aerodynamic center. If this is so, then the dimension d will be negative. When the 
parameter R defined in Egn. 2 is multiplied by -1 the result is similar to the definition of 
the tail volume used in stability and control studies. If R>1 (or —R<-—l) then the 
configuration is stable in pitch in forward flight. 


(b) The purpose of the tail control surface is to create a change in pitch angle. Define control 
reversal as the condition that occurs when we deflect the tail surface, but there is no change in 
aircraft attitude angle. Set up an analysis to determine whether or not tail control reversal g 
occurs. Solve for the tail control reversal dynamic pressure. 


qSeC, 


(c) Solve for divergence dynamic pressure parameter g = = as a function of R. Plot 


is 
this divergence pressure parameter as a function of R. What value of R is required to 
make the divergence pressure parameter equal to 1? 


Helpful hints - The FBDs are set up in Figure 2.10.2. They must be useful. My equilibrium 
equations are also shown below. So, too, is the characteristic equation. 


FBD for pivot (cg) point 


wing shear 
Lw center FBD 


Ms =K,0 


Rw= Lw 


Figure 2.10.2 - Wing/fuselage free body diagram 


=. 4 LY Sau oe = -7{£}( Sa Charait 
it Z 5 |S ree As Ge. 


and ga—[l—g]@=0 so that A =—g[—R)-g)+ q]=—-q(— R+ R@) 
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